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Chapter 

Abstract 



In this Thesis we have investigated some effects appearing in top quark observables, in the framework of 
the Minimal Supersymmetric Standard Model (MSSM). 

The Standard Model (SM) of Strong (QCD) and Electroweak (EW) interactions has had a great 
success in describing the nature at the Electroweak scale, and its validity has been tested up the the 
quantum level in past and present accelerators, such as the LEP at CERN or the Tevatron at Fermilab. 
The last great success of the SM was the discovery in 1994 of its last matter building block, namely the 
top quark, with a mass of rrit — 173.8 ± 3.2 ± 3.9 GeV/c^. However the mechanism by which all the SM 
particles get their mass is still unconfirmed, since no Higgs scalar has been found yet. The fermions couple 
to Higgs particles with a coupling proportional to its mass, so one expects that the large interactions 
between top and Higgs particles give rise to large quantum effects. 

We have focused our work in the MSSM. This is an extension of the SM that incorporates Super- 
symmetry (SUSY). Supersymmetry is an additional transformation that can be added in the action of 
Quantum Field Theory, leaving this action unchanged. The main phenomenological consequence of it is 
that to any SM particle {p) there should exist a partner of it, which we call sparticle {p). This extension of 
the SM provides elegant solutions to some theoretical problems of the SM, such as the hierarchy problem. 

We have computed the radiative corrections to some top quark observables, using the on-shell renor- 
malization scheme, and with a physically motivated definition of the tan/3 parameter, tan/3 is the main 
parameter of the MSSM, and it governs the strength of the couplings between the Higgs bosons and the 
fermion fields (and its superpartners) . 

We have computed the SUSY-EW corrections to the non-standard top decay partial width into a 
charged Higgs particle and a bottom quark r(t ^ H^b). We have found that these corrections are large 
in the moderate and specially in the high regime of tan/3, where they can easily reach values of Sswit ~^ 
H^ b) ~ -1-30% for negative /i and a "light" sparticle spectrum, and Sswit -^ H'^ b) ~ -1-20% for positive 
/i and heavy sparticle spectrum. In both cases we have singled out the domain jiAt < of the parameter 
space, which is the one preferred by the experimental data on radiative B-meson decays (6 — > S7). We 
have singled out the leading contribution to this corrections, which is the supersymmetric contribution 
to the bottom quark mass renormalization constant 6mi,/mh. It is proportional to —jiAt and shows a 
possible non-decoupling effect. The contributions from Higgs particles is tiny, and can be neglected. We 
have added this corrections to the known Strong corrections {5qcd — —60%, 5susy-qcd — +80% and 
dsusY-QCD — —40%) and we have look at its impact on the interpretation of the Tevatron data. The 
standard analysis (using only Sqcd) implies that for a charged Higgs mass of M[j± =110 GeV the values 
of tan/3 > 50 are excluded. If this charged Higgs boson belongs to the MSSM the excluded values are 
tan/3 > 35 or tan/3 > 75 for the two scenarios presented above. So no model independent bound on the 
charged Higgs mass can be put from present experimental data. 

We have looked at the possibility that the top quark could decay via Flavour Changing Neutral Current 
(FCNC) into a neutral Higgs particle and a charm quark. We have computed the EW contributions 
and the QCD contributions, using a mass model motivated by Grand Unification Theories (GUT), but 
not restricted to any specific GUT. We have included the full interaction lagrangian between all the 
particles. The upper theoretical bounds are found to be BR^^^^^^^{t —^ ch) < several x 10^^ and 
^^SUSY-QCD(^^ -> eft,) < several x lO^"* and the typical values for this ratio are 10"* and 10"^ - 10"'' 
for the SUSY-EW and SUSY-QCD induced FCNC decays respectively The Higgs and the purely SUSY 
contributions to the EW induced process are of the same size, and can be of the same or opposite sign. 
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We have found that the SUSY-QCD induced FCNC decay widths are at least two orders of magnitude 
larger than the SUSY-EW ones in most of the parameter space, thus making unnecessary the computation 
of the interference terms. The value of this branching ratio is too small to be measured either at the 
Tevatron or at the Next Linear Collider (LC)^, but there is chance that it could be measured at the Large 
Hadron CoUider (LHC). 

If bottom-like squarks (the superpartners of bottom quarks) are heavy enough they could decay into 
a top quark and a chargino (the superpartners of gauge bosons and Higgs bosons): b —^ tx~ ■ This 
could serve as an unexpected source of top quarks at the Tevatron, at the LHC, or at the LC. We have 
computed the QCD radiative corrections and the leading EW corrections to this partial decay width. The 
QCD corrections are dominant, they are negative in most of the parameter space, and are of the order 
of 5QCD{bi -^ txi) — —60%, SQCD{b2 -^ tXi) — —20% for a wide range of the parameter space. EW 
corrections can be of both signs. These corrections have been computed in the higgsino approximation, 
which gives the leading behaviour of the EW corrections. Our renormalization prescription forces the 
physical region to a narrow range. Within this restricted region the typical corrections vary in the range 
^Ewibi -^ txi) — +25% to — 15% SEw{b2 -^ tXi) — +5% to —5%. However we must recall that these 
limits are qualitative. In the edge of such regions we find the largest EW contributions. We stress that 
in this case it is not possible to narrow down the bulk of the corrections to just the renormalization of 
the bottom quark Yukawa coupling. 

Our general conclusion is that the supersymmetric strong and electroweak radiative corrections can be 
very important in the top/bottom-Higgs super-sector of the MSSM. Therefore, it is necessary to account 
for these corrections in the theoretical computation of the high energy physics observables, otherwise 
highly significant information on the potentially underlying SUSY dynamics could be missed. This is 
true, not only for the future experiments at the LHC and the LC, but also for the present Run I data 
(and the Run II data around the corner) at the Fermilab Tevatron collider. 



^Note Added: See however note on section 5.5 (pg.|76|) 



Chapter 1 

Introduction 



Recently, the Standard Model (SM) of the strong (QCD) and electroweak (EW) interactions 0-§l has 
been crowned with the discovery of the penultimate building block of its theoretical structure: the top 
quark, i y, 1^. At present the best determination of the top-quark mass at the Tevatron reads as 
follows [g|: 

mt = 173.8 ±3.2 (stat.) ±3.9 (sist.) GeV. (1.1) 

While the SM has been a most successful framework to describe the phenomenology of the strong and 
electroweak interactions for the last thirty years, the top quark itself stood, at a purely theoretical level 
-namely, on the grounds of requiring internal consistency, such as gauge invariance and renormalizability- 
as a firm prediction of the SM since the very confirmation of the existence of the bottom quark and the 
measurement of its weak isospin quantum numbers M . With the finding of the top quark, the matter 
content of the SM has been fully accounted for by experiment. Still, the last building block of the SM, viz. 
the fundamental Higgs scalar, has not been found yet, which means that in spite of the great significance 
of the top quark discovery the theoretical mechanism by which all particles acquire their masses in the 
SM remains experimentally unconfirmed. Thus, it is not clear at present whether the SM will remain as 
the last word in the phenomenology of the strong and electroweak interactions around the Fermi's scale 
or whether it will be eventually subsumed within a larger and more fundamental theory. The search 
for physics beyond the SM, therefore, far from been accomplished, must continue with redoubled efforts. 
Fortunately, the peculiar nature of the top quark (in particular its large mass-in fact, perhaps the heaviest 
particle in the SM!- and its characteristic interactions with the scalar particles) may help decisively to 
unearth any vestige of physics beyond the SM. 

We envisage at least four wide avenues of interesting new physics potentially conveyed by top quark 
dynamics and which could offer us the clue to solving the nature of the spontaneous symmetry breaking 
(SSB) mechanism, to wit: 

1. The "Top Mode" realization(s) of the SSB mechanism, i.e. SSB without fundamental Higgs scalars, 
but rather through the existence of tt condensates |^ ; 

2. The extended Technicolour Models; also without Higgs particles, and giving rise to residual non- 
oblique interactions of the top quark with the weak gauge bosons l9l ; 

3. The non-linear (chiral Lagrangian) reahzation of the SU{2)l x U{1)y gauge symmetry [[l0|-[l2|, 
which may either accommodate or dispense with the Higgs scalars. It can also generate additional 
(i.e. non-standard) non-oblique interactions of the top quark with the weak gauge bosons [[l3|,[l4|; 
and 

4. The supersymmetric (SUSY) realization of the SM, such as the Minimal Supersymmetric Standard 
Model (MSSM) ||l5Hl8|| (see also ||l^ for a comprehensive review), where also a lot of potential new 
phenomenology spurred by top and Higgs physics might be creeping in here and there. Hints of this 
new phenomenology may show up either in the form of direct or virtual effects from supersymmetric 
Higgs particles or from the "sparticles" themselves (i.e. the i?-odd jT^-Q partners of the SM 
particles). 

In this Thesis, we shall focus our attention on the fourth large avenue of hypothetical physics beyond 
the SM, namely on the (minimal) SUSY extension of the SM, the MSSM, which is at present the most 
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predictive framework for physics beyond the SM and, in contradistinction to all other approaches, it has 
the virtue of being a fully- fledged Quantum Field Theory. Most important, on the experimental side the 
global fit analyses to all indirect precision data within the MSSM are comparable to those in the SM; in 
particular, the MSSM analysis implies that mt — 172 ± 5 GeV |2y,|2^], a result which is compatible with 
the above mentioned experimental determinations of mt- 

Moreover the SUSY theories are a step forward in the search for an unified theory. On the "light 
energy" point of view a simple supersymmetrized version of the SM yields to a unification of the three 
gauge couplings of the SM, at a scale of A ~ 10-^^ GeV |2|] . On the "high energy" point of view these 
theories can be embedded in a more general framework, Superstring Theories, these theories provide 
a unification of "Classical" Quantum Field Theory with Einstein's Gravitational Theory. Some SUSY 
models of unification have been constructed, that provide unification at large scales, the most hard 
constrain that they must fulfill is to reproduce the SM at the EW scale within the present experimental 
constrains. Maybe the most popular of such model is the so called "Minimal Supergravity" (mSUGRA) 

If i?-parity is conserved SUSY (i.e. i?-odd) particles cannot decay into SM ones (see section |2.1| ), 
thus the Lightest Supersymmetric Particle (LSP) is stable. The LSP is a good candidate for cold dark 
matter, which is a necessity of the present cosmological models. Cold dark matter is necessary to explain 
the flatness of the present universe, and the structure formation mechanism. 

SUSY theories also provide answers to some SM questions, they provide a natural solution to the 
"hierarchy problem" |g^, that is, the impossibility, in the SM, of having two scales (the EW scale, and 
the uniflcation scale) with a large gap between them. This is due to the presence of quadratic divergences 
in the one-loop correction of the boson masses. These divergences appear because of bosonic loops in the 
mass correction. In SUSY theories each bosonic (fermionic) particle has associated a fermionic (bosonic) 
partner, with the same quantum numbers and couplings, thus the fermionic loops cancel the quadratic 
divergence of the bosonic loops, and the two scales remain stable. 

Aside from these facts SUSY theories can be useful also in other subjects, for example they give us 
hints about quark confinement ||2J]. The excitement is so great that a sole event at the Tevatron, not 
expected in the SM, has produced a full analysis of its expectation as a SUSY event |2^ ]. 

All these in a hand, it seems that SUSY could be the solution of all our theoretical problems (or our 
theoretical prejudices) in particle physics, however no supersymmetric particles have been found in the 
high energy physics experiments yet, or, to put it in other words, only "half" of the MSSM spectrum has 
been found (aside from the Higgs sector). Thus SUSY can not be an exact symmetry of nature at the EW 
scale, and we would seem forced to abandon this nice framework. However there exist a mechanism of 
breaking SUSY without losing its most important properties, it is called "Soft-SUSY-Breaking" [g6[. At 
scales lower than the Soft-SUSY-Breaking scale the model can be described by a set of parameters which 
determine the spectrum of the SUSY partners of known particles. One thinks that when the masses of 
the superpartners are large enough the supersymmetric particles eventually decouple pj] , though it has 
not been demonstrated yet. 

Soft-SUSY-Breaking can be realized by means of different mechanisms. Each of these mechanisms 
provides us with a different set of Soft-SUSY-Breaking parameters at the EW scale, determined by a 
small set of parameters at high energies. 

In this Thesis we will take the point of view that the MSSM is the effective theory at the EW scale 
of a more fundamental theory, which we do not know about, thus we will treat the Soft-SUSY-Breaking 
parameters as being arbitrary, within the allowed experimental range. 

Radiative corrections [ES,M have shown to be a powerful tool in particle physics for the last half 
century. Recall the first theoretical and experimental determination of the electron anomalous coupling 
(g — 2) pn] as one of its earliers applications, and the measurement of the radiative corrections to precision 
EW observables (such as the relation Mz/M]^) at present high energy colliders as the most recent one (see 
e.g. [pl[). Radiative corrections are useful also to determine (indirectly) the existence, and the parameters, 
of particles yet to be discovered. As a matter of fact the mass of the top quark was estimated, before its 
direct observation at the Tevatron [H, p| , with the help of its radiative corrections to the Mw — sin^ 9w 
correlation |32]. One could think of estimating also the Higgs mass by this method, unfortunately the 
one-loop Higgs radiative corrections enter this observable as the logarithm of its mass pJ , whereas the 
effect of the top quark grows quadratically with its mass. 

It is a wonderful idea, spread all over the theoretical particle physicist community, the use of radiative 
corrections to determine if there is any physics beyond the SM. One can look at the precision observables, 
taken out of present high energy colliders, and search for deviations of the SM. We must note that present 
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precision data does not present significant deviations from the SM expectations, but this has not been 
always the case. Some time ago there was a a large discrepancy between the SM prediction and the 
experimental measurement of the hadronic fraction of Z decays into bb pairs. This discrepancy could 
be cured by introducing in the theoretical estimate the SUSY radiative corrections |p4,pq| (see also a 



complete study of the Z boson in the MSSM in |p6,37|). Nowadays this mismatch has been brought 
down to non-significant deviation (less than 1 standard deviation), however we learned that using these 
precision measurements a precise prediction on the MSSM spectrum could be found trough global fits to 
electroweak precision data [pOl, pTI, IsM] . 



In this Thesis we will address the important issue of the EW SUSY effects on top quark and Higgs 



boson physics. The top quark presents a privileged laboratory for EW physics, due to its large mass (1.1), 
as the Higgs particle couples to fermions proportionally to its mass. In the case of SUSY theories this 
privilege is enhanced for several reasons. First of all the Higgs sector is extended into a so called "Type 
II Two-Higgs-Doublet Model" (2HDM) p9|, and the Yukawa couplings of the top-bottom weak doublet 
become (normalized with respect to the SU{2) gauge coupling) 

. _ ht mt , _ hb rub 

g V2Mw sin/3 g V2Mw cos/3 

where tan /3 is the ratio between the vacuum expectation values of the two neutral scalar Higgs bosons 
(see chapter 0). Notice that in this extension of t h SM it is not only the top quark that can have 



large Yukawa interactions with Higgs bosons. From (1.2) we see that at large tan/3 ( > 30) the bottom 
quark Yukawa coupling also becomes important. Second, the presence of the superpartners of the top 
and bottom quarks ( "stop" and "sbottom" ) and those of the Higgs bosons ( "higgsinos" ) raise up a very 
interesting top-stop-Higgs-higgsino phenomenology. 

The SUSY radiative corrections to the top quark standard decay mode into a charged gauge boson and 
a bottom quark have been known since some time ago |^,^ (see also [Q for an exhaustive analysis). 
Also the conventional strong (QCD) corrections regarding the phenomenology of top and the charged 
Higgs are well known [|3-47|, and its strong SUSY radiative corrections have been studied too ^,|4F 



Thus the following step is to determine the importance of the Yukawa couplings to the top-Higgs sector 
phenomenology pOf . 

The aim of this Thesis is to study the effects of the radiative corrections to the top-Higgs sector in 
the MSSM, by looking at unconventional decay and production modes. We will show that EW radiative 
corrections are important, and this has an effect both in the interpretation of the present experimental 
data (Tevatron Run I) | pl[ and on the prospects of measurements in future colliders (Tevatron Run II, 
Large Hadron Collider -LHC-, and next Super Linear Collider -LC-) [|2J. 

Moreover one expects that, if SUSY particle exists, they could be an unexpected source of top quarks 
at high energy colliders. The observed top quark production cross section at the Tevatron is equal to the 
Drell-Yan production cross-section convoluted over the parton distribution times the squared branching 
ratio. Schematically 
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Xpp^tt) = / dqdq a{qq-^tt) x \BR{t^ W+ b)\'^ . (1.3) 



However, in the framework of the MSSM, we rather expect a generalization of this formula in the following 
way: 

o-obs. = dqdq a{qq^tt) x\BR{t^W+b)\^ 

+ I dqdqa{qq^gg) x \B R{g ^ t~ii)\'^ x \BR{t -^ W+ b)\'^ 

+ fdqdqaiqq-^baba)x\BRiba~^tXi)fx\BR{t^W+b)f + ..., (1.4) 

where g stand for the gluinos, ti for the lightest stop and ba{a ~ 1,2) for the sbottom quarks. One 
should also include electroweak and QCD radiative corrections to all these production cross-sections 
within the MSSM. For some of these processes calculations already exist in the literature showing that 
one-loop effects can be important on sparticle production [p3|-p6[ as well as on sparticle decays, both the 
QCD [|3,|§ and the EW |5^ MSSM corrections. 
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Thought we have been mainly interested in a scenario where the charged Higgs particle is lighter than 
the top quark, an obvious question is what would happen if this charged Higgs is heavier than the top. We 



have considered this issue in Ref. |60| (see also an exhaustive analysis in ||6l|). The radiative corrections 
in the top- Higgs sector in the MSSM should be compared with those from the generic 2HDM's. We have 
been interested in these extensions of the SM in Ref. p3| and more work is currently in progress. The main 
result is that if a charged Higgs boson is found, one could discriminate to what kind of model it belongs 
by using radiative corrections. These calculations of radiative corrections are in the line of completing 
(within the same order of perturbation theory) our previous studies of the full set of three-body decays 
of the top quark in the MSSM g]. 

The structure of this Thesis is as follows: in chapter we give the basic notation of the MSSM used 
throughout this Thesis; in chapte r p| we explore the renormalization of the MSSM, extending the well 
known formalism used in the SM ]28|,|29|], and using a physically motivated renormalization prescription; 
chapters |j to ^ deal with explicit effects of the one-loop corrections on some physical processes of top 
quark production and decay; and finally in chapter we present the general conclusions. At the end we 
include an appendix with some technical details and notation. 



Chapter 2 

The Minimal Supersymmetric 
Standard Model (MSSM) 

2.1 Introduction 

It goes beyond the scope of this Thesis to study the formal theory of Supersymmetry Q,^, however 
we would like, at least, to give a feeling on what is it. 

Supersymmetry (SUSY) can be introduced in many manners, maybe the most straightforward one 
is adding to the space-time coordinates (t, x) another set of coordinates 6a (a = 1, . . . , n the space-time 
dimension) that are Grassmann variables, i.e. they anti-commute. Now the general "rotations" in this 
space are a superset of the Poincare transformations of space-time. It is clear that being 9a Grassmann 
variables the generators of the rotations that involve these coordinates will behave in a special way, and 
indeed they do. These generators (usually called Qa) anti-commute with themselves, so they do not form 
an Algebra, but a Super-Algebra, and the SUSY transformations do not form a Lie Group. However 
it turns out that it is the only external transformation that can be added to the Poincare Group, and 
leave the Scattering (5) matrix untransformed. One can add as many "supersymmetries" (i.e. sets of 
9 variables) as the dimension of the space-time, thus if we introduce a single set of 9 it is said that we 
have a iV = 1 supersymmetry, and so on. The structure of the full set of coordinates {t, x, 9a) is called 
Superspace. 

The functions defined in the Superspace are polinomic functions of the 9 variables (since 0^ =0). Thus 
we can decompose the functions (superfields) of this Superspace in components of 6''^, 9a, 9a9p, ... each 
of these components will be a function of the space-time coordinates. Analogously to the space-time, 
we can define in the Superspace scalar superfields, vector superfields, . . . For example in a 4-dimensional 
space time with A^ = 1 supersymmetry a scalar superfield has 10 components. 

There can be defined fields with specific properties with respect to the 9 variables. We are interested 
in the chiral fields. A scalar chiral field in a 4_D A^ = 1 Superspace has 4 components, two of them 
(the components of 9a) can be associated to be the components a Weyl spinor, the component of 9^ is 
a scalar field, and the 9a9f3 component is the so called "auxiliary" field. This auxiliary field is not a 
dynamical field since its equations of motion do not involve time derivatives. To this end we are left with 
a superfield, whose components represent an ordinary scalar field and an ordinary chiral spinor. So if 
nature is described by the dynamics of this field we would find a chiral fermion and a scalar with identical 
quantum numbers. That is Supersymmetry relates particles which differ by spin 1/2. Had we started 
with a A^ = 2 SUSY we would end with a set of particles of spin 0, 1/2 and 1 as a part of the same 
scalar superfield, this is called a Supermultiplet. When a SUSY transformation (Q) acts on a superfield 
it transform spin s particles into spin s ± 1/2 particles. 

Thus, for a, N ~ \ SUSY, we find that to any chiral fermion there should be a scalar particle with 
exactly the same properties. This fact is on the basis of the absence of quadratic divergences in boson 
mass renormalization, since for any loop diagram involving a scalar particle there should be a fermionic 
loop diagram, which will cancel quadratic divergences between each other, though logarithmic divergences 
remain. 

Supersymmetric interactions can be introduced by means of generalized gauge transformations, and 
by means of a generalized potential function, the Superpotential, which give rise to masses, Yukawa-type 
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interactions, and a scalar potential. 

As no scalar particles have been found at the electroweak scale we may infer that, if SUSY exists, 
it is broken. We can allow SUSY to be broken maintaining the property that no quadratic divergences 
are allowed: this is the so called Soft-SUSY-Breaking mechanism [^. We can achieve this by only 
introducing a small set of SUSY-Breaking terms in the Lagrangian, to wit: masses for the components of 
lowest spin of a supermultiplet; and triple scalar interactions. However, other terms like explicit fermion 
masses for the matter fields would violate the Soft-SUSY-Breaking condition. 

The MSSM is the minimal Supersymmetric extension of the Standard Model. It is introduced by 
means of a A^ = 1 SUSY, with the minimum number of new particles. Thus for each fermion / of the 
SM there are two scalars related to its chiral components called "sfermions" (/l._r), for each gauge vector 
V there is also a chiral fermion: "gaugino" (w), and for each Higgs scalar H another chiral fermion: 
"higgsino" [h). In the MSSM it turns out that, in order to be able of giving masses to up-type and 
down-type fermions, we must introduce two Higgs doublets with opposite hypercharge, and so the MSSM 



Higgs sector is of the so called Type II (see section |2.4.l| and Ref . [^ ) . 
We can define the following quantum number 

R = (^_l^|^s+L+zB ^ g ^ s]An , L = lepton number , B = barion number , 

called i?-parity, which is 1 for the SM fields and —1 for its supersymmetric partners. In the way the 
MSSM is implemented i?-parity is conserved, this means that _R-odd particles (the superpartners of SM 
particles) can only be created in couples, also that in the final product decay of an i?-odd particle at 
least one SUSY particle exists, and that the Lightest Supersymmetric Particle (LSP) is stable. 



In this chapter we revie w th e MSSM at the tree-level: its field conte nt (i n sec. 2J2); its Lagrangian in 



the Electroweak basis (sec. 2.3); its mass spectrum (sec. 2.4); in section 2.5 the interactions in the mass 



eigenstate basis; and finally we make a short revision of the experimental constraints on the parameters 



in section 2.6 



2.2 Field content 

The field content of the MSSM consist of the fields of the SM plus all their supersymmetric partners, and 
an additional Higgs doublet, so the superfield content of the model will be: 



the matter fields: 



L=[ ,-] R = il Q=[ ,] D = dl U = ul, 



^ \ B-r+ n- ^\ n _ J* fr- .,* (2-1) 



Y = -l Y = 2 Y=\ ^ = 1 ^ = "1 



for each generation of fermions 
• the gauge superfields, which in the Wess-Zumino gauge consist of: 



W^W^W^ W1W2W3 , 

i?o B« , (2.2) 

5m 9 , 



• 



and the two Higgs/higgsino doublets: 



/ 


' Hi \ 
H^ ) 




/ 


' Hi \ 
. H^ ) 


H2 = frQ 
\ ^2 


Y = 


= -1 


Y = 1 . 



(2.3) 



All these fields suffer some mixing, so the physical (mass eigenstates) fields look much different from 
these ones. The gauge fields mix up to give the well known gauge bosons of the SM, W^, Z" A^, the 



2.3 Lagrangian 



gauginos and higgsinos mix up to give the chargino and neutralino fields, and finally the Left- and Right- 
chiral sferniions mix among themselves in sfermions of indefinite chirality. Let aside the intergenerational 
mixing between fermions and sfermions that give rise to the well known Cabibbo-Kobayashi-Maskawa 
(CKM) matrix. For the sake of simplicity in most of our work we will take no intergenerational mixing, 
except in chapter H, where we make an analysis of some FCNC effects. 



2.3 Lagrangian 

The MSSM interactions come from three different kinds of sources: 
• Superpotential: 



W = e, 



fH\VR + hdHlQW + huHiQ'U - ^HIH^^ 



(2.4) 



The superpotential contributes to the interaction Lagrangian (2.11) with two different kind of 
interactions. The first one is the Yukawa interaction, which is obtained from (2.4) just replacing 
two of the superfields by its fermionic field content, whereas the third superfield is replaced by its 
scalar field content: 



Vy = e 



h.c 



JHIL'JR + hdHlQW + KHiQ'U - ^iHlHi 
fH{UR + hdHlQ^b + huHiQ'U 
fHiUR + hdHlQ^D + KHlQ^U 



(2.5) 



The second kind of interactions are obtained by means of taking the derivative of the superpotential: 



Vw^Y. 



dW{ip) 



dipi 



(2.6) 



(p,; being the scalar components of superfields. 
• Interactions related to the gauge symmetry, which contain: 

— the usual gauge interactions 

— the gaugino interactions: 

^GV-^ = *V2.9,^fer (T'^),; i,i + h.c. 



(2.7) 



where (iy9, t/") are the spin and spin 1/2 components of a chiral superfield respectively, T° is 
a generator of the gauge symmetry, Aq is the gaugino field and g'^ its coupling constant. 

— and the D-terms, related to the gauge structure of the theory, but that do not contain neither 
gauge bosons nor gauginos: 



with 

D ^ g iPi{T j,j fj , 

ifi being the scalar components of the superfields. 
Soft-SUSY-Breaking interaction terms: 



(2.9) 



K 



I 

soft 



g 



V2MwC0sP 



miAiHlVR + maAdHlQW - rriuAuHi^QU 



h.c. 



(2.10) 



The trilinear Soft-SUSY-Breaking couplings Af can play an important role, specially for the third 
generation interactions and masses, and they are in the source of the large value of the bottom 
quark mass renormalization effects (see section L4 ) . 



10 



The Minimal Supersymmetric Standard Model (MSSM) 



The full MSSM Lagrangian is then: 

^MSSM = -^Kinetic "1 



i-Gaugo ^G^ip 



Vd-Vy^Y. 



dW{Lp) 



d(fi 



-K 



soft 



ml H\Hi - ml hIH2 - mj^ (H1H2 + H\H. 



-^^u^.^l 



9^9 



-MwiW,- -M' B°B° 



-m\l*L 



mj^R*R-m^^Q*Q 



mjjU*U 



m\b*b 



(2.11) 



where we have also included the Soft-SUSY-breaking masses. 



From the Lagrangian (2.11) we can obtain the full MSSM spectrum, as well as the interactions, 
which contain the usual SM gauge interactions, the fermion-Higgs interactions that correspond to a Type 
II Two-Higgs-Doublet Model g^, and the pure SUSY interactions. A very detailed treatment of this 
Lagrangian, and the process of derivation of the forthcoming results can be found in M] . 



2.4 MSSM spectrum 

2.4.1 Higgs sector 

When a Higgs doublet is added to the SM there exist two possibilities for incorporating it, avoiding 
Flavour Changing Neutral Currents (FCNC) at tree level fS^. The first possibility is not to allow a 
coupling between the second doublet and the fermion fields, this is the so called Type I 2HDM. The 
second possibility is to allow both Higgs doublets to couple with fermions, the first doublet only coupling 
to the Right-handed down-type fermions, and the second one to Right-handed up-type fermions, this is 
the so called Type II 2HDM. 

The Hig gs se ctor of the MSSM is that of a Type II 2HDM ^, with some SUSY restrictions. After 
expanding ( ^.11 ) the Higgs potential reads 



V 



m\\Hi\ 



ml\H2? 



+ k9'+9"){\HA 



\H2 



(ey HI Hi + h.( 



2 5^1^1^21 



(2.12) 



The neutral Higgs bosons fields acquire a vacuum expectation value (VEV), 



(^i)o = 



(^2)0 = 



(2.13) 



We need two physical parameters in order to know their value, which are usually taken to be 



Ml^^\g\vl+vl)=g^\ , tan/3 =^ 
2 Z Vi 



(2.14) 



These VEV's make the Higgs fields to mix up. There are five physical Higgs fields: a couple of charged 
Higgs bosons {H"^); a pseudoscalar Higgs {CP = —1) A^; and two scalar Higgs bosons {CP = 1) H^ (the 
heaviest) and /i° (the lightest). There are also the Goldstone bosons G^ and G^ . The relation between 



the physical Higgs fields and that fields of (2.3) is 



H+ 



H? = 



HI 



-{cos/SG- - sin/3 i/-) , 
sin/?G++cos/3ff+ , 

sin a h^ 



vi + -^ (cos a H° 



i(cos/3G°-sin/3A")) 
W2 + -75 {sinaH^ + cosah° + i{smPG° + cosl3A°)) 



(2.15) 



where a is given in (2.17) 



All the masses of the Higgs sector of the MSSM can be obtained with only two parameters, the first 



one is tan/3 (2.14), and the second one is a mass; usually this second parameter is taken to be either 



the pseudoscalar Higgs mass Af^o or the charged Higgs mass Mh± ■ We will take the last option, as the 
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charged Higgs plays an important role in most of our studies. From ( 2.12| ) one can obtain the tree- level 
mass relations between the different Higgs particles, 



Ml„ = Ml± 



M' 



Mlo 



HO.h" 



W 1 



^ ( Alio + Ml ± ^{M%, + Ml) ' - 4 Ml, Ml cos2 2/3 



(2.16) 



and the mixing angle between the two scalar Higgs is obtained by means of: 



cos 2a 



cos 2/3 



M% - Ml 



M'ho 



Mlo 



sin 2a 



sin 2/3 



Mjj, + M2 



hO 



Mho 



Mlo 



(2.17) 



2.4.2 The SM sector 



In this section we give some expressions to obtain some MSSM parameters as a function of the SM 

parametrization. 

As stated above (sec. ^.4.1 ) the VEV's can be obtained by means of ( 2.14 ), and the Z mass can be 
obtained at tree-level by the relation: 



sin 9w — 1 



Ml 



Fermion masses are obtained from the Yukawa potential ( |2.5| ) letting the neutral Higgs fields acquire 
their VEV ( p. 13 ). The up- type fermions get their masses from the H2 whereas Hi gives masses to 
down- type fermions, so 



m„ = huV2 



huV^Mws\a.(3 



rud = hdVi 



hdV2Mw cos P 



and the Yukawa coupling can be obtained as 



A,„ = ^ = 



g V^Mw sin P 



, hd 

Arf = — = 



rUd 



V2Mw cos P 



(2.18) 



2.4.3 Sfermion sector 



The sfermion mass term is obtained from the derivative of the superpotential (2.6), the D-terms (UM) 
and the Soft-SUSY-Breaking terms (2.11) letting the neutral Higgs fields get their VEV (2.13), and one 
obtain the following mass matrices: 



Mi = 



A/? + m2 + cos2/3(r3«^ - g,s^)M| 



m,Ml^ 



rriqMl^ 



Ml^+mi^ 



cos 2/3 Q^s^Af 



(2.19) 



being Q the corresponding fermion electric charge, Tg^ the third component of weak isospin, Mq^^ ^ the 
Soft-SUSY-Breaking squark masses |p^p^ (by 5f7(2)i-gauge invariance, we must have Mi^ — Mg , 
whereas M^ , Mg are in general independent parameters), se — sin^^y, and 



MIr = A., 



^ cot /3 , M'Ir = Ad — ^i tan /3 



(2.20) 



We define the sfermion mixing matrix as (g'^ = {q' ^ = qL, q' 2 = te} ^-re the weak-eigenstate squarks, 
and qa = {(ji, 92} are the mass-eigenstate squark fields) 






E^i^-"' 



sin On 



ij(9)t^2^(9) ^ diag{m| , ml } {m^, > m^, ) , 



taii29„ — 



'LR 



' Ml - Ml + cos 2/3(T|- - 2Qqsl^)Ml " 



(2.21) 

(2.22) 
(2.23) 
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From eq. ( ^.19 ) we can see that the sfermion mass is dominated by the Soft-SUSY-Breaking parameters 
(Mir ^ TTi/ for / 7^ top) , and that the non-diagonal terms could be neglected, except in the case of the 
top squark (and bottom squark at large tan/3), however we will maintain those terms, the reason is that, 
although the A parameters do not play any role when computing the sfermion masses, they do play 
a role in the Higgs-sfermion-sfermion coupling -see eq. (2.40)-, and thus it has an effect on the Higgs 
self-energies. Moreover these A parameters are constrained by the approximate (necessary) condition of 
absence of colour-breaking minima. 



^2 < 3 (to? -I- to? -h Mfj + M^) , 



(2.24) 



6 13-0. 



where rriq is of the order of the average squark masses for q = t. 

All the Soft-SUSY-Breaking parameters are free in the strict MSSM, however some simplifications 
must be done to be able of making a comprehensive numerical analysis. As the main subject of study 
are the third generation squarks we make a separation between that and the rest of sfermions. This 
separation is justified by the evolution of the squark masses from the (supposed) unification scale down 
to the electroweak scale p^] (see also section 2.6.1 for a more detailed discussion). 



So we will use the following approximations: 



equality of the diagonal elements of eq. (2.19) 



Mh 



qu — -^gll - ---q 



Mf.22 



(2.25) 



for each charged slepton and each squark of the the first and second generation. 



the up and charm type sfermions share the same value of the parameter (2.25) 



the first and second generation squarks share the same value of the A parameter ( 2.20 ) 



sleptons also share the same value for (2.25) and A parameters 



2.4.4 Charginos and neutralinos 

Gauginos and higgsinos develop mixing due to the breaking of the gauge symmetry. To find the mass 
eigenstates we construct the following sets of two-component Weyl spinors 



r+ EE {-iw+,H+), 



(2.26) 



Then from ( |2.5D (higgsino mass parameter /i), the Soft-SUSY-Breaking masses (2.11) (gaugino mass 
terms M, M'), and replacing the Higgs fields by its VEV's in (p. 7]), we obtain the following chargino and 
neutralino mass Lagrangian 



C 



M 



(r+,r-) 







M 




where we have defined 



M = 



M 
\/2MwS\nl3 



;ri,r2,r3,r4)A^° 




\/2Mw cos 13 



X0 = 



-MzCos/3se\ 
Mz cos /3ce 



W MzsinPse 

M -Mz sin Pc0 

Mz sin Psg — Mz sin /3ce 

—Mz cos fise Mz cos Pee —fj. 

We shall assume a grand unification relationship between the gaugino parameters 

M' 5 







/ 



M 



■ tan^ 9 



w 



(2.27) 



(2.28) 



(2.29) 



(2.30) 
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The mass matrices (2.28) and ( 2.29| ) are diagonalized by 



U*MV'^ = Md = diag(Afi,M2) , 



where U, V and N are in general complex matrices that define the mass eigenstates 



*^ 






*; 



C^^ 






and 



*" ^ 






(2.31) 



(2.32) 



(2.33) 



In practice, we have performed the calculation with real matrices U, V and N, so we have been using 
unphysical mass-eigenstates (associated to non-positively definite chargino-neutralino masses). The tran- 
sition from our unphysical mass-eigenstate basis {^} = {"^f, ^!^} into the physical mass-eigenstate basis 
{x} = {Xi 7 Xa} '^^^ ^^ done by introducing a set of e parameters as follows: for every chargino-neutralino 
^I^ whose mass matrix eigenvalue are Mi, Ma, the proper physical state, x, is given by 



X = 



±75* 



if 
if 



(2.34) 



and the physical masses for charginos and neutralinos are fn ± — e Mi and m^i 



eM°, respectively. 
Needless to say, in this real formalism one is supposed to propagate the e parameters accordingly in all 
the relevant couplings, as shown in detail in Ref. [g3|0. This procedure is entirely equivalent [Q to 
use complex diagonalization matrices insuring that physical states are characterized by a set of positive- 
definite mass eigenvalues; and for this reason we have maintained the complex notation in all our formulae. 
Whereas for computations with real sparticles the distinction matters [g3|^; for virtual sparticles the 
e parameters cancel out, and so one could use either basis {*} or {x} without the inclusion of the e 
coefficients. We have stressed here the differences between the two bases just to make clear what are the 
physical chargino-neutralino states, when they are referred to in the text. 



2.5 Interactions in the mass-eigenstate basis 

We need to convert the interaction Lagrangian presented in section E^ to a Lagrangian in the mass- 
eigenstate basis, which is the one used in the computation of the physical quantities. As the expression 
for the full interaction Lagrangian in the MSSM is rather lengthy we quote only the interactions that we 
will need in our studies. Explicit Feynman rules derived from these Lagrangians can be found in [ |7l| . 

• fermion-sfermion~(chargino or neutralino): this int eraction is obtained from the potential ( p.7| ) 
-gauginos-, and form the Yukawa coupling term (2.5) -higgsinos-, in the mass-eigenstate basis: 



-C*.,- = .9 E E [-il^^(A^Z^^^PL+A^-\^PR)h 



a=l,2i=l,2 



,(b) 



(b) 



b:n{A';LPL+Az.e,pR t 



9 

V2 



E E H:*° (4L^^ + ^^L^«pk 

i=1^2a=l,...,4 
+ K^'i{Afao.PL+A^'L^aPB)b 



h.c. 



(2.35) 



where we have introduced the usual chirality projection operators Pl,r = 5 (l T 7^) and the 
matrices 



(t) 



— XbRia Vi2 , 
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4L = R^l*{K^+YLi^newK^) + V2MR^il*N:^, 
aL*L = V2At4*i*7V„3-r^^tan0w4a*Afai, 

^+a% - R-la ^1 ~ ^bn.2a V i2 : 

4(b) _ X R('')*7-r 

4'L = ^la * (K2 - Yl tan 0;^iV:i) - ^/2Ab4f iV:4 , 

aL'L = -\/2Abi?S*iV„4 + r^tan0M/4a*A^ai. (2.36) 

with 11 and Yp' the weak hypercha rges of the left-handed SU{2)l doublet and right-handed singlet 
fermion, and At and A;, are - Cf. eq.(E3q) - the potentially significant Yukawa couplings normalized 
to the SU{2)l gauge coupling constant g. 



quark-squark-gluino: the supersymmetric version of the strong interaction is obtained from (2.7): 
c~,,g = -^qh ^c in., (RifPL - RifPR) q, + h.c. , (2.37) 

where A'^ are the Gell-Mann matrices. 

quark-quark-Higgs: this is the usual Yukawa interaction from Type II 2HDM, in the MSSM it 
follows after replacing in (|2.5|) the mass-eigenstate Higgs fields (2.15|): 



9_ 

4 

grrid 

2MwCoaP 

grriu 



^H+ud = —/= [m (m„ cot/3Pt + "^d tan/3Pfl) diJ^ + h.c] 

[(cosaij" - smah°) dd - i sin (3 d'j5d A°] 

[(sinaiJ° + cosa/i°) UM — icos/?u75uA°] , (2.38) 



2MwsmP 

here we have replaced the Yukawa couplings hi in favour of masses and tan/3. 
• squark-squark-Higgs: the origin of this interaction is twofold, on one side the superpotential deriva- 



tive (2.6), and on the other the Soft-SUSY-Breaking trilinear interactions. 



^^"^^' 'vi^"*"^'^"'''' 



where we have introduced the matrixU 

^ / M^sin2/3- (m^tan/3 + m^cot/3) -m^ (a^ + A^ tan/3) \ -^ ^„^^ 

■' \ — niu ifi + Au cot (3) —m„md (tan /3 4- cot /3) J ' 



• 



chargino-neutralino-charged Higgs: this interaction is obtained from (2.7), we note that in the 



electroweak basis the only interaction present is the Higgs-higgsino-gaugino one 

Ch±^t^o - -gH+i^+ {sin pQ'^*ec.PL+ cos PQ'^*e.PR) ^° , (2.41) 

r Q^, = [/AiV*3 + ^ (iV*2 + tan OwK^) U*, 






rll = V,,N^4--^{N^2 + tan9wN^i)V,2 



• gauge interactions: in this Thesis we only need a small subset of the gauge interactions present 
in the MSSM, so we will only quote the interactions of the W^ boson, and those of QCD. The 
photon interactions are simply those of QED (and scalar QED). For a complete set of the Z boson 
interactions see for example PTII 



^Note that our convention for the fj, parameter in (2.4) is opposite in sign to that of |33 
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quark-W^=' 



squark-14^=' 



Cw+ud = -^ u^PLd W+ + h.c. , (2.42) 



^w+ui = i^Rt^*R[fw+u: d^ dt + h.c. , (2.43) 

chargino-neutralino-iy^ : 

Cw+x+x" - 5^5" 7^ (C^.eae.Pi + C^iPR) ^t W' + h.c. , (2.44) 

Higgs-W^: after SSB there exist three different kind of gauge interactions for the Higgs (and 
Goldstone) bosons p9| , namely triple interactions of a gauge boson and two scalars, triple 
interaction of two gauge bosons and a scalar, interaction between two gauge bosons and two 
scalars. We only need the first one of these interactions to perform the analysis presented here, 
that is 

^w+H-H = ^W+ f- sin {(3- a) H-id^,H" + cos {P- a) H- 1 d^h" 

-H- df, aA + h.c. . (2.45) 



quark strong interactions: this is the usual QCD Lagrangian 

CQCD = ^G';,X'r^q,rqj ■ (2.46) 

squark strong interactions: aside from the well known scalar QCD Lagrangian, the scalar 



potential (2.S) introduces quartic scalar interactions between squarks of order a^, thus we 
have 



^Gqq = -i^ G^ Ay- qli d^ qa.j , 



,5s 
2 



^qqqq — "^ ^(-^la Rib — R2a R2b) [Rlc Rid — R2c R2d) X 

S,i Kj Qb,j q'c*k Ki Qd,i ■ (2.47) 

2.6 MSSM parametrization 

2.6.1 MSSM parameters 

If SUSY were an exact symmetry then only one parameter should be added to the SM ones (tan (3) , but 
we have to deal with a plethora of Soft-SUSY-Breaking parameters, namely 

• masses for Left- and Right-chiral sfermions, 

• a mass for the Higgs sector, 

• gaugino masses, 

• triple scalar couplings for squarks and Higgs. 
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This set of parameters is often simplified to allow a comprehensive study. Most of these simplifications 
are based on some universality assumption at the unification scale. In minimal supergravity (mSUGRA) 
all the parameters of the MSSM are computed from a restricted set of parameters at the Unification scale, 
to wit: tan/9; a common scalar mass ttiq; a common fermion mass for gauginos mi/2', a common trilincar 
coupling for all sfermions Aq; and the higgsino mass parameter /x. Then one computes the running of 
each one of these parameters down to the EW scale, using the Renormalization Group Equations (RGE), 
and the full spectrum of the MSSM is found. 

We will not restrict ourselves to a such simplified model. As stated in the introduction we treat the 
MSSM as an effective Lagrangian, to be embedded in a more general framework that we don't know 
about. This means that essentially all the parameters quoted above are free. However for the kind 
of studies we have performed there is an implicit asymmetry of the different particle generations. We 
are mostly interested in the phenomenology of the third generation, thus we will treat top and bottom 
supermultiplet as distinguished from the rest. This approach is well justified by the great difference of 
the Yukawa couplings of top and bottom with respect to the rest of fcrmions. We are mainly interested 
on effects on the Higgs sector, so the smallness of the Yukawa couplings of the first two generations will 
result on small effects in our final result. We include them, though, in the numerical analysis and the 
numerical dependence is tested. On the other hand, if we suppose that there is unification at some large 
scale, at which all sfermions have the same mass, and then evolve these masses to the EW scale, then 
the RGE have great differences [Q. Slepton RGE are dominated by EW gauge interactions, 1st and 
2nd generation squarks RGE are dominated by QCD, and for the 3rd generation squarks there is an 
interplay between QCD and Yukawa couplings. Also, as a general rule, the gauge contribution to the 
RGE equations of left- and right-handed squark masses are similar, so when Yukawa couplings are not 
important they should be similar at the EW scale. 

With the statement above in mind we can simplify the MSSM spectrum by taking an unified para- 
metrization for 1st and 2nd generation squarks (same for slcptons). We will use: a common massn for 
Ui and qn {rriu); an unified trilinear coupling A^ for 1st and 2nd generation; a common mass for all 1)^ 
and Ij^ (Tif ); and a common trilinear coupling ArQ. 

For the third generation we will use different trilinear couplings At and A},, as these can play an 
important role in the kind of processes we are studying (see chapter H). Stop masses can present a large 
gap (due to its Yukawa couplings), being the right-handed stop the lightest one. We will use a common 
mass for both chiral sbottoms, which we parametrize with the lightest sbottom mass (mj^ ), and the 
lightest stop quark mass (toj ), as the rest of mass inputs in this sector. This parametrization is useful 
in processes where squarks only appear as internal particles in the loops (such as the ones studied in 
chapters ^ and ^, as one-loop corrections to these parameters would appear as two-loop effects in the 
process subject of study. However in chapter we deal with squarks as the main subject of the process 
and in this case a more physical set of inputs must be used. We have chosen to use the physical sbottom 
masses (mg , mj^ ) and the sbottom and stop mixing angles {9^^, %) to be our main inputs. Again one-loop 
effects on other parameters (such as A},) would show up as two- loop corrections to the observables we are 
interested in. 

For the same reasons EW gaugino sector is also supposed to have small effects in our studies. Thus 



the grand unification relation introduced in expression (2.3C). Gluino mass (nig), on the other hand, is 
let free. 

For the Higgs sector two choices are available, we can use the pseudoscalar mass M^o, or the charged 
Higgs mass Mfj± . Both choices are on equal footing. As the charged Higgs particle is a main element for 
most of our studies we shall use its mass as input parameter in most of our work. However in chapter ra 
it is more useful to use M^o. 

Standard model parameters are well known, we will use present determinations of EW observables 

Mz = 91.1867 ± 0.0021 GeV 

Mw = 80.352 ±0.054 GeV 

Gf = (1.16637 ± 0.00001) X 10"^ GeV-2 

nir = 1777.05 ± 0.29 MeV 

a{Mzy^ = 128.896 ± 0.090 . (2.48) 



^Note that aftetrfiagonalization of the squark mass matrix the physical masses will differ slightly. 
•^See section 2.4.3 for the concrete definitions of these parameters. 
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QCD related observables are not so precise. On the other hand as the main results are not affected by 
specific value of these observables we will use the following ones 



nit ~ 


= 175 GeV 


rrib = 


= 5 GeV 


as{mt) = 


= 0.11 



(2.49) 
(the last figure corresponds to as{Mz) — 0.12). 

2.6.2 Constraints 

The MSSM reproduces the behaviour of the SM up to energy scales probed so far [Q . Obviously this is 
not for every point of the full parameter space! 

There exists direct limits on sparticle masses based on direct searches at the high energy colliders 
(LEP II, SLC, Tevatron). Although hadron colliders can achieve larger center of mass energies than 
e^e~ ones, its samples contain large backgrounds that make the analysis more difficult. This drawback 
can be avoided if the ratio signal-to-background is improved, in fact they can be used for precision 
measurements of "known" observables (see e.g. [Q). e+e" colliders samples are more clean, and they 
allow to put absolute limits on particle masses in a model independent way. 

The most stringent bound to the MSSM parameter space is the LEP II bound to the mass of charged 
particles beyond the SM. At present |77| - [79| this hmit is roughly 

Mchargcd > 90 GeV . (2.50) 

Specific searches for Supersymmetric particles are being performed at LEP II, negative neutralino 
searches rise up a limit on neutralino masses of |7|] 

M^o > 30 GeV , (2.51) 

it turns out that after translating this limit to the fi — M parameters it is less restrictive than the one 



obtained for the charginos from (2.5C) 



Actual Higgs searches at LEP II imply that, for the MSSM neutral Higgs sector [gO[ 

Mho > 72.2 GeV , M^o > 76.1 GeV . (2.52) 

Notice that without the MSSM relations there is no model independent bound on M^o from LEP 
pW. Actual fits to the MSSM parameter space project a preferred value for the charged Higgs mass 
oiMH± c^ 120 GeV §|. 

Hadron colliders bounds are not so restrictive as those from e'^e~ machines. Most bounds on squark 
and gluino masses are obtained by supposing squark mass unification in simple models, such as mSUGRA. 
At present the limits on squarks (1st and 2nd generation) and gluino masses are |^^ 

rrig > 176 GeV , rrig > 173 GeV . (2.53) 

From the top quark events at the Tevatron a limit on the branching ratio BR{t — > iJ+ b) can be 
extracted, and thus a limit on the tan/3 — M^± relation. We will treat this limit in detail in chapter 0. 

Finally indirect limits on sparticle masses are obtained from the EW precision data. We apply these 
limits through all our computations by computing the contribution of sparticles to these observables and 
requiring that they satisfy the bounds from EW measurements. We require new contributions to the p 
parameter to be smaller than present experimental error on it, namely 

Spnew < 0.003 . (2.54) 

We notice that as Spncw is also the main contribution from sparticle contributions to Ar p7[ |, new 
contributions to this parameter are also below experimental constrains. Also the corrections in the a- 



and GF-on-shell renormalization schemes will not differ significantly (see section B.l). 

There exist also theoretical constrains to the parameters of the MSSM. As a matter of fact the MSSM 
has a definite prediction: there should exist a light neutral scalar Higgs boson h°. Tree-level analysis put 
this bound to the Z mass, however the existence of large radiative corrections to the Higgs bosons mass 
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relations grow this limit up to ~ 130 GeV. Recently the two-loop radiative corrections to Higgs mass 
relations in the MSSM have been performed ||83|-p5[, and the present upper limit on M^a is 

Mho < 130 - 135 GeV . (2.55) 



The two figures in (2.55|) have been computed by difi^erent groups [p3hB5[ and there is a great interest 



in make them match [ p5[ . It is very important to know as precise as possible this limit, as by means of 
a possible Run III of the Tevatron collider (TEV33, at the same energy, but higher luminosity) either a 
h° should be found, or on the contrary a lower limit to its mass in the ballpark of 130 GeV will be put. 
Thus it is of extreme importance to have both, a very precise prediction for the bound ( 2.55|) , and a very 



precise analysis of the Tevatron data. Of course if the MSSM is extended in some way this limit can be 
evaded, though not to values larger of ~ 200 GeV 



Another theoretical constraint is the necessary condition (2.24) on squark trilinear coupling (A) to 
avoid colour-breaking minima. This constraint is easily implemented when the A parameters are taken 
as inputs, but if we choose a different set of inputs (such as th e mixing angle 6q, as in chapter g) then it 
constrains the parameter space in a non-trivial way -eq. (|6.10 ). 



Whatever the spectrum of the MSSM is, it should comply with the benefits that SUSY introduces 
into the SM which apply the following condition is fulfilled: 

MsvsY < Oil TeV) . (2.56) 

If supersymmetric particles had masses heavier than the TeV scale then problems with GUT's appear. 
This statement does not mean that SUSY would not exist, but that then the SM would not gain practical 
benefit from the inclusion of SUSY. 

A similar upper bound is obtained when making cosmological analyses, in these type of analyses one 
supposes the neutralino to be part of the cold dark matter of the universe, and requires its annihilation 
rate to be sufficiently small to account for the maximum of cold dark matter allowed for cosmological 
models, while at the same time sufficiently large so that its presence does not becomes overwhelming. 
Astronomical observations also restrict the parameters of SUSY models, usually in the lower range of the 
mass parameters (see e.g. Q). 

For the various RGE analysis to hold the couplings of the MSSM should be perturbative all the way 
from the unification scale to the EW scale. This implies, among other restrictions, that top and bottom 
Yukawa couplings should be below certain limits. In terms of tan/3 this amounts it to be confined in the 
approximate interval 

.5 < tan/3 < 70 . (2.57) 

All these restrictions will apply in all our numerical computations. Any deviation from this framework 
of restrictions will only be for demonstrational purpouses, and will be explicitly quoted in the text. 



Chapter 3 

MSSM renormalization 



3.1 Introduction 

In this chapter we perform the renormahzation of the MSSM in the on-shell scheme. We do not pretend 
to make all the renormalization procedure, but just sketching what are the necessary ingredients of this 
renormalization and giving expressions for some non-SM two-point functions. The renormalized three- 
point Green functions are the subject of the forthcoming chapters. We will not give the full expressions for 
the gauge bosons self-energies, or the 5p and Ar parameters, since these have been subject of dedicated 
studies p% pq, |37|, |6q, |89|-|9^ ■ O^i the other hand the various counterterms and self-energies given in this 
chapter are general. We have left some expressions out of this chapter as they are approximations only 
valid in the context where they are used (see chapter o) . 

We address the renormalization of the MSSM extending the SM on-shell procedure described in 
|28,p9,p9-MR. We may use both the a or the Gp parametrizations. At one-loop order, we shall call the 
former the "a-scheme" and the latter the "Gi?-scheme" . In the "a-scheme" , the structure constant a = 
oiemitf' = 0) and the masses of the gauge bosons, fermions and scalars are the renormalized parameters: 
(a, Mw,Mz^ MH,mf,MsuSY, ■ ■ •) -Msusy standing for the collection of renormahzed sparticle masses. 
Similarly, the "GF-scheme" is characterized by the set of inputs {GF,Mw,Mz,MH,mf , Msusy, ■■ ■)■ 
Beyond lowest order, the relation between the two on-shell schemes is given by 



V2 2M, 



(1 + Ar*^^^*'0 , (3.1) 



vVw 

where Ar^^"^"^*^ is the prediction of the parameter Ar [||,|^,||l in the MSSM0. 

Let us sketch the renormalization procedure affecting the parameters and fields related to the various 
processes subject of study. In general, the renormalized MSSM Lagrangian £ ^ C + S£ is obtained 
following a similar pattern as in the SM, i.e. by attaching multiplicative renormalization constants to each 

1/2 

free parameter and field: g^ -^ {l + Sgi/gi)gi, ^i -^ Z^, ^i. As a matter of fact, field renormalization (and 
so Green's functions renormalization) is unessential and can be either omitted or be carried out in many 
different ways without altering physical (S'-matrix) amplitudes. In our case, in the line of Refs. pO| , [4l| , we 
shall basically use minimal field renormalization, i.e. one renormalization constant per gauge symmetry 
multiplet |2^,^,^. In this way the counterterm Lagrangian, 5C, as well as the various Green's functions 
are automatically gauge- invariant. 

The specific sign convention of the various two-point functions used all over this Thesis is based on 
the prescription that the unrenormalized self-energy always add up to the bare mass parameter (or the 
squared mass, depending on the kind of particle), which is equivalent to say that, in the on-shell scheme, 
the mass parameter counterterm is minus the unrenormalized self-energy, that is 



m° 



-I- Re (S(fc2)) = m -f- (5m + Re (^(fc^) = O) , Sm = -Re (i;(fc^)) 



where m is the physical mass parameter -the mass for fermions, the squared mass for bosons- and 6m the 
corresponding counterterm (see next sections for the concrete definition in each case). The convention 
for each kind of particle can be seen in table |3.l| . 

^ Our conventions diflfcr from those of Pg, EQJ] . 

2^j.MSSM j^g^g jjggj^ subject of dedicated studies, see [g§|89||3|. 
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MSSM renormalization 



fermion 



-m-E(fe2) 



^-rn ^-m ^ ' ^ f- 



scalar 



fe^-m^-S(fc^) — fc^-m^ + fe^-m^ V *^(''^ )j fe^-m^ 



gauge boson 






Table 3.1: Self-energies sign conventions for the various kind of particles. The gauge bosons are dealt 
with in the Feynnian gauge. 



For the regularization of the ultraviolet divergent integrals we use the Dimensional Reduction (DRED) 
|£^,Q prescription, as it respects SUSY. As a matter of fact one-loop computations with only i?-even 
external particles yield the same result in DRED and Dimensional Regularization, however this is not 
necessary true for higher loop computation, or for computations with i?-odd external particles. 

3.2 A note on the gauge sector renormalization 



For the sake of fixing notation, in this section we review some well known features of the renormalization of 
the electroweak gauge sector, which is identical to the SM one. We refer to [^,|2^,^ for a comprehensive 
exposition of the subject, and to [p4|,pq,|37|,|66|,|90|-p2|,|94[| for the MSSM expressions of the various self- 
energies. 

For the SU{2)l gauge field we have 



VF,, 



Mw 



(3.2) 



1 -|- 5Z^ is the usual SU{2)l gauge triplet renormalization constant given by the formula 



5ZY 



S^(fc2) 



fc2 



fe2=0 



) Cty S7z(0) c- 
sw Mj. s 



5Ml 
1^ 






and 



SM? 



w 



-E^(fc2 = M^) , 



(5M| 



-Ez(fc2=M|), 



(3.3) 



(3.4) 



are the gauge boson mass counterterms enforced by the usual on-shell mass renormalization conditions. 
Th e S functions denote the (real part of the) unrenormalized two-point Green functions. 6Zhw on 
eq.(3^) is a dimensionless constant associated to the wave-function renormalization mixing among the 
bare H^ and W^ fiel ds; its meaning and value is discussed together with the Higgs renormalization 
procedure (section 3.4). 

For the SU{2)l gauge coupling constant, we have 



5 - (1 + ^).9 = (z^) [z^r"" 



9, 



(3.5) 



where Z|^ refers to the renormalization constant associated to the triple vector boson vertex. Therefore, 
from charge renormalization. 



5a 
a 



S^(fc2) 



Sw ^iziO) 



fc2=0 



cw 



Ml 



(3.6) 



and the bare relation a — g^ s^ /An -^ a + 5a = {g^ + 5g^) (s^ -I- (5s^)/47r, one gets for the counterterm 

to g: 

x„2 x„, „2 /xnj2 xn,r2 \ 

(3.7) 



5a 
a 



c'w f5Ml 



^w 



\Ml 



5M^ 
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and as a by-product 






(3.8) 



3.3 Fermion renormalization 



Following the directives from section 3.1 and references M,Bl| we introduce the fermion wave function 
renormalization constants 



A/2 / tL 



n1/2 



'{zif't,^' 



(Z'n) ' bn, tn^{Z\,y'-ti 



(3.9) 



7t,b\ 



Here Z^ = 1 + 5Zi are the doublet {Zl) and singlet {Z^ ) field renormalization constants for the top 
and bottom quarks. Although in the minimal field renormalization scheme there is only one fundamental 
constant, Zl, per matter doublet, it is useful to work with Z^ = Zl and Z£, where the latter differs from 
the former by a finite renormalization effect |2^, g^, 95|. To fix all these constants one starts from the 
usual on-shell mass renormalization condition for fermions, /, together with the "residue = 1" condition 
on the renormalized propagator. These are completely standard procedures, and in this way one obtains"^ 



5m f 
nif 



— —^ — — + m^f) 



(3.10) 



and 



SZin 



- v/ 



^L,R 



f'f 



-f'( 



(mj) + mjiEi'imj) + E^'(mj) + 2^'imj)] 



where we have decomposed the fermion self-energy according to 



(3.11) 



(3.12) 



and used the notation S'(p) = dYj{p) / dp^ . 

The one- loop Feynman diagrams contributing to these various self-energies can be seen in figures 3.1 
and 3.2 for the bottom and top quarks respectively. To express the various self-energies and vertex 
functions we use the standard one-, two- and three-point one-loop functions from Refs. |100-103| which 
we have collected in Appendix M. Using this notation the bottom quarks self-energies read as 



'{L,R} 



ip' 



= E 



{L,R} 



-ig 



ip') 



i(t) 



Bi{p,M,,mi 



Ab) 



' B^{p,Mlmi^ 



mbT,s{p'^ 



™&S^(P^)|(a 



(a) + (fc) 



«5 



M.Re (4i:AL*L) Bo (p, M„ mJ 
iM°Re (A^^tAfaa) Bo (p, M°, m^^ 



from SUSY-EW particles, and 



^iL,R}ip') 



= S 



{L,R} 

g' 



ip' 



ic) + {d) 



2cos2/3 



{™{t,&} [{cot2/3, tan2/3}Bi {p, nit, Mh± ) + Bi {p, nit , M, 
[cos^a Bi{p, TTif,, Mffo) + sin^a Bi{p, mt,, M^o) 



(3.13) 



We understand that in all formulae defining counterterms we are taking the real part of the corresponding functions. 
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ba 



ta 




H+, G^ 



^ b 



(c) 



C^ 



(e) 



t 
(d) 

ba 




Figure 3.1: One-loop Feynman diagrams contributing to the bottom quark self-energy. 



E|(p2) = Y.'sip^ 



+ sin2/3Bi(p,TO6,M^o)+cos2/3Bi(p,mf„Mz)]} , 



9 



2tM^ 



{mf [Bo{p, mt, Mh±) ~ Bo{p, mt, Mw)] 



- — ^ Icos^ a Bo{p,mb, Mho) + sin^ a Bo{p, rub, Mho) 
icos^p 

- sin^P Bo{p, nib, Mao) -cos'^PBaip, mb,Mz)]} , 



(3.14) 



from Higgs and Goldstone bosons in the Feynman gauge. To obtain the corr es pond ing expressio ns for 
an up-like fermion, t, just perform the label substitutions 6 <-+ t on eqs. ( 3.13 )-( pl^ ); and on eq. ( 3.14 ) 
replace sin a <-^ cos a and sin/3 <-^- cos/3 (which also implies replacing tan/3 <-^ cot/3). 
The "strong" (QCD) self-energies from Figs.BJ^ (e) and (f) are 



^LM}ip' 



s?. .>(/: 



{LM} 

— —iSTrasCp 



Bi{p,mb,X) 



R 



ib) 
{l,2}a 



mg 



{Bq- Bi){p,mr,mg 



?(ft)p(6)*^ 



2 Boip, mb, A) + -^Re{R\">R}Z)B^{p, m^ , rUg) 
nib 



(3.15) 



where Cp = {N^ — l)/2 Nc — 4/3 is a colour factor and we have introduced a small gluon mass A to 
regularize the i nfrare d divergences. The top quark ones from Figs. 3.2 (e) and (f) are easily obtained by 
substituting in ( 3.15 ) the particle indexes b ^ t and b ^ t. 
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ta 



ba 




^ t 

-1 k. 



t 

(c) 



b 
(d) 



t 


o 


t 


^ 




t 

(e) 




9 

(f) 



Figure 3.2: One-loop Feynman diagrams contributing to the top quark self-energy. 



3.4 Higgs sector 



One also assigns doublet renormalization constants to the two Higgs doublets (2^) of the MSSM: 



Hi 
H- 






Hi 
H- 






71/2 






(3.16) 



Following the on-shell procedure we prefer to fix the counterterms using the physical fields. In this 
approach we have decided to take the charged Higgs as the renormalized particle, for the charged Higgs 
mass will be a natural input in most of our computations. This will induce finite renormalization constants 
in the other fields of the sector (^4^^, hP and 7?°). Of course other equivalent choices can be made. For 
example we could have taken A^ to be the renormalized Higgs as in [ 104- 107 1 and in this case the charged 
Higgs sector would have received finite renormalization effects. To this effect we introduce wave function 
renormalization constants for the Higgs particles in the mass-eigenstate basis Z[j±, Zq±, . . ., which are 
only shortcuts for certain combinations of Zhi and ZH2 1 and fix these constants as usual in the on-shell 
scheme by using as input particle the charged Higgs. We will discuss in detai l the charged Higgs sector 
renormalization whereas the neutral sector has been extensively discussed in |I04-107|. We will expose 
two equivalent approaches: the renormalization in the Feynman gauge and in the Unitary gauge. 



3.4.1 Feynman gauge 

As we carry out our calculations in the Feynman gauge, we would also like to perform the renormalization 
of the Higgs sector in that gauge. The Lagrangian is sketched as follows: 

L^Lc + Lgf + Lfp , (3.17) 

where Lq is the classical Lagrangian, Lgf stands for the gauge-fixing term in that gauge, 

Lgf^-F+F- + ... (i^± = a^W+TiMvvG+) , (3.18) 
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H^ 



G^ 




i/^ 



(a) 



(b) 



Figure 3.3: The renormalized mixed blobs W^ — H^ and G+ — H^ at any order of perturbation theory. 



and Lpp = fj"" {dF°' / 89^^ if is the Faddeev-Popov ghost Lagrangian constructed from FP and anti-FP 
Grassmann scalar fields ry, fj. Since we are interested in the charged gauge-Higgs [W^ — H^) and charged 
Gol dston e-Higgs (G^ — H^) mixing terms in that gauge, we have singled out just the relevant term on 
eq.(|3l^). 

As is well-known, although the classical Lagrangian, Lc, also contains a nonvanishing mixing among 
the weak gauge boson fields, W^ , and the Goldstone boson fields, G*, namely 



Caw = iM\ 



w 






h.c. 



(3.19) 



the latter is canceled (in the action) by a piece contained in Lcf- Now, after substituting the renor- 
malization transformation for the Higgs doublets, eq.( 3.16| ), on the Higgs boson kinetic term with 
SU{2)l X U{1)y gauge covariant derivative, one projects out the following relevant counterterms 



5C = 5Zh± df,H+d^'H- + 6Zg± df.G+d^'G- 

+ SZhg {df,H-d^'G+ + h.c.) + 6Zhw {iMw W^d^'H^ 



where 



and 



h.c. 



SZh±\ 
SZg±J 



„2 2 

„2 „2 

■5/3 "^/3 



SZhi 
5Zh2 



(3.20) 
(3.21) 



SZhg = Sf3 Cfs {5Zh.^ - SZh^ ) , 
5Zhw = Sf3 Cfj 






(3.22) 



behave as if they were renormalization constants introduced as 



G± 



/ Zh± SZhg 
\ SZhg ^g± 



G± 



(3.23) 



and 



VF,^ 



Z 



Y) 



1/2 



W^±i 



sz. 



HW 



Mw 



d^H^ 



(3.24) 



with Z^ being the usual 5*17(2)^ gauge triplet renormalization constant (3.3). A note regarding neutral 
Higgs particles is worth here. As the pseudoscalar Higgs and neutral Goldstone boson undergo the same 
mixing procedure as their charged partners, the same procedure above can be used for the A" — G*^ — Z'j^ 
sector, with the only proviso that a factor 1/2 must be put in front of the definition of the mixing terms 
SZag and SZaz in (3.23) and (3.24) to take into account the neutral nature of the particles. 
The renormalization transformation for the VEV's of the Higgs potential ( 2.12 ), 



.1/2. 



Svi 



.H.'i^. + Sv^)= [1 + ^+ SZh, 
' Vi 2 



(3.25) 
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implies that the counterterm to tan /3 is related to the fundamental counterterms in the Higgs potential 

byQ 

^ = ^ o [dZH2 - oZhJ ■ (3.26) 

tan/5 U2 -wi 2 

If one imposes the usual on-shell renormalization conditions for the A"-boson, one has 

tan (3 + cot (3 



SZh., — SZh, — — 



Ml 



■^az{MXo). 



(3.27) 



There exists also another mixing term between H^ and G^ originating from the mass matrix of the 
Higgs sector [104-106, 108 1. This one-loop mixture is contained in: 






H- 
Q- 



^ \/2v'' \/2v*' 

where we have attached a superscript b to bare quantities, and ti are the tadpole counterterms 

io = — sin (/3 — a) t^o + cos (/3 — a) t^o , 
ti = sin (/? — a) t^o + cos (/? — a) 1^° ■ 



(3.28) 



(3.29) 



We are now ready to find an expression for the mixed 2-point Green functions (Figs. 3.3(a) and 3.3(b)). 
For the W^ — H^ mixing (Fig. 3.3(a)) we can write the renormalized 2-point Green function as 



,HW 



fc2 - M^± 



^, ~zW(P) _^^^, 5Znw 



M, 



w 



LW 



M^ 



w 



-Wi^u 



fc2 - M^ 



which allows to define a renormalized self-energy as follows 



(3.30) 



(3.31) 



Now we must impose a renormalization condition on S (fc^); and we choose it in a way that the 
physical Higgs does not mix with the physical W^ : 



^HwiMfj±) = =^ SZhw = 



^Hw{Mjj±) 



M^ 



(3.32) 



Notice also that with this renormalization procedure on-shell W^ do not mix also with H^ since the 
renormalized 2-point Green function ( 3.30| ) is proportional to the external momentum /c". 
We still have another ingredient, the mixed H^ — G* 2-point Green function: 



A 



HG 



k^ - Mfj^ 



^li:HG{k') + lk'6ZHG~l^ 



V2vV k-^-Ml 



This allows to define renormalized self-energy 



(3.33) 



tnaik'') - S^G(fc') - kHZHG + -^ ■ 

The mixed self-energies f^Hwik'^) and Y.HG{k^) obey the following Slavnov- Taylor identity: 

e±Hw{k^) - M^±HG{k^) ^ . 



(3.34) 



(3.35) 



This identity is derived from a BRS transformation involving the Green function constructed with an anti- 
FP field and the charged Higgs field: < 0\Sbrs {v~^ ^''')|0 >= 0. Following the standard procedure 109 
one immediately gets: 



< 0|F+ H+\0 >=< 0\d''W-H+ - iMw G-H+\0 >= , 



(3.36) 



*For a more detailed discussion on the 5tan/3 counterterm, see Sec. 3.4.4 
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which in momcntuni space reads 



fc'^A 



HW 



with 



A"G = 



fc2 



Mi. 



MwA"'^ =0 , 



-tJ^HGik^) 



fc2 



M^ 



(3.37) 



(3.38) 



Clearly, eq.(3.37) implies eq.(3.35). The latter identity guarantees that the contribution from diagrams 
with external charged Higgs particles in Figs. p.3| (a) and 3.3(b) vanishes since no mixing is generated 
among the physical boson H"^ and the renormahzed fields G^ and W^: E//g(M^±) = 12HwiMfj±) — 0. 
There is of course another Slavnov- Taylor identity, derived in a similar manner, which insures that the 
renormalized mixing between G^ and W also vanishes. 

3.4.2 Unitary gauge 

It is useful also to have a look at the renormalization procedure in the Unitary gauge as well, where it is 
straightforward. Introducing the counterterm (3.24) into the Unitary gauge Lagrangian one obtains: 



Lug = ^jFfj.iyF'^'^ 



+W-'" 



M^w;[w 



La = tMwSZHw{W-df'H+ - W+d^'H-) 



In this gauge the corresponding renormalized 2-point Green function reads (Fig. 3.3(a)): 



fc2 - Af^^ 



-iT,Hw{k'^) 



M 



w 



ik" M^ 



M^ 



w 



ydllf 



'M^r 



fc2 - M^ 



(3.39) 



(3.40) 



which is identical to (3.30) but with the W^ propagator in the Unitary gauge. Thus a renormalized self- 
energy can be defined with the same formal expression as ( ^.31 ) and so it applies the same renormalization 
condition to obtain (3.32). 

Thus we have proven that the expression for SZhw is formally the same in both Unitary and Feynman 
gauges, but that in the latter gauge one must take into account the additional renormalization of the mixed 
self-energy T,hg- Moreover, it is possible to use different gauge fixing for particles inside and outside the 
loops | 110| , 111|, so we can use the Unitary gauge renormalization, but still maintain Goldstone bosons 
inside loops. 

3.4.3 Higgs masses and wave functions 

Whether in the Feynman or in the Unitary gauge, the charged Higgs counterterms can be introduced as 

H^ -^ {ZH±y/^H^ , Mh± -^ Mh± +5Mh± , (3.41) 

from which a renormalized self-energy can be defined as follows: 

EH±(fc') = ^H±{k')+5Mh± - (fc2 - MI±)5Zh± , (3.42) 

where I]^±(fc2) is the corresponding unrenormalized self-energy. 

In order to determine the counterterms, we impose the following renormalization conditions: 
i) On-shell mass renormalization condition: 

Sh±(A4±)=0, 

ii) "Residue = 1" condition for the renormalized propagator at the pole mass: 

dtH±{e) 



(3.43) 
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From these conditions one derives 



SMfj± 



5Z 



H± 



S'^±(A/^±)=0. 



^I]h±(Mh±' 
-S^± (A'/^± ) 



(3.44) 



(3.45) 
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Having fixed 5Mu± , 5Zu± , 5Zhw , and 5 tan (3/ tan /3 (see section 3.4.4 below) all the renormalization 
constants of the Higgs sector are fixed, and one can find the value of the original counterterms ( |3.16 ) by 
inverting the set of equations (3.21) and (3.22). 

On the other hand, the renormalization of the neutral Higgs sector has been studied in the series 
of works |104-107|. We can use the expressions for the one-loop neutral Higgs masses given in these 
references by translating the corrections for the charged Higgs mass to corrections to the pseudoscalar 
Higgs mass, that is, from the relations 



M 



H± 



^^ffilxrc 



Ail/4±(M^o,AfsusY) , 



Mlo = m2„(M^o,Msusy) , 
Af|„ = Mfjo{MAO,MsusY) , 

we invert the first equation above, and use the computed value of M'^o as input for the other ones. 



3.4.4 tan/3 renormalization 

At this stage a prescription to renormalize tan/3 = f2/wi, 

tan /3 — > tan (3 + S tan /3 , 



(3.46) 



is still called for. Indeed, eq.(3.2(;) given in the previous section was just a formal expression which was 
unrelated to any physical input. There are many possible strategies. The ambiguity is related to the fact 
that tan (3 is just a Lagrangian parameter and as such it is not a physical observable. Its value beyond 
the tree-level is renormalization scheme dependent. (The situation is similar to the definition of the weak 
mixing angle 6w, or equivalently of sin 6w-) However, even within a given scheme, e.g. the on-shell 
renormalization scheme, there are some ambiguities that must be fixed. For example, we may wish to 
define tan/3 in a process-independent ("universal") way as the ratio W2/W1 between the true VEV's after 
renormalization of the Higgs potential |104- 106| , |108| , 112-116|. In this case a consistent choice (i.e. a 
choice capable of renormalizing away the tadpol e con tributions) is to simultaneously shift the VEV's and 
the mass parameters of the Higgs potential, eq.( ^.12 ). 
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(j = 1, 2) in such a way that Svi/vi 
in that scheme -see eq. (3.26)-: 



= Sv2/v2- This choice generates the following counterterm for tan/3 
(5tan/3 



1 



tan/3 



i6ZH2 - SZhi ) . 



(3.48) 



Nevertheless, this procedure looks very formal and one may eventually like to fix the on-shell renormal- 
ization condition on tan /J in a more physical way, i.e. by relating it to some concrete physical observable, 
so that it is the measured value of this observable that is taken as an input rather than the VEV's of the 
Higgs potential. Following this practical attitude, we choose as a physical observable the decay width 
of the charged Higgs boson into r-lepton and associated neutrino: if+ -^ t^Vt. This should be a good 
choice, because: 

1. When Mh± < mt — rnt, the decay iJ+ -^ r^Vr is the dominant decay of H^ already for tan/3 > 2; 

2. From the experimental point of view there is a well-defined method to separate the final state t's 



originating from iJ+-decay from those coming out of the conventional decay W^ 
7?+ -^ t^Vt should be physically accessible; 



so that 



3. At high tan /3, the charged Higgs decay of the top quark can have a sizeable branching ratio, serving 
as a source of charged Higgs particles; and 

4. If Af//± > mt the branching ratio for H^ — > t^Vt never becomes ne glig ible in a wide range of Higgs 
masses to be explored at the LHC rather than at the Tevatron [|60|,|61|. 



28 



MSSM renormalization 



The interaction Lagrangian describing the decay H~^ —f t^Vt is directly proportional to tan/3 

q niT tan /3 „ , , , 

CHru = ;!., H- fPLyr+ h-C. , (3.49) 



\/2M, 



w 



and the relevant decay width is proportional to tan^ (3. 
it reads: 



Whether in the a-scheme or in the GF-scheme, 
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(3.50) 



where we have used the relation (3.1). By measuring this decay width one obtains a physical definition 
of tan/3 which can be used beyond the tree- level. A combined measurement of A/^± and tan/3 from 



charged Higgs decaying into r-lepton in a hadron collider has been described in the literature |117-120| by 
comparing the size of the various signals for charged Higgs boson production, such as the multijet channels 
accompanied by a r-jet or a large missing p^, and the two-r-jet channel. At the upgraded Tevatron, the 
conventional mechanisms gg{qq) — > tt followed by i ^ H^ b have been studied and compared with the 
usual t -^ W~^ b, and the result is that for Mh± — 100 GeV the charged H iggs production is at least as 
large as the W^ production, apart from a gap around tan/3 ~ 6 117| - 119[| (see also chap. p. 

Insofar as the determination of the counter term S tan (3 in our scheme, it can be fixed unambiguously 
from our Lagrangian definition of tan/3 on eq.( 3.49| ) and the renormalization procedure described above 
(and in chapter ^ for the process-dependent terms). It is straightforward to find: 



S tan P Sv 1 
tan/3 V 2 



SZh± + cot /3 5Zhw + ^T ■ 



Notice the appearance of the vacuum counterterm 



which is associated to v'^ — 
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(3.51) 



(3.52) 



eq.(|3.51|). 
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and whose structure is fixed from eg. ( p. 14 ). The last term on 



A, = 



2 ^ 



-SZJ,, - Fr , 



(3.53) 



is the (finite) process-dependent part of the counterterm (see section 4.3). Here Smr/rrir, SZ'^'^ and 
SZJ^ are obtained from eqs.(3.10) and (B.ll) (with m^^ = ); they represent the contribution from the 
mass and wave- function renormalization of the {vt, T)-doublet, including the finite renormalization of the 
neutrino leg. Finally, F^- on eq.(3.53) is the form factor describing the vertex corrections to the amplitude 
of H~^ — > r"*"!/,-; its value can be inferred from the expressions of the vertex functions in chapter ^ by 
substituting the bottom and top quarks (and squarks) masses and couplings by those of r and Ut leptons 

(and sleptons). 

On comparing eqs.( |3.48 ) and ( 3.51 ) we see that the first definition of tan/3 appears as though it is 
free from process-dependent contributions. In practice, however, process-dependent terms are inevitable, 
irrespective of the definition of tan/3. In fact, the definition of tan/3 where 6vi/vi = 6V2/V2 will also 
develop process-dependent contributions, as can be seen by trying to relate the "universal" value of tan /3 
in that scheme with a physical quantity directly read off some physical observable. For instance, if M^o 
is heavy enough, one may define tan/3 as follows: 
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(3.54) 



where we have neglected ml <C ^I%n and 5V stands for the vertex corrections to the decay processes 
A" — > 6 6 and A" — > i i. Since the sum of the mass and wave- function renormalization terms along with the 
vertex corrections is UV-finite, one can consistently choose 5vi/vi — 5v2/v2 leading to eq.(3.48). Hence, 



deriving tan/3 from eq.(3.54) unavoidably incorporates also some process-dependent contributions. 



Any definition of tan j3 is in principle as good as any other; and in spite of the fact that the corrections 
themselves may show some dependence on the choice of the particular definition, the physical observables 
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Figure 3.4: One-loop Feynman diagrams contributing to the charged Higgs self-energy. 



should not depend at all on that choice. However, it can be a practical matter what definition to use in 



a given situation. For example, our definition of tan/3 given on eq.(3.5C) should be most adequate for 
Mf^± < rrit — 771;, and large tan/3, since then H^ — > t^ v^ is the dominant decay of i?^, whereas the 



definition based on eq.(3.54) requires also a large value of tan/3 (to avoid an impractical suppression of 
the hh mode); moreover, in order to be operative, it also requires a much heavier charged Higgs boson, 
since Mu± ~ M^o > 2mt when the decay ^ ^ tt is kinematically open in the MSSM. (Use of light 
quark final states would, of course, be extremely difficult from the practical point of view.) 



3.4.5 Unrenormalized self-energies 

Now we can write down the expressions for the counterterms necessary in the charged Higgs renormal- 
ization. 



The various Feynman diagrams contributing to the charged Higgs self-energy can be seen in Fig. 3.4 
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+ 2M,M°Re {Q'-^.Q'^*) sin/3cos/?B[,j {Mh±,M°, M,) . (3.55) 



Notice that diagram 3.4 (c) gives a vanishing contribution to 5Zu±. The mixed self-energy diagrams are 
in Fig. 3.5 and their contribution read 
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Figure 3.5: One- loop Feynman diagrams contributing to the mixing H^ — W 
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+ M, (sin/? Qf*C^, + cos/? QL*C«) Bi] (Af^±,MO,M,) , 
A sum is understood over all generations. 



(3.56) 



3.5 Sfermion renormalization 



We follow the renormalization procedure with the scalar superpartners of the matter fields. In fact this 
sector is similar to that of the scalar Higgs, for it involves two scalars that mix between themselves, but 
is different because it does not involve tadpole terms, and on the other hand there exist mixing terms 
between the up- type and the down- type sfermions. In the electroweak basis, and for a doublet (t, 6), the 
counterterms needed are 

{5nn,,5Ab,5mi^,5zl) (3.57) 

for sbottom particles, 

{5mu5At,5mi^,6Zy) (3.58) 

for stop particles, and the common counterterms 



{S Mz,SMw,S tan [3,Sfi,Smqj^,SZ'[) . 



(3.59) 



Of course if we would like to perform a supersymmetric renormalization procedure we should use a single 
wave function renormalization constant for fermions and sfermions, thus we should have 
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However, supersymmetry is explicitly broken and so we may take different renormalization constants for 
fermions and sfermions. 



Though the counterterms (3.57), (B.58) and (3.59) are the fundamental blocks of the sfermion sector 

it is more convenient for the on-shell renormalization scheme to use a different set of counterterms. In 

the EW basis we define 
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where the various SM?. are different combinations of the parameter counterterms in (3.57), ( 3.5^ ) 
and ( 3.59|) , except in the case of SMf:_^_^ and 6Mj^ which are related by SU{2) gauge invariance. Thus 
from ( 2.19 ) one can obtain 
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One can also derive the relation between this new set of counterterms and the original from expres- 



sions (2.19) and (2.20): 
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(3.60) 



where the last expression is valid for both type of sfcrmions, just performing the appropriate substitution 



If we had to deal with observables in which all parameters in the RHS of ( 3.60 ) appear in the tree- 
level expressions, as in the squark decays of Higgs particles, then we should invert this set of equations 
to obtain each counterterm corresponding to the appropriate variables. 

In the observables we will compute in this Thesis only one squark appears as an external particle (see 
chapter |6|), thus all the one- loop contributions to other squarks will be higher order contributions to these 
observables. In this situation it is better to use a different approach which uses the physical particles 
themselves. In this approach we introduce mass counterterms for the physical particles, the mixing 
angle, wave function renormalization constants for each squark, and mixing wave function renormalization 
constants, that is 

{8rai^ , Smi^.bdi, 5Z^ ,5Z^ ,5Z'^^ ,5Z'^^) (3.61) 

for the sbottom squark. The number of parameter counterterms in this set is equal to the one in 
the electroweak basis, thus we are not introducing any new parameter counterterm, but using a new 
combination of the old ones. For the wave function counterterms we have added mixing terms between 
the squarks; its purpose is to make possible the one-particle-irreducible (IPI) renormalization procedure 
and avoid the presence of mixing between physical squarks at one-loop. As noted above (section |3.1| ) 
wave function renormalization is unnecessary, but it allows to renormalize the theory by renormalizing 
at the same time every Green function. 



The definition of the renormalization constants (3.61) is 
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bl - {l + -6Z-)b, + 6Z^'bb (ay^b), 
from which we write the one-loop kinetic Lagrangian 

£0 = C + 6C= [d^bl d,b, ~ ml i: 5J (1 + 6Z^) - Sm^ b*, K 

+ {SZ^'^ + 5Z'^^){dn\dJ)2 + d^b*2dJ>i) 

- 5Z^^ ml {bl k + k b;) - SZ^' ml (bl b, + b, b*^) , 

that allows to obtain the one-loop inverse propagator 
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(3.65) 

Next we follow the on-shell prescription requiring the mass parameters to be the physical masses, the 
"residue=l" condition and the absence of mixing between squarks on-shell 

^abimlj = 0, 

(A-,^)' = 1, (3.66) 

and from that obtain the counterterms 

Sml = -E'^'^fm? ) , 

Da Oa 

E'^^fm? ) 

dZ""^ = —^ V ■ (3-67) 
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For fixing SO^, we require that the renormalized mixing angle (that we use as an input data) does 



not feel a shift from the mixed sbottom bare self-energies E" between the physical states ba and b, 



•b 



(a 7^ b). This is similar to the prescription adopted in Refs. B^,p8, 121|, though it is not identical. In our 



formalism, the 3-point Green functions explicitly incorporate the mixed field renormalization constants 
SZ'^'' (a 7^ b) and are therefore renormalized also in the 0^ parameter. The UV-divergent parts of the 
3-point functions are canceled against 56^^ by defining the latter as follows: 

6er = -{5Z'^-SZ^') = - -^ ^-^. (3.68) 

'' 2 ^ ' 2 m? -m,? ^ ^ 

02 bl 

Of course another equivalent choice could just be 66^^ = SZ^"^ (or —6Z^^), but eq.( ^.68| ) is more sym- 
metrical; the numerical differences among the finite parts of the two choices are negligible p7| . This 
renormalization prescription deviates somewhat from the on-shell philosophy, but, contrary to the tan/3 
case, it is not clear by now how the squark angle will ever be measured; thus it is better to use a generic 



criteria, like eq. (3.68) or the ones in Refs. [p3,p8, 121] 



The various QCD Feynman diagrams contributing to the sbottom self-energies can be seen in Fig. 3.6. 
We have also computed the EW contributions in the Yukawa approximation and they are presented in 
chapter o. From these diagrams, and with the interaction Lagrangians of chapter 0, the unrenormalized 
self-energies can be computed. In the following we describe the contributions corresponding to each 
diagram. 

The gluon graph from Fig. |3.6| (a) only contributes to the diagonal self-energy 

E""(fc2)|^^^ = -i4 7ra,Cf (-v4o(™?J + 2Ao(A2) 

+ (2A:2_A2 + 2m?jSo(fc,A,m^j) , (3.69) 

where we introduced a small gluon mass A to regularize the infrared divergence, and Cp — {Nq — 
l)/2 Nc ~ 4/3 is a colour factor. The wave function renormalization constant derived from this expression 
is 

(5Z°|(„)=-i8 7rasCF(2m?^S^ + So)(m^^,A,mgJ , (3.70) 

which is highly simplified if we use the appropriate limits for the various scalar functions B^ of Ap- 
pendix H, obtaining 

^^1(.)=-2|l(A + log^), (3.71) 
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Figure 3.6: Feynman graphs contributing to the sbottom self-energy from the QCD sector. 



where /i is the scale factor and A represents the UV divergence, as defined in equations (A.4) and ( A.16| ) 
respectively. The other gluonic diagram in Fig. 3.6 (b) is zero as it is proportional to 



Ao(A) -^ . 
The purely gluino contribution depicted in Fig. p^ (c) reads 
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which contributes to the wave function renormalization constants as 
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Finally the squark loop contribution from Fig. 3.6 (d) is 
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+ [{4'^? - {R^2^?] Ao{m~j} , 



(3.72) 



(3.73) 



(3.74) 



which only contributes to the mass counterterm and the mixed self-energy (3.67). From the Lagrangian (2.47) 
it could seem that other squarks could give rise to similar contributions, however these are proportional 
to traces of Gell-Mann matrices, and as a consequence they are identically zero. 



Chapter 4 

Quantum effects on t ^ H^ b in the 
MSSM 

4.1 Introduction 

In this chapter we analyze the one-loop quantum effects (both QCD and EW) to the unconventional 
top quark decay mode t -^ _ff + b in the MSSM, and their effect on the Tevatron Collider physicsQ The 
analysis of this process at the quantum level is useful to unravel the potential supersymmetric nature of 
the charged Higgs emerging from that decay. This decay has been subject of interest since very early in 
the literature |122-129|, we wish to emphasize that this decay is not excluded by present data from the 
Tevatron (see sections |4^ and |4.5D . Therefore, we will analyze both the QCD and the EW corrections to 
that decay. The conventional gluon- mediated strong (QCD) corrections have been computed in p^,H], 
the SUSY-QCD corrections mediated by gluinos, stops and sbottoms have been computed in |}48(|, and 
a very detailed discussion can be found in B3]. Here we will concentrate in the remaining part, to wit: 
the electroweak one- loop quantum corrections mediated by squarks, sleptons, charginos, neutralinos and 
supersymmetric Higgs bosons. These corrections were first computed in |5G], and we will combine all of 
these to obtain the full MSSM quantum corrections. 

In this study we will concentrate on those regions of the parameter space in which the partial decay 
width t -^ H^ b is competitive with the SM decay width t — > W~^ b, and where the one-loop EW quantum 
corrections are important. We will take the convention that the decay t -^ H^ b will be interesting 
whenever its branching ratio is BR{t -^ H'^ b) > 10%. Theoretically this condition is taken for granted 
when tan/3 is large enough ( > 30). Under these conditions the SUSY-QCD corrections can be around 
50%, and the purely SUSY-EW ones, induced by the Yukawa couphngs At and A^ can reach 20%, so both 
effects could be measured at the Tevatron and/or the LHC. 



In section 4.2 we present the tree-level relations, and the status of the charged Higgs decay of the top 



quark in view of the Tevatron data. In section 4.3 we describe the process-dependent renormalization 



procedure and we write down all the analytical formulae for the three point irreducible vertices. In 



section |4J we make an exhaustive numerical analysis of the MSSM (QCD-f-SUSY-QCD-f SUSY-EW) 
corrections. In section |4.5| we make an study of the implications of these corrections for the Tevatron 



data. Finally in section 4.6 we present our conclusions. 



4.2 Tree-level relations and experimental determination of 

BR{t -^ H+ b) 



We recall here the tree-level interaction Lagrangian (2.38) in terms of the mass-eigenstates 



Chm = i- *'' H- b [mt cot (3 PR + rrib tun l3PL]t + h.c., (4.1) 

where Vtb is the corresponding Cabibbo-Kobayashi-Maskawa matrix element. On the phenomenological 
side, one should not dismiss the possibility that the bottom-quark Yukawa coupling could play a mo- 

^ The study in the case of generic 2HDM have been presented in [ p2| , see also Rof . [ pl{ for a comprehensive study. 
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mentous role in the physics of the top quark, to the extend of drasticaUy changing standard expectations 
on top-quark observables, particularly on the top-quark width. Of course, this is possible because of the 
potential tan /^-enhancement of that Yukawa coupling. 

From the Lagrangian ( [4.l[ ) , the tree-level width of the unconventional top quark decay into a charged 
Higgs boson reads: 
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(4.2) 



where 



Al/2(1, a;2, y2) ^ ^[l_(2; + y)2][l_(a;-y)2]. (4.3) 

It is useful to compare eq.(|4.2|) with the tree-level width of the canonical top quark decay in the SM: 
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(4.4) 



The ratio between the two partial widths becomes more transparent upon neglecting the kinematical 
bottom mass contributions, while retaining all the Yukawa coupling effects: 
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(4.5) 



We see from it that if A//j± is not much heavier than Af^/, then there are two regimes, namely a low 
and a high tan/3 regime, where the decay rate of the unconventional top quark decay becomes sizeable as 
compared to the conventional decay. They can be defined approximately as follows: i) Low tan/3 regime: 
tan/3 < 2, and ii) High tan/3 regime: tan/3 > mt/mi, ~ 35. The critical regime of the decay t — > H~^ b 
occurs at the intermediate value tan/3 — y^mt/mi, ^ 6, where the partial width has a pronounced dip. 
Around this value, the canonical decay t -^ W^ b is dominant over the charged Higgs decay; more 
specifically, for 3 < tan/3 < 15 the decay rate of the mode t -^ iJ+ b is basically irrelevant as compared 
to the standard mode: BR{t -^ iJ+ b) < 10%. Therefore, a detailed study of the quantum effects within 
that interval is of no practical interest. 

Even though the approximate perturbative regime for tan f3 extends over the wide range 



0.5 < tan/3 < 70, 



(4.6) 



we shall emphasize the results obtained in the phenomenologically interesting high tan (3 region (typically 
tan/3 > 30). As for the low tan/3 range, while BR{t -^ H^ b) can also be sizeable it turns out that 
the corresponding quantum effects are generally much smaller than in the high tan/5 case (Cf. Sec- 
tion 4.4). Still, we find that in the very low tan/3 segment 0.5 < tan/3 < 1 these effects can be of some 
phenomenological interest and we shall also report on them. 

As a matter of fact, and despite naive expectations, the non-SM branching ratio BR{t -^ iJ+ b) is not 
as severely constrained as apparently dictated by the existing measurements of the SM branching ratio 
at the Tevatron, namely, BR{t -^ W^ b) > 70% [130 . To assess this fact, notice that the former result 



strictly applies only under the assumption that the sole source of top quarks in pp collisions is the standard 
Drell-Yan pair production mechanism qq — * it [131|. Now, as noted in chapter y the observed cross- 
section is equal to the Drell-Yan production cross-section convoluted over the parton distributions times 



the squared branching ratio (1.3), whereas in the MSSM one expects a generalization of the production 
mechanism through the production and subsequent decay of i?-odd particles ( |1.4D . 

It should be clear that the observed cross-section on eq. ( |1.4D refers not only to the standard bW bW 
events, but to all kind of final states that can simulate them. Thus, effectively, we should substitute 
BR{t — > W'^ b) in that formula by BR{t — > Xb), and then sum the cross-section over X, where X is 
any state that leads to an observed pattern of leptons and jets similar to those resulting from VF-decay. 
In particular, X = H^ would contribute (see below) to the r-lepton signature, if tan/3 is large enough. 
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Similarly, there can be direct top quark decays into SUSY particles that could mimic the SM decay of 
the top quark [Q. Notwithstanding, even in the absence of the X contributions, eq.( |l.4| ) shows that 
if there are alternative (non-SM) sources of top quarks subsequently decaying into the SM final state, 
W'^ b, one cannot rigorously place any stringent upper bound on BR{t —> W^ b) from the present data. 
The only restriction being an approximate lower bound BR{t — > W^ 6) > 40 — 50% in order to guarantee 
the purported standard top quark events at the Tevatron 0,0]. Thus, from these considerations it is not 
excluded that the non-SM branching ratio of the top quark, BR{t ^"new"), could be as big as the SM 
one, i.e. ~ 50%. 



Notice that at present one cannot exclude eq.(1.4) since the observed form of the conventional t 



W'^ b final state involves missing energy, as it is also the case for the decays comprising supersymmetric 
particles. A first step to improve this situation would be to compute some of the additional top quark 
production cross-sections in the MSSM under given hypotheses on the SUSY spectrum. For instance, the 



inclusion of the qq —>■ g g mechanism followed by the g —>■ tti decay has been considered in Ref. [132.133], 
where it was claimed that BR{t -^ ti Xi) — 50%. By the same token, one cannot place any compelling 
restriction on BR{t -^ iJ+ b) from the present FNAL data. In particular, if tan/3 is large and there exists 
a relatively light chargino with a non- negligible higgsino component, the third mechanism suggested on 
eq.(|l.4|), namely qq^baba followed hyba^'txi, could also be a rather efficient non-SM source of top 
quarks. Moreover, if 100 GeV < Mfj± < 150 GeV, then a sizeable portion of the top quarks will decay 
into a charged Higgs. Thus, if either mt + m^^ < mg < 300 GeV and/or mt + m - < rn^ < 300 GeV, 
so that at least one of the alternative SUSY sources of top quark final states contributing to eq.( |l.4|) 
is available (and mg , TOj, are not too heavy so that the production cross-section is not too phase-space 
suppressed) , then one may equally argue that a large branching ratio BR{t -^ H^ b) ~ 50% is not 
incompatible with the present measurement of the top quark cross-sections y, g . This could be most 
likely the case if the frequently advocated SUSY decay i ^ Xi ^i is kinematically forbidden. Nonetheless, 
even if it is allowed, it is non-enhanced in our preferential large tan/3 region, in contrast to i — > H'^ b. 

Furthermore, it is worth mentioning that the decay mode t — *■ H^ b has a distinctive signature which 
could greatly help in its detection, viz. the fact that at large tan/3 the emergent charged Higgs would 
seldom decay into a pair of quark jets, but rather into a r-lepton and associated neutrino. This follows 
from inspecting the ratio 

V{H+-^T+Vr) _ I (mr\^ tan2/3 

T{H+ -^ cs) ~ 3 Kj^cJ (m2/m2)tan2/3-f cot2/3 

2 



1 /Tot 



-^ o^ >10 (for tan/3> ^m.c/ms > 2) , (4.7) 

where we see that the identification of the charged Higgs decay of the top quark could be a matter of 
measuring a departure from the universality prediction for all lepton channels. In practice, r-identification 
is possible at the Tevatron; and the feasibility of tagging the excess of events with one isolated r-lepton 
as compared to events with an additional lepton has also been substantiated by studies of the LHC 
collaborations |l2Cf| . The experimental signature for tt — > H^ H~ b b would differ from tt — > W~^ W~ b b 
by an excess of final states with two T-leptons and two b-quarks and large missing transverse energy. 



A study in this direction by the CDF collaboration at the Tevatron |134] has been able to exclude 
a large portion of the (tan/3,MjLf±)-plane characterized by tan/3 > 60 and M^± below a given value 
which varies with tan/3. For extremely high tan/3 > 0(100), the uppermost excluded mass region is 
M[j± < 140 GeV. However, within the interval tan/3 = 60 — 80, the allowed upper limit on Mh± varies 
very fast with tan/3. In particular, the MSSM permissible values Mh± > 110 GeV (compatible with 
Mao > 75 GeV) are not manifestly excluded for tan/3 equal or below the perturbative bound tan/3 = 70, 



eq.(4.6). On the other hand radiative corrections alter this bounds in a significant way (see sec. 4.5) 
We shall nevertheless err on the conservative side and assume that tan/3 < 60 throughout our analysis. 
Thus, as far as the high tan/3 regime is concerned, we will for definiteness optimize our results in the 
safe, and phenomenologically interesting, high tan /3 segment 

30<tan/3<60. (4.8) 

To round off the r-lepton business, it has been shown that it should be fairly easy to discriminate 
between the M^-daughter r's and the iJ^-daughter r's by just taking advantage of the opposite states of 
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Figure 4.1: Tree-level Feynman diagram of the unconventional top quark decay t -^ H^ b. 



T polarization resulting from the W^ and H^ decays; the two polarization states can be distinguished 
by measuring the charged and neutral contributions to the 1-prong r-jet energy (even without identifying 
the individual meson states) | 134 136]. 

In short, there are good prospects for detecting the decay t -^ iJ+ b, if it is kincmatically accessible. 
Unfortunately, on the sole basis of computing tree-level effects we cannot find out whether the charged 
Higgs emerging from that decay is supersymmetric or not. Quantum effects, however, can. 



4.3 One-loop corrected r{t -^ H+ b) in the MSSM 

First of all let us sketch the renormalization procedure for the tbH^-vertex. The full tbH^ bare 
Lagrangian is found by substituting the t, b and H^ fields by bare fields, and g and the various masses 
by bare parameter s in (4.1 ), as defined in chapter ^ There is still another piece to be add, name ly th e 
W^ — H^ mixing ( |3.24 ), which must substitute the W^ field in the gauge interaction Lagrangian ( ^.42 ). 
Proceeding in this way we find 



C%t,t = Cmt + -^ H- b [SCr mt cot /? Pr + SCl mt tan/3 Pl] t + h, 

V 2 Mw 



(4.9) 



with 



SCr 

SCl 



Smt Sv 1 



mt 
Smiy 
nib 



V 

Sv 



SZ 



H+ 



SZ 



H+ 



szi 



5Z\ 



SZ'r 



SZ", 



S tan [3 
tan/3 

5tan/3 
tan/3 



+ SZhw tan/3, 
— SZhw cot /3 , 



(4.10) 



0.999 within ±0.1%, from unitarity of the CKM-matrix under the 



and where we have set Vtb — 1 {Vtb 
assumption of three generations). 

As stated in Section 4.2, the study of the decay t -^ H~^ b is worthwhile in the small (tan/3 < 2), and 
most conspicuously in the high (tan/3 > 30) tan/3 region, where the branching ratio can be comparable 
to the one of the standard decay t -^ W^ b. These are, therefore, the regions on which wc will focus our 
search for potentially significant (strong and electroweak like) SUSY quantum corrections to t ^ H^ b. 
As for the strong effects, they can be rather large and have been evaluated in Refs. |^,Q; here we shall 
not dwell any longer on their detailed structure apart from including them in our numerical analysis and 
recalling some interesting remarks in section 4.4. 

On the electroweak side, one may also expect sizeable quantum corrections from enhanced Yukawa 
couplings of the type (1.2). In the relevant tan/3 regions mentioned above, the latters yield the lead- 
ing electroweak contributions and in these conditions we will neglect the pure gauge corrections from 
transversal gauge bosons in the Feynman gauge. Moreover, as already stressed in section 4.2, the branch- 
ing ratio of the charged Higgs mode in the intermediate tan f3 region is too small to speak of, so that the 
detailed structure of the radiative corrections in this range is irrelevant. 

In the following we will describe the relevant electroweak one-loop supersymmetric diagrams entering 
the amplitude of t ^ i/+ b in the MSSM. At the tree-level, the only Feynman diagram is the one in 
Fig. 4.1. At the one-loop level, we have the vertex diagrams exhibited in Figs. 



3 and the fermion 



and Higgs self-energies of chapter ^. The computation of the one-loop diagrams requires to use the full 
structure of the MSSM Lagrangian. 
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Figure 4.2: One-loop SUSY-EW Feynman diagrams contributing to the decay t -^ i7+ b. 



Specifically, Fig. 4.2 shows the electroweak SUSY vertices involving squarks, charginos and neutralinos. 



In all these diagrams a sum over all indices is taken for granted. The supersymmetric Higgs particles of 
the MSS M a nd Goldstone bosons (in the Feynman gauge) contribute a host of one-loop vertices as well 
(see Fig. |4.3| ). As for the various self-energies, they will be treated as countcrtcrms to the vertices. Their 
structure is dictated by the Lagrangian ( [4.9| ). 

Although we have displayed only the process dependent diagrams, the full analysis should also include 
the SUSY and Higgs/Goldstone boson contributions to the various universal vacuum polarization effects 
comprised in our counterterms. However, the calculation of all these pieces has already been discussed 
in detail long ago in the literature [pq, po|-p2|, M and thus the lengthy formulae accounting for these 
results will not be explicitly quoted here. Their contribution is not tan/3-enhanced, but since we wish to 
compute the full supersymmetric contribution in the relevant regions of the MSSM parameter space, those 
effects will be included in our numerical code. Finally, the smaller -though numerically overwhelming 
- subset of strong supersymmetric one-loop graphs are displayed in Fig.^^ of Ref. Q. We will use 
the formulae from the latter reference in the present analysis to produce the total (electroweak -|-strong) 
SUSY correction to our process. 

Next let us report on the contributions from the various vertex diagrams and counterterms in the on 
shell renormalization scheme. The generic structure of any renormalized vertex function. A, in Figs. 4.2-4.3 
is composed of two form factors Fl, Fh plus the counterterms. Therefore, on making use of (4.9) and 
the formulae of chapter 0, one immediately finds: 



A = 



where 



«.9 



V2Mw 



[mt cotp(l + Aii)PR + mbtanP{l + AL)PL] , 



(4.11) 



mt 2 A 

- -^r- +6Zh+ + {ta.n l3 - cot 13) SZhw , 

Al = F, + ^-^ + \5Zl + \5Zl + K. 
mi, 2 2 



(4.12) 



In the following the analytical contributions to the vertex form factors and counterterms will be specified 
diagram by diagram. 
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4.3.1 SUSY vertex diagrams 

In this section we will make intensive use of the definitions and formulae of chapter g. We refer the reader 
there for questions about notation and conventions. Following the labeling of Feynman graphs in Fig. 4.2 
we write down the terms coming from virtual SUSY particles. 



• Diagram (a): Making use of the coupling matrices of eqs. (2.36) and (g) we introduce the short- 
handsQ 

(4.13) 



A'l. and Af 



A 



(*) 



and define the combinations (omitting indices also for QaiiQai) 

A(i) = cos/3^;g^AL°^ , £;(!) = cos13A*_Q^A^°^ , 

b(i' = COS 13 A\_Q^Af , F(i) = cos pA*_Q^Af , 

C(i) = sin/3A;Q^AL°^ , G(i) = sin/3^1Q^AL°^ , 

D^^^ = sin/3A;Q^vl^°^ , H^^^ = sin [3 A*_Q^ A'^+^ . 



(4.14) 



The contribution from diagram (a) to the form factors Fl and Fn is then 

Fl = Ml [h^^^Co + 

+ mt (mtH^^^ + M°G(i) + mi,A^^^ + M,£;(i)) (Cn - C12) 

+ fmtTObA^i) + mtM^E^^'i + M^mtB'-^'^ + M,M°F'-^A Cq 

Fr = Mfl. [a(i)C'o+ 

+ TOt (mtA^i) + Af^B(i) + mbi/fi) + AhD'^^^^ (Gn - C12) 

+ (mtmbH^^^ + TJitA'hD^^^ + M^mtG^^^ + M,M°C^^A Co 



where the overall coefficients Ml and Mr are the following: 

ig'^Mw ,, ig^M^ 



Mr = — 



mi, tan /3 



Affl, 



Jw 



mt cot /3 



(4.15) 



(4.16) 



The notation for the various 3-point functions is summarized in Appendix W. On eq. (4.15) they 
must be evaluated with arguments: 



a^a{p,p',mi^,Mli,M,) 



(4.17) 



• Diagra m (b): For t his d iagram -which in contrast to the others is finite- we also use the matrices 
on eqs. ( 2.36 ) and ( 2.40 ), and introduce the shorthands 



A 



(fc) 



Ail and A^) 



A 



(t) 

ztaa. 1 



(4.18) 



to define the products of coupling matrices 



A(^) = CaA^l^A^^ , C(2) = G,aA^^>A^^ , 



i?(2) = G,,4^>4) , 



Z?(2) = G,aAL'^*A^) 



(4.19) 



^ Lower indices are summed over, whereas upper indices (some of them within parenthesis) are just for notational 
convenience. 
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The contribution to the form factors Fl and Fj^ from this diagram is 

Ml 



Fl - 
Fr = 



2Mw 
Mr 

2Mw 



rmB^^^Cu + mtC^^^ {Cn - C12) - M^D^^^Co 
mbC^^^Ci2 + mtB^^^ {Cn ~ C12) ~ M^A^^'^Co 



(4.20) 



the coefficients M^, Mr being those of eq. (4.16) and the scalar 3-point functions now evaluated 
with arguments 

C* = C* [p, p', M°, mi^ ,m.i^j . 



• Diagram (c): For this diagram we will need 

l(6) 



A± ^ A^± and Al°) ^ Ai'^ 



(4.21) 
(4.22) 



± — ■^±aa ' 

and again omitting indices we shall use 

A(3) = cos/3Af *g^A_ , £;(3) = cos/3AL°)*g^A_ , 

5(3) = cos PA^°'>*Q^A+ , f(3) = cos (3 A^°'>*Q^A+ , 

C(3) = sin/3Af *Q^A_ , G^^) = sin/3AL°^*Q^A_ , 

i?(3) ^ sin/34°)*Q«A+ , i/(3) ^ sin/3v4L°^*Q^A+ . 

From these definitions the contribution of diagram (c) to the form factors can be obtained by 
performing the following changes in that of diagram (a), eq. (4.15): 

— Everywhere on eqs. (4.15) and (4-17) replace Mi <-> M° and m^ ^-> m^ . 



(4.23) 



— Replace on eq. (4-15) couplings from (4.14) with those of ( 4.23| ). 

— Include a global minus sign. 

4.3.2 Higgs vertex diagrams 

Now we consider contributions arising from the exchange of virtual Higgs particles and Goldstone bosons 
in the Feynman gauge, as shown in Fig. 4.3. We follow the vertex formula for the form factors by the 
value of the overall coefficient N and by the arguments of the corresponding 3-point functions. 



Diagram (a): 



Fl = N[ml{Ci2-Ca) + m^tC0t^PiCii-Ci2)] 
Fr = Nml[Ci2-Co+t&n'l3{Cn-Ci2)], 



N 



T- 



.W 



1 



{Mfja,Mf^o}\ {cos a, sin a} 



2 V 2M^ J cos/3 

a = a {p, p, mb, Mh± , {Mho , A^o}) . 

• Diagram (b): 

Fl - Ncotp[mUCii-Ci2)+mliCo-Ci2)], 
Fr = NmltanP{2Ci2-Cn~Co), 
ig^ {cos a, sin a} 



{cos(/3 — a), sin(/3 — a)} , 



N 
• Diagram (c): 



-{sin(/3-a),cos(/3-a)} 
4 cosp 

= C** {P,p', mb, Mw, {MHO,Mho}) . 



M^± {M]^o,Mlo} 



M^ 



M^ 



Fl = 7Vm?[cot2/3Ci2+Cn-Ci2-Co], 

Fr = N[mltan'pCi2+m^t{Cn-Ci2~Co)]. 



N = 
C, = 



ig^ {sin a, cos a} 



2 sin/3 

C,{p,p',mt,{ Mho , Mho } , Mh 



{cos(/3 — a), sin(/3 — a)} 1 1 



{M],„,Ml„} 



2M2, 
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Figure 4.3: One-loop Higgs Feynman diagrams contributing to the decay t -^ H^ b. 
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• Diagram (d): 

Fl 
Fr 

N 

• Diagram (e): 



Diagram (f): 



Nmt{2Ci2 - Cii + Co) cot (3 , 
N i-mlCu + m2(Cii - C12 - Co)] tan/3 , 
ig^ {sin a, cos a} 



T 



4 sin/3 

a (p, p' , mt , { Mho , A4o } , Mi^ 



{sin(/3 — a), cos(/3 — a)} 



Af|± {M^o,M2o} 



A^M^ 



A/2 



Fl 
Fr 

N 

Fr 

N 



iV [m2(Ci2 + Co 
Nml tan2/3(Cii 



' 4 



Af^+ A/2 



-m2(Cii 
Co), 

A" 



A^W 



A^W 



= C* (p,p', mfc, Afw, Ma«) . 



Nmlcoi^l3{Cii+Co), 
7V[mgCi2+m2(Cii-Ci2 



4 









C,ip,p',m,t,MAo,Mw) 



C12)] 



Co)] 



• 


Diagram (g) 






Fl 






Fr 






N 






C, 


• 


Diagram (h) 



N [{2mlCn + Co + 2{mj - ml){Cn - C12)) cot^f3 + 2ml{Cn + 2Co)]m^t , 

N [{2mlCn + Co + ^m^t - ^DiCn ~ C^)) tan2/3 + 2m?(Cn + 2Co)]m2 , 
^2 



± 



tg- 



sm a cos a 



4Af2^ sin /3 cos /3 ' 
C* (p,p',{Af^o,Af,jo},mt,mf,) . 



N 



Nmf cot2/3 Co , 
Nml tan2/3 Co , 



= T 



«r 



4Af2, 



C* {p,p',{MAO,Mz},mt,mb) 



In the equations above, it is understood that the CP-even mixing angle, a, is renormahzed into acs 
by the one- loop Higgs mass relations 1 112 -116]. 

The evaluation of A^-, the process dependent term of the tan/3 renormalization, on eq.( |3.53 ), yields 
similar bulky analytical formulae, which follow after computing diagrams akin to those in Figs. 4.2-4.3 
and the corresponding counterterms to the r, v^ and iJ+ external legs from chapter y for the MSSM 
corrections to H^ —^ t^ Vt- We refrain from quoting them explicitly. The numerical effect, though, will 



be explicitly used, and isolated (Fig. 4.17 ), in our computation. 

We are now ready to furnish the corrected width of t ^ if"*" h in the MSSM. It just follows after 
computing the interference between the tree-level amplitude and the one-loop amplitude. It is convenient 
to express the result as a relative correction with respect to the tree-level width both in the a-scheme 
and in the GF-scheme. In the former we obtain the relative MSSM correction 



srMSSM 



r(o) 

J- a 



r(o) 

J- a 



Nl 
D 



[2Re{AL)] 



Nr 
D 



[2Re{AR)] 



Nlr 
D 



2Re{Ai 



Ar 



(4.24) 
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where the corresponding lowest-order width is 






with 



D 



16 M^ mt 



^ (.-^'Z:2'~i:2~ 



mf 



(4.25) 



D 

Nl 

Nlr 



= (mj +ml- Mfj± ) (m^ cot^ (3 + ml tan^ (3) + ^mlml , 
= [ni^ + ml — Mff± ) mf, tan /? , 
= (mj +ml- Mfj± ) mj cot^ /? , 



(4.26) 



From these equations it is obvious that at low tan (3 the relevant quantum effects basically come from the 
contributions to the form factor A^ whereas at high tan/3 they come from A^. 

Using eq. (|3.l[) we find that the relative MSSM correction in the GF-parametrization reads 



cMSSM 
Gf 



r-r 



.(0) 



(0) 

Gf _ gMSSM _ ^^MSSM 



(4.27) 



.(0) 



where the tree- level width in the G^-scheme, F^ , is given by eq.(F2) and is related to eq.( 4.25 ) through 



F(«) = f[?^ (1 - a 



MSSM\ 



(4.28) 



4.4 Numerical analysis and discussion 

Quantum effects should be able to discriminate whether the charged Higgs emerging from the decay 
t —f H^ b is supersymmetric or not, for the MSSM provides a well defined prediction of the size of these 
effects for given values of the sparticle masses. Some work on radiative corrections to the decay width of 
t —f H'^ b has already appeared in the literature. In particular, the conventional QCD corrections have 
been evaluated |E3,E4] and found to significantly reduce the partial width. The SUSY-QCD corrections 



are also substantial and have been analyzed, only in part in Refs. [137, 138 1, and in more detail in 
Refs. p2,Eq]. The clcctroweak corrections produced by the roster of genuine (i?-odd) sparticles was first 
computed in Rc f. pO| . As for the virtual effects mediated by the Higgs bosons, a first treatment is given in 
Refs. [139| and ]140| |. However, these references disagree in several parts of the calculation, and moreover 
they are both incomplete calculations on their own, for they fully ignore the Higgs effects associated 
to the bottom quark Yukawa coupling, which could in principle be significant in the large tan/3 region. 
On the other hand, even though the latter kind of Higgs effects have been discussed in the literatur e 
in other renormalization schemes based on alternative definitions of tan/3 [104- 106 , 108| , 115, 141- 145 [, 
a detailed analysis including the genuine SUSY effects themselves has never been attempted. Thus, if 
only for completeness, we are providing here not only a dedicated treatment of the i?-odd contributions 
mediated by the sparticles of the MSSM, but also the fully-fledged pay-off of the supersymmetric Higgs 
effects. 

Before presenting the results of the complete numerical analysis, it should be clear that the bulk of the 
high tan/3 corrections to the decay rate of i ^ H^ b in the MSSM is expected to come from SUSY-QCD. 



This could already be foreseen from what is known in SUSY GUT models [146-148[; in fact, in this 



context a non-vanishing sbottom mixing (which we also assume in our analysis) may lead to important 
SUSY-QCD quantum effects on the bottom mass, rrih = m'^^^ + Ami,, where Ami, is proportional to 
M^j^ — > — /itan/3 at sufficiently high tan/3. These are finite threshold effects that one has to include when 
matching the SM and MSSM renormalization group equations (RGE) at the effective supersymmetric 
threshol d sca le, Tsusy, above which the RGE evolve according to the MSSM /3- functions in the MS 
scheme 149 . In our case, since the bottom mass is an input para meter for the on-shell scheme, these 



effects are just fed into the mass countcrtcrm Snib/mb on eq.(4.12) and contribute to it with opposite 
sign {Smt/mt = —Amb/mt, -I- . . .). 

Explicitly, when viewed in terms of diagrams of the electroweak-eigenstate basis, the relevant finite 
corrections from the bottom mass counterterm are generated by mixed LR-sbottoms and gluino loops 
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Figure 4.4: Finite Feynman diagrams contributing to the renormalization of the bottom mass in the EW 
basis. 
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— mg/i tan /i;/(m^^,m^^,TOg), 



(4.29) 



where the last result holds for sufficiently large tan/3 and for /i not too small as compared to Ai,. 
have introduced the positive-definite function (Cf. Appendix Kf) 



We 



/(mi,TO2,m3) = 16 7r^iCo(0,0, mi, 7712,7713) 



■In. 



it z J 



+ m\ rn^ In ■ 



(77I1 — 7772) ("^2 ~ "^3) (™i ~ "^3) 



(4.30) 



In addition, we could also foresee potentially large (finite) SUSY electroweak effects from dmb/mb- They 
are induced by tan /3-enhanced Yukawa couplings of the type (1.2). Of course, these effects have already 
been fully included in the calculation presented in Section 4.3 that we have performed in the mass- 



eigenstate basis, but it is illustrative of the origin of the leading contributions to pick them up again 
directly from the diagrams in the electroweak-eigenstate basis. In this case, from loops involving mixed 
LR-stops and mixed charged higgsinos (Cf. Fig.|4.4|(b)), one finds: 



(5771b 
'^'' / SUSY-Yukawa 



■ Y^ — "^* M^i^Iimi^ , 777^^ , /i) 



167r^ 



^ tan (3 At I{mi^ , mi^ , /i) , 



(4.31) 



where again the last expression holds for large enough tan [3. 



Notice that, at variance with eq.( 4.2£ ), the Yukawa coupling correction ( 4.31 ) dies away with increasing 
/z. Setting ftt ~ 1 at high tan /?, and assuming that there is no large hierarchy between the sparticle masses, 
the ratio between (4.29) and (4.31) is given, in good approximation, by Amg/At times a slowly varying 
function of the masses of order 1, where the (approximate) proportionality to the gluino mass reflects the 
very slow decoupling rate of the latter E2,Ea]. 



In view of the present bounds on the gluino mass [150|, and since At (as well as Ab) cannot increase 
arbitrarily, we expect that the SUSY-QCD effects can be dominant and even overwhelming for sufficiently 
heavy gluinos. Unfortunately, in contradistinction to the SUSY-QCD case, there are also plenty of 
additional vertex contributions both from the Higgs sector and from the stop-sbottom/gaugino-higgsino 
sector where those Yukawa couplings enter once again the game. So if one wishes to trace the origin of 
the leading contributions in the electroweak-eigenstate basis, a similar though somewhat more involved 
exercise has to be carried out also for vertex functions. Of course, all of these effects are automatically 
included in our calculation of section |4.3| within the framework of the mass-eigenstate basisn. It is worth 
noticing that these effects do not decouple when we send the scale of the sparticle masses to infinity. 

•^Thc inass-cigenstate basis is extremely convenient to carry out the numerical analysis, but it does not immediately 
provide a "physical interpretation" of the results. The electroweak-eigenstate basis, in contrast, is a better bookkeeping 
device to trace the origin of the most relevant effects, but as a drawback the intricacies of the full analytical calculation can 
be (in general) abhorrent. 
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Indeed, if we scale all SUSY parameters in cqs. (4.29) and ( [1.31 ) by the same factor, this factor drops off 
from these formulas; and of course this factor could be sent to infinity! Therefore, we have here a dramatic 
example of non-decoupling phenomenon which is intrinsically associated to the breaking of SUSY. Indeed, 
in a truly supersymmetric theory these threshold corrections would (obviously) be zero! Another matter 
would be to judge what sort of weird fine-tunings among the parameters would entail letting the scale 
factor to infinity, but what cannot be denied is the bare fact that it really drops out automatically! One 
could work in an alternative framework (as in Ref. p8[ ) assuming no mixing in the sbottom and/or stop 
mass matrix. Nonetheless, the typical size of the radiative corrections does not change as compared to the 
present approach (in which we do assume a non-diagonal sbottom and stop matrix) the reason being that 



{b,t} 

LR 



0, the contribution 6mi,/mi, ex — /xtan/? at large tan/3 is no longer 



in the absence of mixing, i.e. M 
possible but, in contrast, the vertex correction (Cf. diagram [4.2|(b) and its SUSY-QCD analog) does 
precisely inherits this dependence and compensates for it (see eq. ( p.20[ )). The drawback of an scenario 
based on Af]^^ = 0, however, is that when it is combined with a large value of tan ( 3 it may lead to a 
value of Ab which overshoots the natural range expected for this parameter -eq. ( 2.24 ) . 

We may now pass on to the numerical analysis of the over-all quantum effects. After explicit compu- 
tation of the various loop diagrams, the results are conveniently cast in terms of the relative correction 
with respect to the tree-level width: 



^H - T^H _ r(i ^H+b)- r(o) {t -^ H+ b) 



.(0) 
H 



T(o){t-.H+b) 



(4.32) 



In what follows we understand that 6 defined by eq.(4.32) is 6a ~ Cf. eq.(4.24) - i.e. we shall always give 
our corrections with respect to the tree-level width Fj^ in the a-scheme. The corresponding correction 



with respect to the tree-level width in the GF-scheme is simply given by eq.( 4.27 ), where A' 



^MSSM 



was 



object of a particular study [g6|,g9|,^ and therefore it can be easily incorporated, if necessary. Notice, 
is already tight ly bound by the experimental data on Mz at LEP and the ratio 



however, that Ar^^'^'^*^ 

Mw/Mz in pp, which lead to Mw ( ^.48]) . Therefore, even without doing the exact theoretical calculation 

of Ar within the MSSM, we already know from 



Ar = 1 



1 



V2Gf M2, (1 - Af2,/M| 



(4.33) 



that Ar^^"^"^*^ must lie in the experimental interval Ar 



0.037 ±0.014. 



larger than A 



MSSM 



Now, since the corrections computed in Section 4.3 can typically be about one order of magnitude 



the bulk of the quantum effects on t — )■ iJ+ b is already comprised in the relative 
correction (4.32) in the a-schemerl Furthermore, in the conditions under study, only a small fraction 
of Ar*^'^'^*^ is supersymmetric pfl,p9, 93 , and we should not be dependent on isolating this universal, 
relatively small, part of the total SUSY correction to 6. To put in a nutshell: if there is to be any hope 
to measure supersymmetric quantum effects on the charged Higgs decay of the top quark, they should 
better come from the potentially large, non-oblique, corrections computed in Section 4.3. The SUSY 
effects contained in l^r^^'^ssM ^ 



M,p3], instead, will be measured in a much more efficient way from a 
high precision {5M^^ = ±40 MeV) determination of Mw at LEPII. 



Another useful quantity is the branching ratio 



Bh = BR{t -^ H+b) = 



F 



H 



^W + T/f + VsuSY 



(4.34) 



where Tw = r{t -^ W^ b) and Vsusy stands for decays of the top quark into SUSY particles. In 
particular, the potentially important SUSY-QCD mode t ^^ tig is kinematically forbidden in most 
part of our analysis where we usually assume nig — 0(300) GeV. There may also be the competing 
electroweak SUSY decays t ^> iix^a ^-^d t ^ bixt ^'^^ some a = 1, . . . , 4 and some i = 1,2. The latter, 
however, is also phase space obstructed in most of our explored parameter space, since we typically assume 
TTijj — 150 GeV. The decay t ^>- ti Xa, instead, is almost always open, but it is not tan /3-enhanced in our 



"^For the standard decay t 
the maximal SUSY electroweak corrections 



W^ b, the situation is quite different since the SM electroweak corrections 151- 158 1 and 
:0| in the Q-scheme are much smaller than for the decay t — > H^ b , namely 
they are of the order of Ar. Therefore, for the standard decay t — > VK+ b there is a significant cancellation between the 
corrections in the a-scheme and Ar in most of the tan /3 range resulting in a substantially diminished correction in the 
Gj?-scheme. 
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favourite segment ( |4.8| ). However, when studying the branching ratio (4.34) as a function of the squark 
and gluino masses, we do include the effects from ah these supersymmetric channels whenever they are 
kinematically open. Thus in general Tsusy on eq.( |4.34 ) is given by 



TsusY = T{t ^ i, ~9) + J2 r(^ ^ *i Xa) + E r(* ^ ^1 X»^) ■ 



(4.35) 



The various terms contributing to this equation are computed at the tree-level. Recently, the SUSY- 
QCD corrections to some of these supersymmetric modes have been evaluated and in some cases may 
be important |159, 160 1. Similarly, we treat the computation of the partial width of the standard mode 
t — > W^ b at the tree-level. This is justified since, as shown in Refs. pO|-^, 161-163{, this decay cannot 



in general develop large supersymmetric radiative corrections, or at least as large as to be comparable to 
those affecting the charged Higgs mode (for the same value of the input parameters). The reason for it 
stems from the very different structure of the counterterms for both decays; in particular, the standard 
decay mode of the top quark does not involve the mass renormalization counterterms for the external 
fermion lines, and as a consequence the aforementioned large quantum effects associated to the bottom 
quark self-energy at high tan/3 are not possible. 

Figures 4.5-4.17 and 4.19- 4.2C| display a clear-cut resume of our numerical results. We wish to point 
out that they have been thoroughly checked. Scale independence of 6, eq.(4.32), and cancellation of 
UV-divergences have been explicitly verified, both analytically and numerically. In all our numerical 



evaluations we have imposed the various restrictions of sec. 2.6, 

To start with, we co ncentrate on the case /i < 0, which we study in Figs. 4.5- 4.17| . (The case fi > 



is studied apart in Figs. 4.1E - 4.20 and will be commented later on.) We observe that, for negative fi, the 
leading SUSY-QCD effects on S are positive. This means that in these circumstances the potentially large 
strong supersymmetric effects are in frank competition with the conventional QCD corrections, which 
are also very large and stay always negative as will be discussed later on. 

Needless to say, a crucial parameter to be investigated is tan/3. In Fig. 4.5 we plot the tree-level width, 
ro(i — > H~^b), and the total partial width, TMSSAiit -^ H'^b), comprising all the MSSM effects, as a 
function of tan /3. A typical set of parameters is chosen wel l w ithi n can onical expectations (see below); the 
individual influence of each one of them is tested in Figs. 4/7 to 4.17 , Also shown in Fig. ^^ is the (tree- 
level) partial width of the standard top quark decay t -^ W^ b, whi ch i s (as noted above) far less sensitive 
to quantum corrections. For convenience, we have included in Fig. 4^ a plot of rQCD(^ —^ H^ b), i.e. the 
partial width that would be obtained in the presence of only the standard QCD corrections. Recently 
it has been shown that EW corrections in the absence of SUSY particles, i.e. 2IIDM effects, can also be 
large, though it happens for a Hi ggs spectrum very different of the MSSM one |6|,^. On the other 
hand by simple inspection of Fig. 4.6 we see that the Higgs sector corrections are indeed small, so the 
QCD corrected curve is in fact the partial width that would be obtained in a non-SUSY 2HDM under 
the conditions of the present study. 



From eq.(4.24) it is clear that, for large (resp. small) tan/3, the renormalized form factor yielding the 
bulk of the SUSY contribution is A]^ (resp. Aj i). To app raise the relative importance of the various 
types of MSSM effects on r(i -^ H^ b), in Figs. |4.6| (a) -[4.6|( b) we provide plots for the correction to the 



partial width, eq.( 4.32| ), and to the branching ratio, eq.( 4.34 ), as a function of tan /3, reflecting the various 
individual contributions. Specifically, we show in Fig. |4.6| (a): 



• (i) The supersymmetric electroweak contribution from genuine (i?-odd) sparticles (denoted 5susy- 
i.e. from sfermions (squarks and sleptons), charginos and neutralinos; 



EW) 



(ii) The electroweak contribution from non-supersymmetric (i?-even) particles {Sew)- It is com- 
posed of two distinct types of effects, namely, those from Higgs and Goldstone bosons (collectively 
called "Higgs" contribution, and denoted ^Higgs) plus the leading SM effects |^,^,|9^ from con- 
ventional fermions ((5sm)- 

Sew = feiggs + SsM ; (4.36) 

The remaining non-supersymmetric electroweak effects are subleading and are neglected. 

(iii) The strong supersymmetric contribution (denoted by ^susy-qcd) from squarks and gluinos; 

(iv) The strong contribution from conventional quarks and gluons (labeled (5qcd); and 
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Figure 4.5: The total partial width, TMSSM{t -^ H^b), including all MSSM effects, versus tan/3, as 
compared to the tree-level width and the QCD-corrected width. Also plotted is the tree-level partial 
width of the standard top quark decay, t -^ W^ b. The masses of the top and bottom quarks are 
mt — 175 GeV and nib ~ 5 GeV, respectively, and the rest of the inputs are explicitly given. We remark 
that At = Ab = ... = A is a common value of the trilinear coupling for all squark and slepton generations. 
Unless explicitly stated otherwise, the inputs staying at fixed values in the remaining figures are common 
to the values stated here. 
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Figure 4.6: (a) The relative corrections S, eq.( 4.32| ), as a function of tan/3. Shown are the SUSY-EW, 
standard EW (i.e^ non-supersymmetric electrowe ak) , SUSY-QCD, standard QCD, and total MSSM 



contribution, eg. ( 4. 37 ); (b) The branching ratio ( 4.34 ), as a function of tan/3; separately shown are 
the values of this observable after including standard QCD corrections, full MSSM corrections, and the 
tree- level value. 



• (v) The total MSSM contribution, (5mssm, namely, the net sum of all the previous contributions: 

^MSSM — (^SUSY-EW + ^EW + ^SUSY-QCD + ^QCD- (4.37) 



In Fig. 4.6(b) we reflect the impact of the MSSM on the branching ratio, as a function of tan/3; also 
shown are the tree-level value of the branching ratio and the latter quantity after including the (non- 
supersymmetric) QCD corrections. A typical common set of inputs has been chosen in Figs. 4.6(a)-E^b) 
such that the supersymmetric electroweak corrections reinforce the strong supersymmetric effects (SUSY- 
QCD). For this set of inputs, the total MSSM correction to the partial width of t ^ iJ+ b is positive 
for tan/3 > 20 (approx.). Remarkably enough, this is so in spite of the huge negative effects induced 
by QCD. In fact, we see that the gluon effects are overridden by the gluino effects provided tan/3 is 
sufficiently large, to be concrete for tan/3 > 30. Beyond this value, the strength of the supersymmetric 
loops becomes rapidly overwhelming; e.g. at the representative value tan/3 = mt/mf, = 35 we find 
^MSSM — -1-27%; an d at tan /? ~ 50, which is the preferred value claimed by 5*0(10) Yukawa coupling 
unification models [146-148|, the correction is already (5mssm — +55%. Quite in contrast, at that tan/3 
one would expect, in the absence of SUSY effects, a (QCD) correction of about —57%, i.e. virtually of 
the same size but opposite in sign! 

Coming back to Fig. 4.5, we see that, after including the SUSY effects, the partial width of i ^ _ff+ b 
equals the partial width of the standard decay t -^ W^ b near the "5*0(10)" point tan/3 = 50. (The 
meeting point is actually a bit earlier in tan/3, after taking into account the known [BO,Blt|, negative, 
SUSY corrections to t -^ W^ b, but this effect is not shown in the figure since it is relatively small.) 
Now, for the typical set of parameter values introduced in Fig. 4.5, the top quark decay width into SUSY 



particles, eq.(4.35), is rather tiny. Thus it is not surprising that in these conditions the branching ratio 
of the charged Higgs mode can be remarkably high: BR{t — > H^ b) ~ 50%, i.e. basically 50% — 50% 
versus the standard decay mode. In contrast, the branching ratio without SUSY effects (i.e. essentially 
the QCD-corrected branching ratio) is much smaller: at the characteristic 50(10) value, tan/3 = 50, it 
barely reaches 20%. Clearly, if the SUSY quantum effects are there, they could hardly be missed! 

As noted before, even though the dominant MSSM effects are, by far, the QCD and SUSY-QCD ones, 
they have opposite signs. Therefore, there is a crossover point of the two strongly interacting dynamics, 
where the conventional QCD loops are fully counterbalanced by the SUSY-QCD loops. This leads to a 
funny situation, namely, that at the vicinity of that point the total MSSM correction is given by just 
the subleading, albeit non-negligible, electroweak supersymmetric contribution: (5mssm — ^susy-ew- 
The crossover point occurs at tan/3 > 32 ~ mt/mi,^ where ^susy-ew ^ 20. For larger and larger tan/3 
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beyond mt/rri},, the total (and positive) MSSM correction grows very fast, as we have said, since the 
SUSY-QCD loops largely over-compensate the standard QCD corrections. As a result, the net effect on 
the partial width appears to be opposite in sign to what might naively be "expected" (i.e. the QCD sign). 
Of course, this is not a general result since it depends on the actual values of the MSSM parameters. 
In the following we wish to explore the various parameter dependences and in particular we want to 
assess whether a favourable situation as the one just described is likely to happen in an ample portion 
of the MSSM parameter space. In particular, the value tan/3 = mt/rrib = 35 will be chosen in all our 
plots where that parameter must be fixed. We consider it as representative of the low end of the high 



tan/? segment, eq.(4.8). Thus tan/3 = mt/mi, = 35 behaves as a sort of threshold point beyond which 
the MSSM quantum effects on t ^ H^ b take off so fast that they should have indelible experimental 
consequences on top quark physics. 



As regards to the non-supersymmetric electroweak corrections, few, it is apparent from Fig. 4.6(a) 
that they are very small, especially in the high tan/3 segment. Also in the very low tan/3 segment, 
0.5 < tan/3 < 1, Sew is relatively small; and this is so not only because both feiggs and 6sm become 
never too large in absolute value, but also because in that region there is a cancellation between feiggs < 
and ^SM > 0. As it happens, we end up with the fact that the complicated Higgs effects result in a very 
tiny contribution, except in the very low tan/3 end, where e.g. they can reach —15% at tan/3 ~ 0.5. In 
this corner of the parameter space, feiggs becomes the dominant part of ^mssm, being even larger than 
the QCD effects, which stay at the level of -8%, and also larger than the SUSY-QCD and SUSY-EW 
corrections, which remain below -1-4% and —1%, respectively. 

We have treated in detail the very low tan/3 segment by including the one-loop renormalization of 



the Higgs masses [112-11(:]. This is necessary in order to avoid that the lightest CP-even Higgs mass 
either vanishes at tan/3 = 1 or becomes lighter than the phenomenological bounds near that value. In 
passing, we have checked that the one-loop shift of the masses, as well as of the CP-even mixing angle. 



a, has little impact on the partial width of i — > H~^ b in the entire range of tan/3, eq.(4.6). They entail 
at most an additional 5% negative shift of iJmssm in the very low tan/3 region^. It is precisely in this 
region where the Higgs effects could be expected of some relevance, and thus where the renormalization of 
the CP-even mixing angle could have introduced some noticeable change in the neutral Higgs couplings. 
Quite on the contrary, at high tan /3 the corresponding effect is found to be of order one per mil and is 



thus negligibly small. On the other hand, a simple inspection of Figs. 4.5 and [4.6| (b) shows that even 
in the very low tan/3 ballpark, where there may be some ten percent effect from the Higgs sector, the 
rising of the tree-level width is so fast that it becomes very hard to isolate these corrections. We conclude 
that, despite the rather large number of diagrams involved, the over-all yield from the Higgs sector of the 



MSSM on i ^ H~^ b is rather meagre in the whole tan/3 range (4.6). This fact is somewhat surprising and 



was not obvious a priori, due to the presence of enhanced Yukawa couplings (1.2) in the whole plethora 



of Higgs diagrams. The cancellations involved are reminiscent of the scanty SUSY Higgs effects obtained 



for the standard top quark decay t — > W'^ b |161, 162| 



We come now to briefly discuss the standard QCD effects up to 0{as), which involve one-loop gluon 



corrections and gluon bremsstrahlung |^, ^. As it is plain from Fig. 16(a), Sqcb is negative-definite 
and very important in the high tan/3 segment. It quickly saturates for tan/3 > 10 at a large value of 
order —60%. Therefore, the QCD effects need to be considered in order to isolate the virtual SUSY 
signature |^,Q. The leading behaviour of the standard QCD component in the relative correction 
( 4.32| ) can be easily assessed by considering the following asymptotic formula 



2 ^^ 14^ {ml tan2 /3 -t- m^ cot^ /3) -f 3(4 -I- tan^ /3 - 2^ cot^ f3)ml In (^) 

Sqcd ^ -^ 27 — 2~5~; 2 — I2I3 ' ^^ ' (4.38) 

37r m^ tan-^ (3 + mf cot"' (3 

which is obtained by expanding the exact one-loop formula up to 0{m1/mf,M'^^/m'^). Here as = 
asim^), normalized as asiM"^) ~ 0.12. The big log factor ln(m(/TO^) originates from the running b- 
quark mass evaluated at the top quark scale. The correction is seen to be always negative. We point 
out that while we have used the exact 0{as) formula for the numerical evaluation, the approximate 
expression given above is sufficiently accurate to convey the general features to be expected both at low 
and at high tan/3. In particular, for nib 7^ and tan/3 in the relevant high segment ([4.^), the QCD 



^To perform that check, we have included both the stop and sbottoin contributions to the one-loop Higgs mass relations. 
A set of 7 independent parameters has been used to full y charac terize these effects, viz. {Mjj±,fj,,ts,nf3,m^ , m^ ,Af,,At). 



We refrain from writing out the cumbersome formulae [112-116 
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Figure 4.7: (a) The relative corrections S, eq.( 4.32| ), as a function of the gluino mass 
SUSY-EW, standard EW, SUSY-QCD, standard QCD contributions, and total MSSM contribution, 
(b) As in (a), but for the tree- level and full MSSM-corrected branching ratios ( 4.34 ). We have set 
tan/3 — mt/rrii, — 35. This value is maintained in all figures where tan/3 is fixed. 



correction becomes very large and saturates at the value 



2 a, 



'JQCD 



ilT^ - 15 

36 



In 4 



-62% (tan/3 >> ^/^mjm, ^ 6) . 



(4.39) 



(The exact 0{as) formula gives slightly below —60%.) At low values of tan /3, the corrections are much 
smaller, as it follows from the approximate expression Sqqd — (— cks/7'')(87r^ — 15)/18 ~ —12%. We 

, so that one would 
The limit mi, = 0, 



remark that for mi, = the dependence on tan/3 totally disappears from eq.(4.3 
never be able to suspect the large contribution (4.39) in the high tan/3 regime 
nevertheless, has been considered for the standard QCD corrections in some places of the literature but, 
as we have seen, it is untenable unless one concentrates on values of tan /3 of order 1 , in which case the 
relevance of our decay for SUSY is doomed to oblivion. This situation is similar to the one mentioned 
above concerning the SUSY-QCD corrections in the limit mi, = 0, which leads to an scenario totally blind 
to the outstanding supersymmetric quantum effects obtained for m,t, 7^ at high tan/3 |4^,Q. We stress 
that in spite of the respectable size of the standard QCD effects, they become fast stuck at the saturation 
value ( 4.3E ), which is independent of tan/3. On the contrary, the SUSY-QCD effects grow endlessly with 
tan/3 and thus rapidly overtake the standard QCD prediction. 



Worth notici ng i s the evolution of the quantities ( 4.32 ) and ( 4.34 ) as a function of the gluino mass 
(Cf. Figs. 4.7(a)-4.7(b)). Of course, only the SUSY-QCD component is sensitive to nig. Although the 
SUSY-QCD effects have been object of a particular study in Ref. [^,Q, we find it convenient, to ease 
comparison, to display the corresponding results in the very same conditions in which the electroweak 
supersymmetric corrections are presented. The steep falls in Fig. 4.7(a) are associated to the presence of 
threshold effects occurring at points satisfying rrig + mi_^ ~ m,t ■ An analogous situation was observed in 
Ref. 1 40 -42 1 for the SUSY corrections to the standard top-quark decay. Away from the threshold points, 
the behaviour of Ssusy-qcd is smooth and perfectly consistent with perturbation theory. In Fig. [4.7| (b), 
where the branching ratio (4.34) is plotted, the steep falls at the threshold points are no longer present 
since they are compensated for by the simultaneous opening of the two-body supersymmetric mode 
t ^ tig, for TTig < m.t — m^^ . 

We emphasize that the relevant gluino mass region for the decay t —>■ H^ b is not the light gluino region, 



but the heavy one, the reason being that the important self-energy correction mentioned above, eq.(4.29), 
involves a gluino mass insertion. As a consequence, virtually for any set of MSSM parameters, there is a 
well sustained SUSY correction for any gluino mass above a certain value, in our case rrig > 250—300 GeV. 
The correction raises with the gluino mass up to a long flat maximum before bending -very gently - into 
the decoupling regime (far beyond 1 TeV) El, Esf . The fact that the decoupling rate of the gluinos 
appears to be so slow has an obvious phenomenological interest. 
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Figure 4.8: (a) The relative corrections S, eq.( 4.32| ), as a function of the supersymmetric Higgs mixing 
parameter, /i, for the various contributions as in Fig j4.6| (a); (b) As in (a), but for the same branching 
ratios as in Fig. [4.7Kb). 



Next we consider in detail the sensitivity of our decay on the higgsino-gaugino parameters (/i, M) 
characterizing the chargino-neutralino mass matrices (Cf. section 2.3). We start with the supersymmetric 



Higgs mixing mass, fj,. As already stated above, we will largely concentrate on the ^ < case. Together 



with At > this yields At /i < 0, which is a sufficient condition [148, 164 170| (see also section 4.5) for 
the MSSM prediction on BRih — > s 7) to be compatible with experiment in the presence of a relatively 
light charged Higgs boson (as the one participating in the top decay under study). In fact, it is known 
that charged Higgs bosons of 0(100) GeV interfere constructively with the SM amplitude and would 
render a final value of BR[b — > s 7) exceedingly high. Fortunately, this situation can be remedied in the 
MSSM since the alternative contribution from charginos and stops tends to cancel the Higgs contribution 
provided that Aj /x < 0. Furthermore, one must also require relatively lig ht values for the masses of the 
lightest representatives of these sparticles, as well as high values of tan/? [14J, 164-17C]; hence one is led 
to a set of conditions which fit in with nicely to build up a favourable scenario for the decay t — > H'^ b. 

The evolution of the in divi dual contributions ( 4.37 ), together with the total MSSM yield, as a function 
of /i < 0, is shown in Figs. 18(a)-^^(b) for given values of the other parameters. We immediately gather 
from these figures that the SUSY-QCD correction is extremely sensitive to fj,. In fact, <^susy-qcd grow s 
rather fast with |/i|. This is already patent at the level of the leading Smh/rrn, effect given by eq.(09). 
In all figures where a definite ^ < is to be chosen, we have taken the moderate value /i = —150 GeV. 



Concerning the electroweak contribution, we noted above that the component Srrib/mb, eq.(4.31) 
actually decreases with fi. However, the 11 dep endences in the full (5sus y-ew are more complicated than 
in (^SUSY-QCD and cannot be read oif eq.( 4.31 ). This is evident from Fig. 4 
is fairly insensitive to fi 
QCD component 



a) where the total (Ssusy- 



EW 

(^MSSM, therefore, inherits its marked /i-dependence basically from the SUSY- 
As for the sensitivity of the corrections on the S'J7(2)i-gaugino soft SUSY-breaking 
parameter, M, Fig. 4^ conveys immediately that it is virtually non-existent. 

There is some slight evolution of the corrections with Ai, (Fig. 4.10(a)), mainly on the SUSY-QCD 
side. We realize that (5susy-qcd is not perfectly symmetric with respect to the sign of Ab- Once the sign 
/i < is chosen, the correction is larger for negative values of Ab than for positive values. We have erred 
on the conservative side by choosing Ab = -1-300 GeV wherever this parameter is fixed. As far as At is 
concerned, <^susy-qcd can only evolve as a func tion of that parameter through vertex corrections, which 
are proportional to At cot/3 (Cf. section 2.4.3 ); however, at large tan/3 these are very depressed. The 
electroweak correction i^susy-eWj instea d, is very much dependent on At; indeed, a typical component 
exhibiting this behaviour is given by e q.( 4.31 ), which is linear in At- The full dependence, however, is 
not linear and is recorded in Fi g.|4.10| (b). We realize that (^susy-ew and (5mssm chan ge s ign with At- 



The shaded vertical band in Fig. 4.10 (b) is excluded by our choice of parameters in Fig. 4.5 

Another very crucial parameter to be investigated is the value of mi 
QCD correction hinges a great deal on the value of the sbottom masses, as it is plain from eq.( 4.29 ) 



This is because the SUSY- 
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Figure 4.9: Dependence of 5, eq.( 4.32| ), on the S'C/(2)i-gaugino soft SUSY-breaking mass, M, assuming 

, is related to A 
and total contributions as in Fig.[4T6|(a) are shown. 



that the U{1)y gaugino mass, M', is related to M through M'/M — | ta.n^9w- The same individual 
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Figure 4.10: (a) The relative corrections 6, eg. ( 4. 32 ), as a function of the trilinear soft SUSY-breaking 
parameter Af, in the bottom sector. The other trilinear couplings are kept as in Fig. 8; (b) As in (a), but 
for the trilinear soft SUSY-breaking parameter At in the top sector. Shown are the same individual and 
total contributions as in Fig. 4.6 (a). 
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Figure 4.11: (a) The relative correctio ns 5 , eq.( 4.32| ), as a function of the lightest sbottom mass, 
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for the various contributions as in Fig. 4.6 (a); (b) As in (a), but for the same branching ratios as in 
Fig.0(b). 
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Figure 4.12: The relative corrections S, eq.(4.32), as a function of the lightest stop 
various contributions as in Figj4.6|(a). 



mass, mz , for the 



As a matter of fact, a too large a value of m^ may upside down the leadership of the SUSY-QCD 
effects. As a typical mass value for all squarks other than the stop we use rriq > 150 — 200 GeV(g 7^ i). 
From Fig. 4.11(a) we see that provided m^ < 300 GeV the SUSY-QCD effects remain dominant, but 



they steadily go down the larger is m^^ . The electroweak correction Ssusy-ew, on the other hand, is 
quite sustained with increasing rrif^ and there are parameter configurations where for sufhciently heavy 
sbottoms the supersymmetric electroweak effects are larger than the SUSY-QCD effects. However, this 



is not the most likely situation. The behaviour of the branching ratio is plotted in Fig. 4.11(b). 

Obviously, the evolution of the SUSY-QCD corrections with the stop masses is basically flat (Fig.ETa) 
since the leading contribution is independent of m^ . Therefore, for definiteness we fix TOj ~ 100 GeV. 
Nonetheless, if we wish to keep (5mssm > 0, we cannot go too far with mj:_^, for the electroweak correction 
is also seen to decrease with m^^. Indeed, whereas for mj^ = 65 — 100 GeV one has (5susy-ew ^ 20%, 
for m^^ > 250 GeV one finds (5susy-ew ^ 10%. For heavier stop masses, (5mssm becomes zero or slightly 
negative. In this situation, the imprint of SUSY lies in the fact that the total quantum effect is not as 



negative as predicted by standard QCD, eq.(4.39). 

T he in flu ence from the sleptons and the other squarks is practically irrelevant as it is borne out by 
Figs. 4.13 (a)- 4.13 (b). They enter the correction through oblique (universal) quantum corrections. The 
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Figure 4.13: (a) The relative corrections S, eq.(4.32), as a function of the up-squark masses m.u 



771^2 • The c-squarks are assumed to be degenerate with the up-squarks; (b) (5 as a function of the 
sneutrino masses, assumed to be degenerate. 



only exception are the r-sleptons Tq ("staus"), since they are involved in the process-dependent (non- 
oblique) contribution eq.(3.53), where the r-lepton Yukawa coupling becomes enhanced at large tan/3. 
For this reason, ^susy-ew in Fig- 4.13| (b) is somewhat larger the smaller is the r-sneutrino mass (assumed 
to be degenerate with the other sneutrinos). In all our calculation we have fixed the common sneutrino 
mass at 200 GeV. 

We have also tested the variation of our results as a function of the external particle masses, namely 
the masses of the top quark, bottom quark and charged Higgs. As for the external fermion masses, the 
corrections themselves are not very sensitive (see Figs. 4.14(a) and 4.15| (a)). Among the SUSY corrections, 
the most sensitive one on rut (resp ective ly on m i,) is (5susy-ew (resp. (5susy-qcd)- The branching ratios 
also show some dependence (Figs. 4.14(b) and 4.15(b)), especially on mi,. This effect is mainly due to 
the variation of the tree- level partial widths as a function of mt and mi,. As for the charged Higgs mass, 
Mh+, up to now it has been fixed at Mjj+ = 120 GeV, which is the preferred value for this mass at 
large tan/3 [p2|. We confirm from Fig. 4.16 (a) that there is nothing special in the chosen value for that 
parameter since the sensitivity of the correction is generally low, except near the uninteresting boundary 
of the phase space where the branching ratio (Fig. 4.1(;(b)) boils down to zero. 

We close our study of the corrections in the /i < case by plotting 5^ as a function of tan/3 (see 
Fig. 4.17 ). By definition, St is that part of (5mssm originati ng fro m the full process-dependent term A^-, 
eq.( 3.53 ), which stems from our definition of tan/3 on eg. ( 3. 50 ). This piece of information is relevant 
enough. In fact, it should be recalled that the quantum corrections described in the previous figures 
are scheme dependent. In particular, they rely on our definition of tan/3 given on eq.( 3.5Q| ). What is 
710^ scheme dependent, of course, is the predicted value of the width and branching ratio (Figs. |4.5| and 
[4.6| (b)) after including all the radiative corrections. Now, from Fig. 4.17 it is clear that the A7--term is 
not negligible, and so there is a process-dependence in our definition of tan/3, as it was announced in 
chapter |3[ At first sight, the (57--effects are not dramatic since they are small as compared to ^susy-qcd, 
but since the latter is canceled out by standard QCD we end up with 5^ being of the order (roughly half 
the size) of the electroweak correction (5susy-ew- 

The main source of process-dependent (5,--effects lies in the corrections generated by the r-mass coun- 
terterm, SmT/m,T, and can be easily picked out in the electroweak-eigenstate basis (see Fig. 4.18) much 
in the same way as we did for the 6-mass counterterm. There are, however, some differences, as can 
be appraised by comparing the diagrams in Figs. 4.4 and 4.18| , where we see that in the latter case the 
effect derives from diagrams involving r-sleptons with gauginos or mixed gaugino-higgsinos. An explicit 
computation of the diagrams (a) + (b) in Fig. 4.18 yields 
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Figure 4.14: (a) The relative corrections (5, eq.(4.32), as a function of the top quark mass within about 
2 (7 of the present experimental range at the Tevatron. (b) As in (a), but for the same branching ratios 
as inFigJO(b). 
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Figure 4.15: (a) The relative correction (5, eq.( tl.32|) , as a function of the bottom quark mass, (b) As in 
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Figure 4.16: (a) The relative corrections S, eq.( 4.32| ), as a function of the charged Higgs mass; (b) As in 
(a), but for the same branching ratios as in Fig. 4. 7(b). 
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Figure 4.18: Leading supersymmetric electroweak contributions to 6771^/771^- in the electroweak-eigenstate 
basis. 
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Figure 4.19: (a) The relative corrections S, eg. ( 4. 32 ), as a function of the hghtest sbottom quark mass, 
for positive ^ = +150 GeV, negative At and given values of the other parameters. In this case, huge 
values of toj; are needed in order to damp the absolu te v alue of the total correction down to 100%; (b) 
As in (a), but for the same branching ratios as in Fig. 4. 7(b). 



+ 



IGtt' 



■^M tan/3/(/i, m£;^,M), 



(4.40) 



where g' = gsyy/cw and M',M (Cf. section 2.4.4 ) are the Soft-SUSY-Breaking Majorana masses 
associated to the bino B and winos W^ ^ respectively, and the function /(t?t.i, 7712,77^ 3) is again given 
by eq.(4.3C). In the formula above we have projected, from the bino diagram in Fig. l.lS(a), only the 
leading piece which is proportional to tan/?. Even so, the contribution from the wino-higgsino diagram 
in Fig. 4.18 (b) is much larger. Numerical evaluation of the sum of the two contributi ons o n eq.( |4.4(]| ) 
indeed shows that it reproduces to within few percent the full numerical result (Cf. Fig. 4.17) previously 
obtained in the mass-eigenstate basis, thus confirmi ng t hat eq.( 4.40 ) gives the leading contribution. In 
practice, for a typical choice of parameters as in Fig. 4.5, this contribution is approximately canceled out 
by part of the electroweak supersymmetric corrections associated to the original process t — > H^ b, and 
one is effectively left with eq.(4.31) as being the main source of electroweak supersymmetric quantum 
effects at high tan/3. 

Finally, the corrections corresponding to the case /i > are studied in Figs. 4.19(a) and 4.19(b). The 
problem with /i > is that, then, the large SUSY-QCD corrections have the same (negative) sign as the 
conventional QCD effects, and as a consequence the total MSSM correction can easily blow up above 100%, 
the branching ratio becoming negative! To avoid this disaster (from the point of view of perturbation 
theory), we enforce the SUSY-QCD correction to be smaller than in the fi < case by assuming an 
"obese SUSY scenario" characterized by very large values for the sbottom mass (to^ — 600 GeV) and 
the gluino mass (rrig = 1000 GeV). We also choose At < so that the electroweak SUSY correction 
becomes opposite in sign to the SUSY-QCD correction (a feature that also applies in the /i < case, see 



58 



Quantum effects on t —f H^ b in the MSSM 



z,.u 


nij =600 GeV 


1.3 
1.0 


\ 






/ 


/ 


r„(t->H*b) 

r,cD('-> H* b) 

---r„,,„(t->H*b) 
r„(t->w*b) 




0.5 
n n 


^^ 


y 



10 20 30 40 50 60 

tanCP) 



Figure 4.20: The total partial width, Y(t -^ H^ b), including all MSSM effects, versus tan/3, for the same 
inputs as in Fig. 4. 19 (a), as compared to the t ree-level width and the QCD-corrected width. A typical 
value of mi within the range used in Fig. 4.19| (a) is selected. Also plotted is the tree-level partial width 
of the standard top quark decay t — > W^ b. 



Fig. [4.1G| (b)) and in this way the total SUSY correction is further lessened in absolute value. Incidentally, 
we remark that the simultaneous sign change of both /x and At is also necessary in order to keep At ^j. < 0; 
as noted above, this is required in order that the MSSM can be compatible with BR{b -^ 57) in the 
presence of a relatively light charged Higgs. In Fig. 4.20 we bring forward the effect of the new situation 



on the total partial width. In the present instance, the physical signature would be to measure a partial 
width significantly smaller than the one predicted by QCD. Clearly, the /i > {At /i < 0) scenario is not 
as appealing as the n < {At ^ < 0) one. 



4.5 Implications for the Tevatron data 



The results presented in previous sections indicate that the SUSY corrections to the decay under study 
could have a great impact on the search for charged Higgs bosons at the Tevatron. The CDF Col- 
laboration at the Tevatron has carried out direct searches for charged Higgs production in pp colli- 
sions at -^s = 1.8 TeV |134, 171, 172| ]. In these studies one is concerned with the final configurations 
tt -^ H+ H- bb,W+ H- bb, H+ W- bb. The latter would differ from that of the standard model, 
tt -^ W^ W~ b b, by an excess of states with one (or two) r-lepton "jets" (i.e. usually tagged in the 
hadronic decay mode) and two b-quarks and large missing transverse energy associated to the decays 
H~^ — > r+ Vr and/or H^ —^t^Vt- 

As stated in section [l.2| , if Mh± is not too close to the phase space limit, then there are two regimes, 
namely a low and a high tan/3 regime, where the partial width of the unconventional top quark decay 
becomes sizeable as compared to the standard decay t — > W'^ b. Nevertheless we shall focus only on 
the high tan/3 regime as it is this case that is correlated with the Higgs maximum rate into the r-mode 
versus the hadronic mode (Cf. eq.(4/7)). Clearly, the identification of the charged Higgs decay of the 
top quark could be a matter of observing a departure from the universality prediction for all the lepton 
channels through the measurement of an excess of inclusive (hadronic) r-events. However, from the 
non-observation of any r-lepton surplus, one may determine an exclusion region in the (tan/3, Mjj)-plane 
[134, 171-173| for any (Type II) 2HDM |39). The region highlighted in this plane consists of a sharply 
edged area forbidding too high values of tan /3 in correlation with Mh± . In the relevant SUSY region 
Mh± > 100 GeV (see below) the most recent analysis would imply that values in the range tan/3 > 40 
would be excluded |l7|,[l7|. 

In spite of its foreseeable importance, the impact of the SUSY quantum corrections on the dynamics 
oit^ ff+ b was not included in any of the aforementioned analysis [134, 171-173|. And this is especially 
significant in a decay like t — > iJ+ b whose sole exis tence could, in a sense, already be an indirect sign of 
SUSY. For, as is well-known, the CLEO data |174| on the radiative decays i?° ^ X^ 7 (viz. 6^57) set 
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Figure 4.21: The allowed region (shaded area) in the (Af , tan /3)-plane by the & — > 57 decay within the 
framework of the MSSM, and for a given set of inputs. 



a rather stringent lower bound on the mass of any charged Higgs scalar belonging to a generic 2HDM, to 
wit: Mh± > 240 GeV. Therefore, with only the W^ and H^ electroweak corrections, the charged Higgs 
mass is forced to lie in a range where the decay t -^ H^ b becomes kinematically blocked up. Of course, 
this is so because the virtual Higgs effects go in the same direction as the SM contribution. Fortunately, 
this situation can be remedied in t he MSSM wh ere the complete formula for the 6 — > s 7 branching ratio 
reads (see the extensive literature |148, 164- 170 1 for details): 



BR{b-y s-f) ~BR{b^ceV) 



I{mc/mb) [1 - ^as{mb)fQCDimc/mb)] 



(4.41) 



with 



^7 — ^SM 



A„- 



Ayc-q 



(4.42) 

being the sum of the SM, charged Higgs and chargino-squark amplitudes, respectively. (The contributions 
from the neutralino and gluino amplitudes are in this case generally smaller as they enter through FCNC.) 
Now, the important feature here is that the unwanted charged Higgs effects could to a large extent be 
compensated for by the chargino-stop contributions. And in this case a relatively light charged Higgs 
particle would perfectly be allowed in the MSSM for the decay t -^ H^ b to occur. 

In our renormalization framework, we use H'^ — > t~^Vt to define the parameter tan/3 through 



eqs. (3.5C) and (3.51), which allows to renormalize the tb H^-yertex in perhaps the most convenient 
way to deal with our physical process t -^ i/+ b ^ t^ i^r b. Apart from the full set of electroweak and 
strong corrections from the roster of SUSY particles (squarks, gluinos, chargino-neutralinos and higgses), 
we of course include the standard QC D correct ion Q. 



The results are presented in Figs. 4.21-4.24 



We point out that in the present work we have locked 
together the MSSM parameter space regions for the two decays b ^ sj and t —^ H^ b in order to find 
compatible solutions. In doing so we have used the full structure involved in eqs. (4.41), ( 4.42| ). Notice 
that recently the NLO QCD effects in the SM amplitude have been computed and the total error has 
diminished from roughly 30% to 15% (including the error in mt/mc) |175| ]. Also NLO order corrections 
to 6 ^ 57 in extensions of the SM have become available |17q-178|, but, on the other hand, new data 
from the CLEO collaboration makes the upper limit on this decay grow up [179|. The inclusion of these 
new information would be complicated and would not change our results. 



In Fig. 4.21 we determine the permitted region in the (tan/3, ^t)-plane in accordance with the CLEO 
data on radiative i?" decays at 2 a. For fixed /^ < 0, we find that At 11 < in the allowed region. This 
piece of information could be important since, as it is patent in that figure, the trilinear coupling At - 
entering the SUSY electroweak corrections - becomes strongly correlated with tan/3. This correlation 
depends slightly upon the value of the charged Higgs mass, Mh± , and it is built-in for the rest of the 
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Eff. 


i.T 

Pj, & geom. 




jets, Ht & ^T 


Process 




'^id^ 'Hd 




WW 


.16 (.13) 


.93 X .9 


.64 (.54) 


WH{80) 


.19 


X .87 X .5 


.61 


WH{100) 


.21 


= .36 


.62 


WH{120) 


.22 




.64 


WH{UO) 


.22 




.65 



Table 4.1: The efficiency factors for the It channel in WW and WH decay of tt at the CDF |T73|. For 
the WW process, the corresponding efficiencies from the CDF simulation are shown in parenthesis. The 
middle column shows the triggering, isolation and identification efficiencies from the CDF simulation. 
The total efficiencies £{1,2} are obtained by multiplying the three columns. 



plots (Figs. 4.22 - 4.24 ). From the SM result mentioned above, we have made allowance for an uncertainty 
of order 30% stemming from the non-included NLO corrections within the MSSM. 

For definiteness, and to ea se co mparison with the non-supersymmetric results, we will normalize our 
analysis with respect to Ref. [173|. Here the (?,T)-channel, with I a light lepton, is used to search for 
an excess of r-events. This should suffice to illustrate the potential impact of the MSSM effects on this 
kind of physics. To be precise, we are interested in the ti cross-section for the {I, r)-channel, aw, i.e. for 
the final states caused by the decay sequences ii — > iJ+ &, W~ b and H~^ -^ t^ i^t, W~ -^ I vi (and vice 
versa). The relevant quantity can be easily derived from the measured value of the canonical cross-section 
att for the standard channel t ^> blvi^t ^r bqq' , after inserting appropriate branching fractions, namely 



CTZt 



81^^ 



4 T{t^ Hb) 
9 r{t^Wb) 



£2 



f^ti 



(4.43) 



where the first term in the bracket comes from the SM decay, and for the second term we assume (at high 
tan/3) 1 00% branching fraction of H~^ into r-lepton, as explained before. Finally, e^ are detector efficiency 
factors [173 1, which we quote in table 4T. Notice that the use of t he measu red value of atf; |l80|] , instead 



of the predicted value within the standard NLO QCD approach |181, 182], allows a model-independent 
treatment of the result. In thi s re spect, we note that there could be MSSM effects on the standard 
mechanisms for ti production |183| (viz. Drell-Yan qq annihilation and gluon-gluon fusion) as well as 
corrections in the subsequent top qua rk d ecays |4y,^. To be concrete, we use the following value for 
the top pair production cross-section |18C]: 



ffft = 7.5 ± 1.5 pb . 



The number of events found in the (Z,T)-channel up to an integrated luminosity of llOpb^ is 4 |74,184], 
with an expected background of ^^ 2 events and ^ 1 event expected in the SM. This implies an upper 
limit of 7.7 events at 95% C.L., that is 

air < 70 fb (95% C.L.) . 

Therefore, proceeding in this way the bulk of the MSSM pay-off stems from the t -^ 7J+ b contribution 
in eq.(4.43). Specifically, in Fig. 4.22 we determine, as a function of tan/3 and for a given Higgs mass and 
fixed set of SUSY parameters, the cross-section for the {I, r) final state. There we show the tree-level (ao), 
QCD-corrected {(tqcd) and fully MSSM-corrected {(Tmssm) results. Of course, ctmssm includes both 
the SUSY-QCD and standard QCD effects, plus the MSSM electroweak corrections. Note that the QCD 
curve is similar to the one in Ref. ]173| ] M, but as it is also patent the full MSSM curve is quite different 
from the QCD one: in fact, the two curves lie mostly on opposite sides with respect to the tree-level 



curve!. This is the same effect we have seen in Fig. 4.5 translated to the cross-section determination 



The horizontal line in Fig. 4.22 gives the cross-section for the number of events expected in the (/, r)- 
configuration at the 95% C.L. after correcting for the detector efficiencies. Hence the crossover points 
of the three curves with this line determine (at 95% C.L.) the maximum allowed value of tan/3 for the 

^There is, however, a small difference due to the fact that we use the top quark scale, instead of the Higgs mass scale, 
to compute the QCD corrections. 
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Figure 4.22: The cross-section for the (r, ^)-channel (in fb) for the tree-level (gp ), QC D-corrected {(Tqcd) 
and fully MSSM-corrected {(Tmssm) cases, for the same parameters as in Fig. 4.21 . The horizontal line 
gives the 95% C.L. cross-section for the observation of the (t, /) final state. 



given set of parameters. It is plain that the MSSM curve crosses that line much earlier than the QCD 
curve, so that the tan/3 bound is significantly tighter than in the non-supersymmetric case. Notice that 
for this particular set of parameters the MSSM and tree-level curves turn out to meet the horizontal line 
at about the same point, which means that the SUSY effects fully counterbalance the standard QCD 
correction. We remark that this feature may occur for negative values of the higgsino mixing parameter 
(in Fig. 4.22, /i — —90 GeV). The situation with ^ > 0, with its very different corrections, is different 
and it will be discussed below. 



In Fig. 4.23 we present our results in the (tan/3, A//f±)-plane, by iterating the procedure followed in 
Fig. 4.22 for /I < and for charged Higgs masses comprised in the relevant kinematical range 100 GeV < 
A/^± < mt- The lower bound from the LEP constraint Mao > 75 GeV and the SUSY Higgs mass 
relations implies Mh± > 110 GeV. We also show the three exclusion curves for the tree-level, QCD and 
MSSM corrected cross -sections. The excluded region in each case is the one below the curves. By simple 
inspection of Fig. 4.23 , it can hardly be overemphasized that the MSSM quantum effects can be dramatic. 
Thus e.g. while for Mh± =110 GeV the maximum allowed value of tan/3 is about 50 according to the 
QCD contour, it is only about 35 according to the MSSM contour. We have also checked that, after 
all, the modulation of the latter by the 6 — *■ s 7 constraint is not too significant even when including the 
30% uncertainty mentioned above. For, it turns out that although the branch ing r atio for the b ^ sj 
decay severely limits the set of possible values of At for each tan/3 (Cf. Fig. 4.21 ), the corresponding 
impact on i — > H^ b is really minor. This is due in part to the fact that the supersymmetric electroweak 
corrections are not the dominant component in f — > H^ b, and also in part to the observed stabilization 
of its contribution within the region of parameter space allowed by 6 — > s 7. 

The above picture may undergo a significant qualitative change when we move to the /i > scenario, as 



can be easily guessed from the very different corrections obtained in it. This can be appraised in Fig. 4.24 



where we plot the excluded region in the (tan/3, M^±)-plane again for the same cases as before. Although 
not shown, we have also determined the portion of the (tan/3, ^t)-plane permitted by 6 — > 57 for n > 
(implying that At < 0), and checked that also in this case the influence on our top quark analysis is not 
dramatic. The point with the /i > scenario is that the MSSM curve is, in contradistinction to the /i < 
case, the less restrictive one. As a matter of fact it is even less restrictive than the original CDF curve for 



the inclusive t channel! (Cf. Ref. |134]). The reason being that for /i > 0, the SUSY corrections have the 
same (negative) sign as the standard QCD corrections (Fig. 4.19) and, therefore, the cross-section for the 
T-lepton signal becomes extremely depleted. In Fig. 4.24 we have chosen a heavier SUSY spectrum than 



in the previous figures in order to keep the total correction within the limits of perturbation theory. We 
see that for squark masses of several hundred GeV and a gluino mass of 1 TeV the excluded area can be 
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Figure 4.23: The 95% C.L. exclusion plot in the (tan/3, M/f)-plane for /i < 0. Shown are the tree-level 
(dashed), QCD-corrected (dotted) and fully MSSM-corrected (continuous) contour lines. The excluded 
region in each case is the one lying below these curves. The set of parameters is as in Fig. 4.21 , with At 
within the allowed region. 
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Figure 4.24: As in Fig. [4.23| , but for a /^ > scenario characterized by a heavier SUSY spectrum. 
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enforced to withdraw to a corner of parameter space. However, in this corner one cannot be precise any 
more since further reduction would make also the Higgs sector non-perturbative (see below). Hence the 
conclusion emerging for the case /i > is quite remarkable, to wit: relatively light ( > 110 — 120 GeV) 
charged Higgs masses within the kinematical range of i ^ iJ+ b could be allowed for essentially any 
admissible value of tan/3 within perturbation theory (i.e. tan/3 < 60 — 70). In other words, within this 
scenario one could not disprove the existence of relatively light supersymmetric charged higgses by the 
current methods of r-lepton analysis at the Tevatron. 

It is also interesting to compare our results with the bounds obtained from semileptonic and semi- 



tauonic i?- meson decays. In Ref. [185 



186] the excluded region in the (tan/3, M^±)-plane is computed 



for a general 2HDM whereas in Ref. |187] the corresponding MSSM analysis is performed and it is also 



confronted with the (uncorrected) top quark decay exclusion region. However, in the presence of the 



corrected results, we may compare Fig. 4.22 of the present work with Fig. 3 of Ref. [187| (both for /i < 0). 
We realize that the supersymmetric results on the top quark decay greatly improve the bound from semi- 
tauonic B decays across the crucial region defined by 30 < tan /3 < 65 and Higgs masses ranging between 
100 — 150 GeV. Even though for tan/3 > 65 the semitauonic S-meson decays are more restrictive, it 
should be pointed out that this range is already ruled out on sound theoretical grounds, namely by the 
breakdown of perturbation theory; for instance, the top quark Yukawa coupling with the CP-odd Higgs 
boson j4° would become grrib ta,nf3/2 Mw > 1- On the other hand, the /i > region is not so favoured 



by B-meson decays, but it is still compatible with experimental data at the la level for tan/3 < 40 [187| 



4.6 Conclusions 



From the explicit numerical analysis (section 4^), we have confirmed our expectations that the SUSY- 
QCD contribution to T(t — > H~^ b) is generally dominant. This conclusion would not hold only in some 
(unlikely) cases, e.g. if the gluino is very light and/or the lightest bottom squark is "obese" as compared 
to lightest top squark, i.e. if the former is unusually much heavier than expected. Furthermore, by 
restricting ourselves to the case /i < {At /i < 0) we confirm that at large tan /3 and for typical values of 
the parameters the total (standard plus supersymmetric) QCD correction largely cancels out, leaving a 
remainder on the SUSY-QCD side (Figs. p~5| - [4.6[ ). In all circumstances the virtual Higgs effects remain 
comparatively very small. Around tan/3 = mj/mf, ~ 35, one is left with basically the electroweak 
supersymmetric correction, Ssusy-ew^ which can be sizeable enough to be pinned down by experiment. 
However, as stated above, there is in general a strong remainder, Ssusy-qcd + Sqcd > 0, which grows 
very fast with tan/3 and it has the same sign as Ssusy-ew- In this favourable scenario, the virtual SUSY 
effects could be spectacular. This is true not only because in the relevant window of parameter space the 
SUSY quantum corrections are by themselves rather large, but also because they push into the opposite 
direction than the "expected" standard QCD corrections. As a result, the relative deviation between the 
MSSM prediction and the QCD prediction effectively "doubles" the size of the observable effect, a fact 
which is definitely welcome from the experimental point of view. 

From all the previous discussion there is one fact standing out which can be hardly overemphasized: If 
the charged Higgs decay mode of the top quark, t — > H'^ b, does show up with a branching ratio of order 
10% or above (perhaps even as big as 50%), a fairly rich event statistics will be collected at the Tevatron 



and especially at the LHC e.g. by making use of the identification methods described in Section 4.2. 
If, in addition, it comes out that the dynamics underlying that decay is truly supersymmetric, then 
the valuable quantum signatures that our calculation has unveiled over an ample portion of the MSSM 
parameter space should eventually become manifest and, for sure, we could not miss them. 

At present all the collected event statistics basically relies on our experimental ability to recognize the 
top quark decays originating from standard patterns (angular distribution, energy spectrum, jet topology 
etc.) associated to the usual Drell-Yan production mechanism. Notwithstanding, we wish to point out 
that it should in principle be possible to clutch at the supersymmetric virtual corrections associated to the 
vertex t b H^ also through an accurate measurement of the various inclusive top quark and Higgs boson 



production cross sections in hadron colliders. As an example, in Fig. 4.25 we sketch a few alternative 



mechanisms which would generate top quark production patterns heavily hinging on the properties of 



the interaction i6iJ*-vertex ^,18S]. Thus, while this vertex could be responsible in part for the decay 



of the top quark once it is produced, it might as well be at the root of the production process itself at 



LHC energies, where it could take over from Drell-Yan production |131 



We observe that in some of these mechanisms a Higgs boson is produced in association, but in some 
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Figure 4.25: Typical diagrams for top quark and charged Higgs production in hadron colliders involving 
the relevant tbH^-veitex. 



others (fusion processes) the Higgs boson enters as a virtual particle. Now, however different these 
production processes might be, all of them are sensitive to the effective structure of the tbH^-veitex. 
The top quark-charged Higgs associated production can be used to search for a charged Higgs boson at 
hadron colliders, using the single top quark production as a signature, and searching for an excess of rvr 
lepton pairs in this processes, in a kind of analysis similar to that of section 4.5. Work is currently in 
progress to compute these effects, and the discovery limit of a charged Higgs particle at the Tevatron 
Run II and at the LHC. Similar mechanisms can of course be depicted involving the neutral Higgs bosons 
of the MSSM interacting with ii and bb via enhanced Yukawa couphngs |^, 188 1. While it goes beyond 
the scope of this Thesis to compute the SUSY corrections to the production processes themselves, we 
have at least faced the detailed analysis of a partial decay width which involves one of the relevant 
production vertices. In this way, a definite prediction is made on the properties of a physical observable 
and, moreover, this should suffice both to exhibit the relevance of the SUSY quantum effects and to 
demonstrate the necessity to incorporate these corrections in a future, truly comprehensive, analysis of 
the cross-sections, namely, an analysis where one would include the quantum effects on all the relevant 
production mechanisms within the framework of the MSSM. For this reason we think that in the future 



a precise measurement of the various (single and double) top quark production cross-sections [189, 19C] 
should be able to detect or to exclude the tbH^-vertex as well as the vertices qqA^{h°,H°) involving 
the neutral Higgs particles of the MSSM and the third generation quarks q = t,b. Notice the fact that 
finding a light v4° either at LEP or at the Tevatron and not finding a H'^ below rrit through t -^ H^ b 
would not exclude the MSSM at all, provided {At < 0,/i > 0). In this latter scenario the bbh and bh 
channels production cross-sections will suffer from the same large negative corrections than the t b H^ 
vertex itself (see below). Then the most plausible production process of A^ , H^ and hP should be the 
associated production with a gauge boson (e.g. W'^) JL91[ . In fact, the observation of this latter channel 
and non-observation of the bottom-Higgs coupling would point toward the MSSM nature of this Higgs 
boson, with a definite prediction on the sign of the /x parameter: u> 



The finite threshold effects to the bottom quark mass -eqs. ( 4.29 ) and ( 4.31 )- will be on top of any 
observable proportional to the bottom quark Yukawa coupling. At one-loop, and maintaining only the 
leading finite contributions, this Yukawa coupling reads -see eq. ( 2.1S| )- 



grub 



V2Mw cos/3 



5mi, 
nil, 



(4.44) 



where the factor Smb/mi, is the net sum of expr. ( 4.29 ) and ( 4.31 ). Thus we could define an "effective' 
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bottom Yukawa coupling (4.44) which would give us a first estimation of the supersymmctric radiative 



corrections to any process in the bottom-Higgs supersector, such as A^ -^ hh [105, 188| or the production 
processes quoted above. However a full computation will always be needed to be sure this estimation 
really gives the bulk of the correction (see chapter ^. 

Whereas, on the one hand, one expects that some top quark partial widths will be determined with 



an accuracy of 10% at the upgraded Tevatron and perhaps better than 10% at LHC |150|, on the other 



hand we believe that from the point of view of an inclusive model-independent measurement of the total 



top-quark width, Fj, the future e"*" e supercollider should be a better suited machine [192,193]. For 



in an inclusive measurement, all possible non-SM effects will appear on top of the corresponding SM 



effects already computed in the literature [151-158|. Moreover, as shown in Ref. [192|, one hopes to be 



able to measure the total top-quark width in e+ e^ supercolliders at an unmatched precision of ~ 4% 
on the basis of a detailed analysis of the threshold effects in the cross-section, in particular of the top 
momentum distribution and the resonance contributions to the forward-backward asymmetry in the tt 
threshold region. Under the assumption that Th — ^w, and that the SUSY effects on Ft are purely 
virtual effects, it follows that a large SUSY correction of, say 50%, to t -^ iJ+ b translates into a 20% 
correction to F4. This effect could not escape detection. Thus, the combined information from a future 
e"*" e~ supercollider and from present and medium term hadron machines can be extremely useful to pin 
down the nature of the observed effects. 

From the study of the quantum effects on the top quark decay channel into charged Higgs particles 
we arrive at the conclusion that the recently presented r-lepton analyses by the CDF Collaboration at 
Fermilab are in general model-dependent and could be significantly altered by potentially underlying new 
physics. In particular, since in the absence of new interactions the results from radiative i3-meson decays 
generally preclude the existence of charged Higgs bosons below the top quark mass, it is reasonable to 
link the existence of the decay t -^ H^ b to the viability of the leading candidate for physics beyond 
the SM, viz. the MSSM. In this framework we find that, depending on the sign of the higgsino mixing 
parameter, /i, the recent r-lepton exclusion plots in (tan/3, Mj:/±)-space presented by CDF could either be 
further strengthened or on the contrary be greatly weakened. This dual situation could only be decided 
from additional experimental information unambiguously favouring a given sign of fi in other physical 
processes. 

In mSUGRA models in the literature one usually claims ^t < 0. This would imply a positive value 
for the fi parameter in order to comply with the b —> sj constraint. Thus it seems that the prediction 
from these models would be that there is no actual limit to the charged Higgs boson mass, and that there 
is no hope of seeing the t —> H^ b at collider experiments even if tan (3 were large. However the explicit 



prediction of mSUGRA models is [194, 195 



At = (l-r)^o-2mi/2 , 

r being the ratio of the Yukawa coupling squared with respect its value at the RGE fixed point. Its value 
at large tan/3 is r ~ 3/4. We can see that for mi/2 — 100 — 200 GeV and Aq > 1 TeV these models 
already predict At > 0. Together with /i < we already see that a large branching ratio for t — > iJ+ b is 
not excluded at all. 

We remark that although for brevity sake we have presented our numerical analysis for a given 
choice of the MSSM parameters, we have checked that our conclusions hold basically unaltered in ample 
regions of parameter space involving typical sparticle masses of a few hundred GeV. While the details of 
the exclusion plot in (tan /3, M^± )-space may depend on the particular channel used to tag a potential 
excess of r-leptons, all of these plots (and of course also the one from the inclusive measurement) should 
undergo significant changes. Finally, it is clear that similar considerations apply to experiments of the 
same nature being planned for the future at the LHC and at the LC. Thereby a general conclusion seems 



to consolidate [196]: In contrast to gauge boson observables, the MSSM quantum effects on Higgs boson 
physics can be rather large and should not be ignored in future searches at the Tevatron, at the LHC 
and at the LC. 



Chapter 5 

FCNC top decays into Higgs bosons 
in the MSSM 

5.1 Introduction 

In this chapter we perform the computation of the Flavour Changing Neutral Current (FCNC) decay of 
the top quark into a charm quark and a neutral Higgs particle in the framework of the MSSM, t —f ch 
where h is any of the neutral Higgs particles of the MSSM. We compute the contributions from the SUSY 
electroweak, Higgs, and SUSY-QCD sectors, in a sparticle mass model motivated by model building and 
Renormalization Group Equations (RGE). However, we neither restrict ourselves to a spectrum of any 
SUSY-GUT model (such as SUGRA) -which would constrain the masses in a narrow range-, nor to a 
generic, phenomenological motivated, spectrum -which would have too many parameters to play with. 



There exist some computations of FCNC top quark decays, both in the SM and in the MSSM 1 197 -206]. 
The Standard Model branching ratio BR{t -^ cH) is ~ 10^^"^ for Higgs boson mass around 80 GeV, and 
it decreases with the Higgs mass |207|] . There has been some work concerning the decay channel into 



gauge bosons {t ^> cV , V = ^, Z,g), see for example Refs. |201-205] for some works on the subject. The 
conclusion of these works is that the branching ratio of this decay is at most 10^^, maybe a bit larger 
in the gluon channel. However, to our knowledge, there are not so many works on the FCNC decay of 



the top quark into Higgs in the MSSM |206], and they are not so complete as in the case of the gauge 
bosons. For example in [206| it is concluded that the branching ratio for the decay channel i —> c/i in 
the MSSM is at most of 10"^, for the SUSY electroweak contributions, and 10"^ for the SUSY-QCD 



contributions. However we think that the work of [206| is not complete. They do not include effects 
of the Higgs particles in the loops, and they do not take into account the q^ Qr h vertices, so they miss 
the potentially large contributions coming from the trilinear soft-SUSY-brcaking terms At^i, ( p.2C| ), and 
from the higgsino mass parameter fi. We find that a full treatment of the SUSY-QCD contributions may 
greatly enhance the FCNC width by some orders of magnitude. Therefore, a more general and rigorous 
computation of the decay i — > c /i is mandatory. 



5.4 we 



In section 5.2 we make a summary of the technics of the computation. In sections 5.2 and 
present our results for the SUSY electroweak and the SUSY-QCD contributions to the decay width t — > c /i 
respectively. Finally we present the conclusions. 

5.2 One-loop FCNC decays 

The computation of FCNC processes at one-loop, unlike the other calculations presented in this Thesis, 
does not involve renormalization of parameters or wave functions, so one is left only with the computation 
of the different diagrams that contribute to the process. The generic type one-loop Feynman diagrams 



contributing to the decay under study are in Fig. 5.1. The vertex diagram V follows after a straightforward 



calculation. As for the diagrams St and Sc we define a mixed self-energy, 

j:tcik) = ltJ:L{k^)PL+lt^R{k^)PR + mt {i:Ls{k'')PL + ^Rs{k')PR) (5.1) 

-where the rrit factor multiplying the scalar part is arbitrary, put there only to maintain the same units 
between the different S^. 
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Figure 5.1: Generic one-loop Feynnian Diagrams contributing to t — > ch. 



To present the expressions of this computation we shall introduce a notation that allows to treat the 
three possible decays in an unified way. To this end we introduce a vector of neutral Higgs fields 



as a function of which interaction Lagrangian with up-type quarks reads 



2MH/sin/3 



r=l,3 



the K matrix being 



if°" = 




Now we are ready to give a general expression of the effects of S^ to the amplitude t ^ c<J>51: 
—igrut 1 



in 



-iT. 



St 



2Mv^sin/3m2-m2""^^^ ^ 

Pl Kf [ml^R{ml) + m,mt {^rA^D + ^^^1)) + m^ ^Ls{ml)] 

+ PB.iK°y [L^R]}utik) 



2MH'sin/3TO^ - 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



After adding up the vertex contributions from diagram V (Fig. 5.1) to the expressions (5.5) we can 
define an "effective" vertex 



i T = ^tgu.ip) (Fl Pl + Fr P^) ut{k) 



(5.6) 



We have taken into account all three generations of quarks and squarks, and have performed the usual 
checks of the computation, in particular that the form factors Fl and Ffi are free of divergences before 
adding up the three quark generations, both analytically and numerically in the implementation of the 
code. 
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Figure 5.2: One-loop electroweak vertex diagrams contributing to the decay t ^ ch. d ((i{a,fc}) represent 
mass-eigenstate down type quarks (squarks) of any generation. 



After squaring the matrix element (g_^), and multiplying by the phase space factor, we can compute 
the decay width, 



r{t^ch) = 



and define the ratio 



32 7rm? 



\l/2r 2 2 2 



[(- 



ml){\FL\^ + \Fr\^) + 2mtm, (Fl F^ + Fl Fr)] 



B(t 



cri) = 



ch) 



T{t^hW+) 



(5.7) 



(5.8) 



which will be the main object under study. This ratio is not the total branching fraction of this decay 
mode, as there are many other channels that should be added up to the denominator of (5.8) in the 
MSSM, such as the two a nd th ree body decays of the top quark into SUSY particles, and also the decay 
channel t -^ H^ b |5(]|,|65, p09 |. For the mass spectrum used in the numerical analysis (see sections 5^ 
and |5.4|) the former are phase space closed, whereas the latter could have a large branching ratio. 



5.3 SUSY-EW contributions 



For the electroweak contributions to the decay channel t -^ chwe work in the so called Super-CKM basis, 
that is, we take the simplification that the squark mass matrix diagonalizes as the quark mass matrix, so 
that FCNC processes appear at one-loop through the charged sector (charged Higgs and charginos) with 
the same mixing matrix elements as in the Standard Model (the CKM matrix). 

We have taken into account the contributions from charginos (xt) s-nd down type squarks {da, 
a = 1, 2, . . . , 6 = di,d2, ■ ■ ■ ,b2, the mass eigenstates down squarks), and from charged Higgs and Gold- 
stone bosons {H^, G ^) and down type quarks {d, s, b). The various diagrams contributing to this decay 
can be seen in Figs. 5_^ and |]^. We have not included the diagrams with gauge bosons (W"'') as the 
largest contributions will come from the Yukawa couplings of the top and (at large tan/3) bottom quarks. 
However, the leading terms from longitudinal W^ are included through the inclusion of Goldstone bosons. 



5.3 SUSY-EW contributions 
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H- , G- 




Figure 5.3: One-loop electroweak diagrams contributing to mixed t — c self-energy, d {da) represent 
mass-eigenstate down type quarks (squarks) of any generation. 



5.3.1 Vertex and self-energy functions 



To write down the concrete form of the various contributions to the form factors (5.6) we generalize 
the compact notation introduced in section 5.2, thus we define a vector of charged Higgs and Goldstone 
particles 

$+ = (i/+,G+), (5.9) 

and write down the interaction Lagrangian of down-type squarks with Higgs particles analogously to ( ^.3| ) 



C 



*d 



+ 



2Mwcos(3 
9 



Y, <^''rd{K"/PL + {K'^'rPR)d 



r=l,3 






h.C. 



where the K matrices are 







J^+ud _ I cot P 



(5.10) 



(5.11) 



The interaction Lagrangian between up-type quarks, down-type squarks and charginos can be read di- 
rectly from the expressions (2.36), multiplying by the appropriate element of the CKM matrix. In a 
compact notation we write it as 



C 



udx^ 



-gVuddl^t {At't 



^PL + A^'al^Pn) u + h.c. 



with u (d) up-type quarks (down-type squarks) of any generation, the coupling matrices being 



A 



{d.u) 



j(d) 



id)* 



■Kla ^il - ^d ^2a K2 7 ^ 



.{d.u) 



'^la ^uUi2 



We similarly define a 3-dimensional matrix containing the triple Higgs interactions asn 



(5.12) 



and the chargino couplings to neutral Higgs 

^x+x+* = -5 E <xt(Wl^L Pl + Wl^j, Pr)xI 



r,i,] 



(5.13) 



These Brst and W/" matrices are the corresponding Feynman rules (divided my ~ig) of the respective 
processes and can be found in |3^ -they can also be found in ||71|] where there is also a detailed explanation 
of how to obtain themg. 

^Note that elements -Bijs are complex and Birj, = 0. 

^We have generated all the Feynma n rul es derived the scalar potential (Higgs particles self-couplings and squark-Higgs 
couplings) by means of a Mathematica [21C| code. 
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The contributions from the various diagrams in Fig. 5.2 can be written generically as, 



Na 



{Ci2 - C,i)mt 4^) Af - Ci2 m, 4^^ A^^'^ + Co tua 4^ 4'^ 



Fr - FLiA 



(*) 



4*)) 



(5.14) 



for diagrams |5.2| (a) and (d). Diagrams 5.2(b) and (c) contributions can be written as 



Fr 



Nr 



i(i) n(2) n(3), 



i(i) n(2) r)(3), 



Co {D'l' D'^' D'^> rric rut + D^'' D'£> D'£' m^ moi 



+D^j^'> i^f D^j^'> mt mD2 + o'^rt ^f ^f ^Di ^02) 
+C12 m, (i^W ^(2) ^(3) ^^ ^ ^W ^(2) ^(3) ^^ 

, n(i) r)(2) r)(3) ,-„ , n^^) n(2) n^^) ^ ^ 
HCii - C12) mt [D^l'^ Df^ D^i^ m, + d'^^^ of of mt 

+i?« i^i?) Df m,, + D^^ D^^ D^^ m^2) + Co D^^ i^g) i^f ) 



Fr = Fl{D 



(*) 



D%^) 



(5.15) 



Now we are ready to write down the different contributions, by giving the values of the different matrices 
and masses appearing in the expressions above, thus to obtain the vertex functions of the decay t ^r c^^ 
we must apply the following rules 



diagram 5.2(a): take (5.14) and assign 



A 



(1) 



A{d,c) .(1) 
^+bi ' ^R 



.{d,c) 



l(2) 



Ad-t) A2) 



A 



{d.t) 



ruA = M, , Na = i g^VtdVcdRi^J (^IV G 
C* =C^{k,-p' ,Mi,mj, ,tot) , 



/e ' 



where G^^ is the Feynman rule for $" — > dlj d'* divided by —ig, d\ 2 are the electroweak eigenstates 
of down type squarks 1 39 , 7l|p| . 



diagram 5.2(b): substitute in (5.15) 



I?« = A 



(d^c) 
+aj 



l(l) 



.(d,c 



(2) 






D 



(3) 



Aid,t) r)(3) 
^+ai ' ^R 



A 



-a. J 
(d,t) 



ijL I ^R. 



m_Di = Afj , mD2 = Mj , No ^ i g Vtd Vcd 
a ^C,{k,-p',m^ ,M,,M,), 



d1?^ = w,":- 



ijR 



diagram 5.2(c): the following assignations must be performed to (5.15) 



d'l' - K^l' , D^r' 



(2) 



D 



(3) 



^iR. ' ^R 



(3) 



^iR ' ^L 
^iR 



K. 



Od 

i 1 ^R 



D 



(2) 



t)d\* 



(Kn 



rriDi = ■mD2 = mi, , Nd = « 



■TOd 



AML cos/3 



T4dKd 



C* = C* (fc, -p' , m^+ , mrf, md) , 



• diagram 5.2(d): make the following substitutions to ( 5.14 ) 

^2 



li^^ 



JV^; _ T^+cd /.(I) _ 7^+cd 4(2) _ T^+td A^) _ T^+td 



^iR 



mA = md , Na ^ i 



9 



2Af2,-^ 



Bijk Vtd Vcd 



C* (fc, -p', ma, m^+ , m^+ ) 



We recall that our convention for the ji parameter is opposite in sign to that of 
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As can be noted from above expressions the form factors induced by Higgs mediated diagrams -Fig. 5.2(c) 
and (d)- have the property F^ = Fn for iJ° and /i°, and Fl = —Fji for A'^. This only form factor for 
each one of the different Higgs mediated diagrams can be written in a more convenient form, but then 
we should write down 16 different expressions! 

The one-loop mixing Feynman diagrams between the two mass-eigenstates quarks t and c can be seen 
in Fig. 5.3, their contribution to the mixing self-energies (pj) can be written as follows: 



Sfl(fc2)| 

rrit SLs(fc^)| 
Sii(fc') 



(a) 
(a) 
(a) 
(a) 



mt Sfls(fc^) 
mt Y,Ls{P) 



1(6) 



(6) 



(b) 



= ig^Vt,V,aAf:^ Af^} B,{k,M,,m^J 
= I g^ Nh Vtd Vcd Af^^ A^t*^ Bo(fc, M^, m^^ ) 
- ^ 5^ M, Vtd Vcd A^t? ^+;*? So(fc, M„ mj ) 



^g TOc Tnt 



2 Ml 



Vtd V,d [coi" 13 {Bo + Bi) (fc, Mh± , md) 



9 9 

^g m^d 
2Ml 



-{Bo + Bi){k,Mw,md)] 

Vtd Vcd [tan2/3 (So + Bi) (fc, Mh± , md) 



{Ba + Bi){k,Mw,md)] 



I g mt m^ 

« (?^ m,cm^ 
2M2, 



^td Kd [So(fc, Mjj± , md) - Bo(fc, Mh/, TOd)] 



Vtd Vcd [Bo{k, Mh± , md) - So(fc, Mvy, m^)] 



(5.16) 



The compact form of these self-energies allows to avoid the use of the cumbersome notation we used for 
the vertex factors. 



5.3.2 Numerical analysis 

With all these expressions we are now ready to look at the numerical results. We plug in all these 
contributions in {').(:) and (5. J;) and e valuate numerically the expression (5.8). The input parameters 
chosen to illustrate the results in Figs. 5.4-5.5 are: 



tan/3 = 35 , /i = -200 GeV , M = 150 GeV , M^o = 80 GeV 



150 GeV , mr 



At ^ An 



300 GeV , Ab 



= 200 GeV , 
-300 GeV 



(5.17) 



where mi^ , m^ are the lightest t and h mass, and all the masses are above present experimental bounds. 
This s ome what light value of the pseudoscalar Higgs mass is not essential in the results, as can bee seen 
in Fig. 5^ (d). We have chosen a SUSY mass spectrum around 200 GeV, which is not too light, so the 
results will not be artificially optimized. We have also checked all through the numerical analysis that 
other bounds on experimental parameters (such as 5p) are fu lfille d. 

In Fig. 5.4 we have plotted the different form factors of (5.6) as a function of tan/3 for the channel 
with the lightest scalar Higgs {hP). We can see that the contributions from the Higgs sector and the 
contributions from the chargino sector are of the same order. It turns out that they can be either of the 
same sign, or of opposite sign. The chosen negative value for A^ is to make the two contributions of the 
same sign. It is also clear that in both cases Fr ^ F^. This can be easily understood by looking at the 
interaction Lagrangians involving higgsino-sbottom-charm and Higgs-bottom-charm: 



^hbc 
C-Hbc 



-g Vcb c (i?fj \c Pl + 4a At Pr) X+foa + h.c. 



g 



V2M 



-Vcbc{mc cot(3PL + m,b ta.n(3 P^) bH^ + h.c. 



(5.18) 



w 



we can see that in both of them the contribution to the right-handed form factor will be enhanced by the 
Yukawa coupling of the bottom quark, compared with the charm Yukawa coupling that will contribute 
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t->ch 
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b 

< 0.0 
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Fl Higgs 
Fr Higgs 
F^ Charg 
Fr Charg 
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tan(P) 



50 



60 



Figure 5.4: Differe nt for m factors (5^) for the channel t 
set of inputs of eg. ( 5. 17 ). 



chP as a function of tan /?, with the typical 



to the left-handed form factor. On the other hand we have checked that the inclusion of the first two 
generations of quarks and squarks only has an effect of a few percent on the total result. 

In Fig. 5.5 we can see the evolution of the ratio ( |5.8| ) with various parameters of the MSSM, by taking 
into account only the electroweak contributions. The growing of the width with tan/3 (Fig. |5.5| (a)) shows 
that the bottom Yukawa coupling plays a central role in these contributions. The evolution with the 
trilinear coupling A^, and the higgsino mass parameter ^ -the two parameters that appear in the trilinear 
coupling 5l bRh- displayed in Figs. |5.5| (b) and (c) shows that these parameters can enhance the width 
some orders of magnitude. We have artificially let A}, grow up to large scales (that are not allowed if one 
wants that squarks do not develop vacuum expectation values) in order to emphasize the dependence on 
A}y. The various spikes in these figures refiect the points where the form factors change sign, whereas the 
shaded region in Fig. |5.5| (c) reflects the exclusion region of /i by present LEP bounds on the chargino 
mass. 

is smaller for the heaviest scalar Higgs [H'^) because with the 
-it/2, so making the couplings of H^ with 



In all these flgures the ratio (5. 
parameters ([5.17) the CP-even Higgs mixing angle a is near 



down quarks and squarks much weaker, but in fig. |5^ 

Higgs mass grows (and this shifts a far away from — tt/ 

We conclude that the typical value of the ratio ( |5J 

around 200 GeV, is 



(d) it can be seen that when the pseudoscalar 
2) the two scalar Higgs bosons change roles. 
1), at large tan/3 < 50 and for a SUSY spectrum 



B 



SUSY-EW 



{t ^ ch) ~ 0{lQ' 



-«) . (5.19) 

^. The maximum value is found to be in 



In favorable regions of the parameter space it can grow up to 10 

the ballpark of several 10^^, for sufficiently large Af,. This is an improvement of the previous result |2G6| 

specially in the A^ channel, by two orders of magnitude. 



5.4 SUSY-QCD contributions 

The gluino-mediated supersymmetric strong interactions in the MSSM can also produce FCNC processes. 
This occurs when the squark mass matrix does not diagonalize with the same matrix as the one for the 
quarks. We introduce then intergenerational mass terms for the squarks, but in order to prevent the 
number of parameters from being too large, we have allowed (symmetric) mixing mass terms only for the 



left-handed squarks. This simplification is often used in the MSSM, and is justified by RGE analysis |211| 
The mixing terms are introduced through the parameters 5ij defined as 



{M\ 



LLJij 



Dlj Tfli TTlj , 



(5.20) 



where rrii is the mass of the left-handed i squark, and to? is the mixing mass matrix element between 
the generations i and j. Thus we must diagonalize two 6x6 mass matrices in order to obtain the 
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Figure 5.5: Evolution of the ratio (WM) with (a) tan/?, (b) the trihnear couphng Ab, (c) the higgsino 
mass parameter fi, and (d) the pseudoscalar Higgs mass M^o. The rest of inputs are given in eq.(pTT7). 



mass-eigenstates squark field s. Fo llowing a notation similar to the standard one we introduce the mixing 
matrices as follows (Cf. Sec. 2.4.3 ) 



i?(«)^Al?i? 






^a!3<lp 



M 



qD 



diagjm? . . , m? } , q = u, d 



(5.21) 



where A4^^ r. is the 6x6 square mass matrix for up-type (or down-type) squarks in the EW basis, with 

indices a = 1, 2, 3, . . . , 6 = ul-, uu, cl, . . . ,tji for up-type squarks, and an equivalent choice for down-type 
squarks. In this study we are only interested in the up-type quarks-squarks system, so we will drop out 
the (q) super-index in the forthcoming expressions. The rotation matrix R introduces gluino mediated 
tree-level FCNC between quarks and squarks, the corresponding interaction Lagrangian can be deduced 
using the very same formalism of the "ordinary" SUSY-QCD (2.37) interactions, but using the more 
general rotation matrix (5.21), 



-SUSY-QCD 



yEi'c [R*5a Pl - Rla Pr] Qa.i Kj ^J 

-^ ^? [Rla Pl - Rlc Pb] q*a.i Kj C] 

^ re [Rio. Pl ~ RL Pb] £,. Kj ^, ■ 



(5.22) 
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Figure 5.6: One-loop SUSY-QCD vertex diagram contributing to the decay t 
mass-eigenstate up type squarks of any generation. 



ch. U{a,i3} represent 




Figure 5.7: One-loop SUSY-QCD diagrams contributing to mixed t — c self-energy. Ua represent mass- 
eigenstate up type squarks of any generation. 



5.4.1 Vertex and self-energy functions 



Using the Lagrangian (5.22) one can find the SUSY-QCD one-loop contributions to the process under 
study, which Feynman diagrams are depicted in Figs, b.t and 5.7. The vertex diagram contributions to 
the form factors (5_^) in Fig. |5.6| can be written as 



+mc Rz[3 Rla C'i2 + rug R^p RlaCo] 
Fr = Fi (3 <-> 4 , 5 ^ 6) 
N ^ i%T:asCFR*pGljRjo. 



(5.23) 



where G[- is the Feynman rule of <&" 
up- type squarks [p9|,|71[. 



^u'* divided by —ig, with u[ ^ the clcctroweak eigenstates of 



The one-loop mixing self-energy in Fig. 5.7 takes the following form 



SL(fc') 

mtT,Rs{k^) 



-i 2 77 as Cp i?3Q Rla Bi (-fc, rUg , ma^ ) 

-i 2 TT Cks Cf Ria Rla ^1 (-fc, TOg , m£i^ ) 

-i 2 TT cks C'f nrig R^a Rla Boi-k, rrig, m^^ ) 
-i 2 TT a^ Cf rrig R^^ Rla Bo{~k, nig, nii,^ ) 



(5.24) 



where Cf = {N^ — ^)/2Nc — 4/3 is a well known colour factor. Now we can introduce these expressions 
in (5.5) and (p^) to obtain the relevant branching ratio under study (|5.8|). 



5.4.2 Numerical Analysis 

For the numerical analysis we must provide as input parameters, apart from that of the Higgs sector, the 



various squark masses and mix ings , i.e. the 5 parameters of (5.20). These 5 parameters are constrained 
by low energy data on FCNC | 212| , 213 ]. The bounds have been computed using some approximations. 
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Figure 5.8: Evolution of the ratio (5.8) with (a) the mixing parameter between the 2nd and 3rd squark 
generations (523, (b) the higgsino mass parameter /i. (c) the ghiino mass m~g^ and (d) the pseudoscalar 
Higgs mass M^o, the rest of inputs are given in eqs.(p3.17) and (5.26). 



so they must be taken as order of magnitude hmits. We use the foUowing bounds |212, 21 



1^121 
I'5l3| 
1^231 



< .lV"ifim£/500 GeV 

< .098 ^mu mt /500 GeV 

< 8.2m£mj/(500GeV)2 , 



(5.25) 



For the various parameters that are common to the EW analysis we use the same input parameters ( |5.17| ) 
to which we must add the specific parameters of the SUSY-QCD sector, namely 



150 GeV 

0.03 0.03 
0.6 




(5.26) 



A comment is in order for the present set of inputs: we have introduced in (5.17) the lightest stop 
mass as an input, and this stop is mostly a ijj. However, in this new parametrization we introduce this 
mass as the lightest Ua mass, which will be mostly a t^. 



Again the largest contribution comes from the right-handed form factor of (5.6), but this is only 
because we have chosen not to introduce mixing between right-handed squarks. 



We have plotted the evolution of the ratio ( 5.8 ) with some parameters of the MSSM in Fig. 5.8 . As can 
be easily guessed, the most important parameter for these contributions is the mixing mass parameter 
between the 2nd and 3rd generation of left-handed squarks, the less restricted one of the three (eq. ( p3.25 ) ) . 
In Fig. 5.8 (a) it is shown that changing Si-^ by 3 orders of magnitude, the ratio (5.8) can increase by 7 
orders of magnitude! We can see in F ig.|5.8| (b) that the [i parameter also plays an importan t ro le, like 
in the electroweak contributions (Fig. |5.5| (c)), and for the same reasons, bringing the ratio (5^) up to 
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values of 10~^. Notice that the central region of |/i| < 90 GeV is excluded by present LEP bounds on 
the chargino mass. 

The evolution with the gluino mass (Fig.|5.8| (c)) is asymptotically quite stable, showing a slow de- 



coupling. Finally in Fig. 5.8 (d) we have plotted the evolution with the pseudoscalar Higgs mass, it is 



also quite stable, until near the kinematic limit for A^ and H'^ . 

We conclude that the typical value of the SUSY-QCD contributions to (5.S), with a SUSY spectrum 
around 200 GeV, is 

B^^^^-^c^{t^ch)~O{l0-^') , (5.27) 

but in favourable regions of the parameter space (i.e. large /i, or relatively light gluino) it can easily reach 
values of 10~*. The upper bound is at several 10"*. This is 1-2 orders of magnitude larger than the 
previous estimate |206|. 



5.5 Conclusions 

We have computed the SUSY-EW, Higgs, and SUSY-QCD contributions to the FCNC top quark decay 
t ^> ch {h = h^^H^^A^) in the MSSM, using a mass spectrum motivated, but not fully restricted, by 
model building and Renormalization Group Equations. 

We have found that with a SUSY mass spectrum around 200 GeV, which is well above present bounds, 
the different contributions to this decay are typically of the order 

5SUSY-QCD(^_^^^) ~ 10-5-10-4. (5.28) 

The difference of at least two orders of magnitude between the two contributions makes unnecessary 



to compute the interference between the two contributions, but if the limits on 523 (eq. (5.25)) improve, 
it should be necessary to make the full computation. 



The results (5.2?) are an improvement of the previous estimate [206|, specially in the A^ channel, 
thanks to the inclusion of the ql Qr h vertex. 

It would probably be difficult that this decay can be measured either at the Tevatron, or at the NLC, 
but there exists a possibility for LHC. As an example to assess the discovery reach of these accelerators 



the FCNC top quark decays into a vector boson are [214 



LHC: B{t^cV) > 5 x 10"^ , 

NLC: B{t~^cV) > 10"=^ - lO"" , (5.29) 

where the lack of sensitivity of the NLC is due to the lower luminositj^. So, if the discovery reach for 
FCNC Higgs processes are not very different from that of the gauge bosons, there is a possibility to 
measure this decay channel at the LHC even if SUSY particles are not seen at the LEPH. 



*Note added: This estimate 



Collider expect to reach 500 fb^^ |222 



s fo r a 50fb integrated luminosity |214|. Present studies for TF. ST.A future Linear 



then it would be possible to measure this ratio at the LC [E23 



Chapter 6 

One-loop corrections to scalar quark 
decays 



6.1 Introduction 

Sparticles not much heavier than a few hundred GeV could be produced in significant numbers already 
at the Tevatron. For instance, selectron production was advocated in Ref. ||25|| to explain a purported 
non-SM event in the CoUider Detector at Fermilab (CDF). Subsequently, in Refs. | 132| , 159| it was argued 
that half of the top quarks at the Tevatron might come from gluino decays into top and stop, g — > tti. 
Similarly, as discussed in the introduction, we may envision the possibility that sbottom squarks are pair 
produced by the usual Drell-Yan mechanism and then decay into top quark and charginos: ba -^ tx7- 
Indeed, this would be the leading two-body decay if gluinos are heavy enough that the strong decay mode 
ba ^ gb is kinematically blocked up^. The observed cross-section would then be the one of eq. (L4). We 
shall assume thorough present study that gluinos are much heavier than squarks, so that their contribution 
to this cross-section through qq ^ gg foll owed by g -^ tii is negligible. We could have non-SM top 
quark decay modes, such as e.g. t ^^ iaXa |f32| , 159[ and t -^ i/+ b |^,Q, that could serve, pictorially, 
as a "sinkhole" to compensate (at least in part) for the unseen source of extra top quarks produced at 
the Tevatron from sbottom pair production (Cf. eq.( |l.4| )). As stated in section 12 one cannot exclude 
that then non-SM branching ratio BR{t -^ "new"), could be as big as the SM one, i.e. ~ 50%. 

If tan/3 is large and there exists a relatively light chargino with a non-negligible higgsino component, 
the alternative mechanism suggested in eq.( |l.4| ) could be a rather efficient non-SM source of top quarks 
that could compensate for the depletion in the SM branching ratio. 

While the squark production cross-section has already received some attention in the literature at 
the level of NLO radiative corrections [218- ^20 |, an accurate treatment of the decay mechanisms is 
also very important to provide a solid basis for experimental analysis of the top quark production in 
the MSSM. Thus in this chapter we consider the computation of the QCD and leading supersymmetric 
electroweak (SUSY-EW) quantum effects on ba -^ txj , namely the ones induced by potentially large 



Yukawa-couplings from the top and bottom quarks (2.18) 



6.2 Vertex renormalization 

The one-loop Lagrangian of the btx~ interaction follows after substituting the one-loop countcrterms, 
that of &, t, and x~ from chapter O into the bare Lagrangian (2.35), 



with 



X tb 



SCI' 



-9 Kxt 



A 



(b) 



5C^' ] Pl 



A 



-\~ai 



Sg 
9 



'-6Z'' 



-5Z 



R 



A 



':5Z\, 



(b) 



(5C" P, 



t-l-h.c. 



(6.1) 



5 A 



(6) 



6Z'^A% , 



^Squark decays have been discussed at the tree-level in several places of the literature. See e.g. Refs. 1 215- 217 1 for some 
relatively recent references on the subject. 
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ha{k) 



tip) 




Xiip') 
Figure 6.1: Tree- level Feynman diagram of the process ba — > ^Xi~ 



SC!: 



A'' 



(b) f^9 ^'^ 



9 2 
where we have introduced the shortcuts 



SZ'' 



SZl 



-5Z% 



5 A 



[b) 

— ai 



ha ^ (b) 
— ai ' 



5Z'"'A 



(6.2) 



6A^'''> 



SA 



(b) 






(6.3) 



We have not introduced the mixing self-energies between the two charginos {5Z^^) as they do not 
contribute. In the case of the QCD corrections the chargino gets no correction, whereas the Yukawa 
coupling approximation implies no mixing between charginos (see section 6.5). For the very same reasons 
we do not show the shift that the chargino mixing matrices U and V would develop at one-loop. 

The full structure of the four on-shell renormalizcd decay amplitudes for ba —> t Xi (a = Ij 2; i = 1,2) 
follows from the previous Lagrangian after including the contributions from the (LH and RH) one-loop 



vertex form factors -Ff*^: 



iTiba^tX- 



igut 



(^^'L+Ar) 



^t ) Pl 



A 



where 



Kf = (5C" + e,; F\ 



i^ L I 



Afl - 5C'\ 



(6) 
-\-ai 



Fi 



^R ) Pr 



(6.4) 



(6.5) 



Let now Fg* be the tree-level partial width of the decay ba — > tXi ^ the only Feynman diagram 
contributing to this process in depicted in figure 6.1, and from the Lagrangian (2.35) we can obtain: 



^0 — 



16 7rm~ 

br 



■\''^a^^,t){[{AfJ^ + {A^^lf] {rnl 
-^A^:iA^^lm,eah}, 



Mf 



(6.6) 



with \{a,i,t) — \{m? , Mf, mf) the usual Kallcn function for the given arguments. The quantum correc- 
tion to Fg' can be described in terms of the quantities 5'" — (F"^* — Fg')/FQ*, where F°* is the corresponding 
one-loop corrected width. From the previous formulae, 5°"^ can be worked out as follows: 



2 (mi - Mf 



f) {a^'I K' + Afa^ K) - 4 ™t 6, AU (aL'L A?,^ + A^+l A20 



'{Afj^ + [A^y] (. 



M? 



m^,)-AAf^^A(^lm,e,M, 



(6.7) 



As mentioned in chapter S we use an on-shell renormalization procedure, so the input parameters for 
the sbottom sector will be the masses as well as the mixing angle 



(m^^,m^^,6l^), 
whereas for the stop sector we just have in addition 



(6.8) 



(6.9) 
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since by SU{2)l gauge invariance the value of the other stop mass mi^ is aheady determined. Of course, 
tan 13 and the SUSY Higgs mixing parameter /i are also additional independent inputs for our calculation. 
Similarly, the sbottom and stop trilinear terms Ab and At are fixed by the previous parameters as follows: 



A = 



/i tan (3 - 



m? — 



bi 



2 mb 



sin 2 0jj ; At = fi cot (3 - 



m? — m 



t2 



'ti 



sin2 6lr, 

2mt *' 



(6.10) 



of course the ^2 niass, as well as the A parameters form eq. ( 6.10| ) will receive radiative corrections, 
but these would be second order correctio ns to our process. We must be careful with the value of 
these A parameters, as explained in section 2.4.3 they are bounded by the condition of colour-breaking 
vacua (2.24), as well as the condition of perturbativity of the Higgs-squark-squark couplings. These 
bounds will translate into a forbidden region in the parameter space defined by (6.8) and (S.9). As the 



conditions mentioned above are fairly qualitative we will present our results in all of the parameter space. 



but at the same time we will single out the regions where the conditions (2.24) are fulfilled 



6.3 Tree-level results 



In this section we make a simple tree-level analysis of the MSSM parameter space, trying to single out 
the regions where the process under study should be interesting. In order to achieve this we focus on the 
sbottom decay itself, but also on the sbottom production mechanism and on the top quark decay. 

As we are interested in the effects that this decay could have at the Tevatron we should use relatively 
light sbottom masses (a few hundred GeV). On the other hand, as this processes would imply a growing 
of the top quark production cross-section (pT^), it is necessary a mechanism that would keep this cross- 
section at its measured value, so the top quark should have available non-standard decay channels, e.g. 
t — ^ iJ+ 6, thus the charged Higgs mass should comply with 



Mf{± < mt — nib 



(6.11) 



another possibility could be the supersymmetric channel t ~^ iiXiJ however this latter channel cannot 



be under control, for example, in the Yukawa approximation performed in section 6.5 it is necessarily 
blocked up. 

It is clear that the radiative corrections to the process ba —<■ t Xi will only be interesting in the region 
where it also has a large tree-level branching ratio. Apart from the already stated gluino decay channel 
there are also other channels {ha — > ^Xai ^2 -^ bih^ , . . .) that will contribute to this decay width. To 
have an appreciable branching ratio ba —^ t Xi we start out supposing that the gluino is much heavier 
than the squarks 



rUg > mi^ 



(a= 1, 2) 



(6.12) 



neutralino masses, on the other hand, are related to chargino ones, so no additional conditions can be 
put on this side. Let us define the branching ratio in which we are interested: 



BRoiba-^tx^) 



^o{ba -^ Xl t) 

Y^Toiba^bxD+Toiba^tXi 

a 

+ J2^o(ba^bb^°), 



Toiba -^ h H-) 



(6.13) 



(where ^^ = h'^, H'^, A'^). To maximize this branching ratio we should work in an scenario where the 
lightest chargino is higgsino-dominated, and tan/3 should have a low- moderate value, if tan/3 is large 



( > 40) then Fq is dominated by the neutral higgsino contribution (first summand of expr. (6.13)) 



In figure 3.2 we have plot the value of the branching ratio (|6.13|) as a function of tan/3, mj^ and Of^ for 



given values of the other parameters. From the figure it is clear that low tan/3 enhances this branching 
ratio. From now on we will concentrate in the region of tan/3 ~ 20, with this typical value the branching 
ratio still has an appreciable branching ratio, whereas the electroweak corrections can be enhanced by 
means of the bottom Yukawa coupling (2.18). In Fig. |6.2| (b) we can see the opening of the Higgs channels, 
namely 62 -^ 61 $1 (at the left end of the figure) and 61 -^ ii H^ (at its right end), it is clear that when 
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Figure 6.2: (a) The branching ratio of ba ^ Xi t ^-s a function of tan/3 for the various decays a 



with TTijj < mj^ ; (b) As in (a), but as a function of mj^ ; (c) As in (a), but as a function of 
marked parts of the abscissa in both figures are excluded by the condition (p. 2 
for (a) and (b) are given in the frame. 



-- 1,2 

The 

The fixed parameters 



this channels are open they tend to take the branching ratio ( 6.13| ) to undetectable value s. T his large 
value of the Higgs decay width is due to the fact that in t his r egions the A parameters ( S.IO ) acquire 
large values. Of course one could fix the input parameters ( |6.8|) in such a way that ^{t.h} are small in 
one of this regions (say at m^^ light), but at the price of making them large at its central value and even 
larger a t the other end. This effect is also seen in Fig. S.2(c), as the A parameters are related to the angle 
trough ( 6.10| ). Another possibility wo uld be to push the Higgs masses to a high value, which would be 



in contradiction with the requirement (3.11). Note that the allowed range of 0^ is rather narrow, so that 



the physical sbottom masses basically coincide with the LH and RH electroweak eigenstates. 



6.4 QCD corrections 

The first step to the computation of the quantum effects on 6a ^ i x~ is to compute the QCD corrections, 
as they are expected to be larger than the EW ones. 

The evaluation of the one-loop QCD corrections to the decay 6q — > i x7 comprise the computation of 
the gluon and the gluino mediated diagrams. These corrections were computed in [p7|,|5q|, however we 
use slightly different renormalization conditions, and we present our numerical results in a completely 
different way. We have checked analytically our results with that of |^,Q, and also numerically in the 
case of Q . 

In this process it is not possible to separate between the gluon-mediated and the gluino-mediated 
contributions, this is so because there are supersymmetric (i.e. i?-odd) particles in the external lines, and 
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(a) 




(b) 



Figure 6.3: One-loop QCD vertex diagrams contributing to the process ba ^ tXi 



thus the supersymmetric theory must be taken as a whole to be renormalizable. 

The one- loop vertex Feynman diagram s co ntributing to the one-loop form factors (|6J|) are depicted 
in figure 6^. The gluon contribution (Fig. p.3| (a)) can be written as 






(6.14) 



where Cp = {N^ — l)/2 Nc — 4/3 is a colour factor and 



Sg2 = 2MimtCi2 , 



(6.15) 



where we have introduced a small gluon mass A to regularize the infrared divergences. The gluino 
contribution from Fig. p^(b) is far more complicated 



Fr = 



—iSnas Cp 



-{ml Cn+Co^{ml ~ Mf)C^2)Ali; R'tr U, 
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^-, TTi r< aW * dC*)* p(*)* 
m^mfc Co A+b, ^2a R2b 

fUb mtiCo + Cn - Ci2)A%* 4a * i?[f 
m,M,iCo + Cn)A^VRf^*R['^* 



"-26 



Fr 



= Fl[a+^a^,r[*J^4*J 

= C*{k,-p,mb,mg,mfJ . 



(6.16) 



The i nfrar ed divergences from (6.15) and from (3.15), (3.71) cancel with the real corrections from the 
diagrams 6.4. The gluon bremsstrahlung contribution to the ba decay width is 



r? 



9^ as 
Gtt^ mr 



{{A^^lr + {A^!:lf) [2m?Jm[ - m? - Mf)Io. + 2 (m[ - Aff - mf) I, 
-2 m1li-2 m1{m^t + Mf - m?^ )/ii + I^ 

+ Ii+ mlhi + {mj - Mf + mljloi + Iq] } 



+2 {mliml -Mf) + M^{Mf - mlj 



SAf^^A^m^M. 
I^{mi^,mt,Mi) , 



(6.17) 
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Figure 6.4: Real QCD corrections to the process ba -^ Xi 
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Figure 6.5: The correction Sqqj^ to the decay wid th b g ~* tx as a function of (a) tan/? and (b) the 
higgsino mass parameter /i. Set of inputs as in Fig. 3.2, 



where we have used the bremsstrahlung functions defined in [9^]p|. We have checked exphcitly (analytically 



and numerically) that after adding up the one-loop (6.7) and the real (6.17) corrections the final result 



^QCD 



is free of ultraviolet and infrared divergences. 



n 






(6.18) 



In Figs. 6.5-3.7 we present the evolution of the corrections (3.18) with some of the parameters. For 



the numerical evaluation we use as{mi ), using the one-loop MSSM /3-function, but, for the m^ we use, 
it is basically the 4-flavour SM /3-functio n, as the scale is almost always below the threshold of coloured 
SUSY particles (and top quark). In Fig. |6.5| we can see the evolution with tan/3 and /i, which are the 
most interesting ones. The corrections are large (> 10%) and present a weak evolution for large values of 
tan/3 ( > 20). We remark that for /i < —120 GeV and tan/3 > 20 the corrections can be very large near 
the phase space limit of the lightest sbottom decay. However, this effect has nothing to do with the phase 
space exhaustion, which is described by the kinematic function X{a,i,t) on the RHS of the tree-level 
expression (3.6), but rather with the presence of the dynamical factor in brackets on that equation which 
also goes to the denominator of 5 in eq.(6.7). That factor is fixed by the structure of the interaction 



Lagrangian of the sbottom decay into charginos and top; and, for the parameters in Fig. 6.5, it turns 
out to vanish near (actually past) the phase space limit in the case of the lightest sbottom (&i) decay. 
However, this is not so either for the heaviest sbottom (62) or for /i > 120 GeV as it is patent in the 
same figure. The different evolution that present the corrections of the two sbottoms has more relation 
with the electroweak nature of the process than with the purely QCD loops, it is due to the fact that. 



We have corrected a typo present in expressions D.ll and D.12 of Ref. 
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Figure 6.6: The correction ^gp^ as a function of d^. Inputs as in Fig. |6.2 
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being 0^ and Qi so small, the squarks are mostly chiral, namely 



6i ~ hR , 62 ~ 6l 



U 



tR , ^2 — tL 



(6.19) 



SO its very different couplings to charginos ( 2.36| ) provides a very different evolution of (6.7), even if Arg 
were constant. In fact the sbottom mixing angle plays a crucial role in this corrections as seen in Fig. 3.6, 
however we also see that its value is highly constrained by the condition ( ^.24| ). Finall y w e would like 
to comment on the effect of the gaugino mass parameter AI and the gluino mass in Fig. |6.7| . The gluino 
evolution is rather flat once the pseudo-thresholds ofba^bg are passed, so even though the gluino could 
not be produced at the Tevatron it would have an effect on the sbottom decay^. As for the gaugino mass 
parameter the correction is saturated for M > 200 GeV, so the corrections computed in this section can 
be compared with the ones obtained in the higgsino approximation discussed in the next section. 

The other parameters of the model present a rather mild effect on the corrections for squark masses 
in the ballpark of several hundreds of GeV. In summary the QCD corrections on the decay ba ^r txj 
are large (~ —20% for 62, — —60% for bi) and negative for values of the parameter space relevant to the 
Tevatron energies, with a higgsino- like chargino and moderate or large tan/3. 



6.5 Yukawa corrections 

At large (> 20) or small (< 1) tan/3 these effects could be competitive with the QCD corrections of 
the previous section. Since in these conditions the full MSSM quantum effects can be rather large, 

^In |57[ it is shown that there exist anpn-decoupling effect at large gluino masses, however this effect is numerically 
small and is not the one reflected in Fig. p.7|(a). 
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(a) (b) 

Figure 6.8: Feynman graphs contributing to the chargino self-energy in the Yukawa approximation. 



their calculation is indispensable to account for the observed top quark production cross-section (1.4) 
in the MSSM or, alternatively, to better assess how much the determination of the SM branching ratio 
BR{t — > W^ b) is affected in the MSSM context after plugging in the experimental number on the LHS 



of eq.(y^. 

The analytical formulation developed so far in chapter |^ and in section 5^ is well suited to tackle 
the general problem of the SUSY-EW corrections to squark decays. Since the dominant part is from the 
Yukawa sector we wish to pursue our calculation in the following within the Yukawa coupling approxi- 
mation. This means that we are going to compute the leading electroweak effects of O(A^) and 0{Xl) 
that emerge for large values of the Yukawa couplings ( |l.2| ) when the remaining gauge contributions - of 
0{g^)~ are subdominant. In practice we shall only explore the large tan/3 regime, typically tan/3 > 20; 
the possibility tan /? < 1 is not so appealing from the theoretical point of view. Thus within our approx- 
imation we will include the correction A,- in leading order 0{X^) of the r Yukawa-coupling, A,-. Notice 
furthermore that for Af, ^ 1 the tree- level decay rate, eq. ( |6.6D , is maximized. Therefore, the large tan /? 
range is expected to be the most relevant one for the decay under consideration. 

In our approach, we set the SU{2)l gaugino mass parameter Af » |^|,Af^ in the chargino mass 
matrix (see section ^.4.4 ), and therefore the chargino Xi is mainly higgsino, whereas the chargino X2 is 
mainly gaugino and does not contribute to our decays. It is only in this case that the Yukawa- coupling 
approximation makes sense. Thus, since m^^ > 80 — 90 GeV, in this approach the decay into stop and 
neutralino t ^ taXa is kinematically forbidden. In this approximation the relevant counterterms SA^j^^ 



in eq.(6.2) boil down to 



SA 
6A 



(b) 
(b) 



?(b) 



SR': Xb - R\: sx, 



'■2a 

Ab) 



(b) 



SRY'X.-RY'SX.. 



2a 
{b} 



.20) 



In the higgsino approximation only two Feynman graphs contribute to the Xi self-energy, and they 
are depicted in Fig. p.q . We can write the chargino Xi self-energies, defined like the fermion self-energies 
in chapter ||, as 

Ei(fc2) = -tg^Nc{A^'l,)'B^{k,mt,m-J 
^Uk') = ~i9'NciAf^,fB,ik,mum,r) 



s^(fc') 



+t g' Nc^e, Af^, A^'ii i?o(fc, m, m~, ) , 



M,: 



(6.21) 



for diagram 6.8(a), note that in this approximation the chargino Xi coupling matrices are simply 



A 



(b) 



y{b) 



~XbRZ' , A 



(b) 
al 






we prefer however to maintain a more general notation. The bottom-stop contribution (diagram 3.8(b)) 
can be obtained performing the following substitutions to (6.21): mt — > mi,, mj^ — > mi , AjJ^^ — > 



A-J^^. From these self-energies we compute the wave function renormalization constants with the help 



of expressi on ( 3.11). A brief comment is mandat ory r especting the relation between the renormalization 
constants ( 6.21 ) and that of the charged Higgs ( [3.55 ). The chargino definition ( 2.32|) , in the case of the 
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Figure 6.9: Feynman graphs contributing to the sbottom self-energy from the Electroweak sector. 



higgsino approximation, is 






(6.22) 



thus in a supersymmetric renormalization we should obtain that the wave function renormalization con- 
stants obtained from (6.21) are 

SZl - SZh, , SZ}^ = SZh, , (6.23) 

as we are not dealing with a supersymmetric renormalization procedure this is not the case, however we 
have checked that the divergent part of this different renormalization constants is the same. This fact is 
crucial in the cancellation of the divergences of SZK ^^i with the term 6Zh± appearing in the definition 
ofJtan/3 (^!5ll). 



From diagrams of figure 6.9 we can obtain the sbottom self-energies. The cliargino top contribution 
from Fig. 3.9(a) is 



KAk') 



A^.^. + A^.AnMBo + k^B^ 



+e. M, mMTM A'^-l + ^^ <'L)SoJ {k, M„ m,) 
and the Higgs contributions from diagrams p^(c) and (d) is simply 



(6.24) 



(6.25) 



with H = {H+, G+; /i", i/", A^, G°),q= {i; h), and Gf^ is the Feynman rule for H -^ q* q' , with gi,2 the 



weak eigenstates squarks, divided by —ig p^UMn- On the other hand diagram 3.£(e) only contributes 



*Note that with this convention G^. = —G^,. 
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Figure 6.10: Feynman diagrams contributing to the one- loop vertex form factors of the decay ba ^ tXi 
in the Yukawa approximation. 



to the mass counterterm and mixing wave function renormalization constants, 



Ef,^ ^^9' Y. ^» <* E?,A,{Mn) + ^ Y. ^ <* E^^AoiM^) , (6.26) 

H+,G+ h°,H°,AO,G° 



with E^^j the corresponding Feynman rule of H H* -^ b'* 6' , 61^2 the weak eigenstate squarks, divided by 



-ig 



Then the diagonal wave function renormalization constant is 



sz^. 



-i2g^ 



2e,M,mtA%A'^^lB', 



\ ({A%f + (A^^lr) ((m? + 5Mf+ mljB', Bo 



(6.27) 
(6.28) 



The neutralino-bottom contribution form Fig. |6.9| (b) can be easily found by performing the following 
transformations in expressions (6.24) and ( 3.27| ): substitute nit -^ mi,, the chargino indices by neutralino 
ones i — > a, and divide the expressions by 2. 

The rest of the renormalization constants are computed using the very same expressions of chapter 
by taking the Yukawa approximation, i.e. by removing the interactions with the gauge bosons (but 
maintaining that of the Goldstone bosons) and with the gauginos. 



The one-loop Feynman diagrams contributing to the vertex form factors can be seen in Fig. 6.1C. 
There exists other possible diagrams, but they do not contribute in the Yukawa approximation. From 
this diagrams we can compute the corresponding -F/l.a} form factors. 
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(6.29) 
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Figure 6.11: (a) The SUSY-EW corrections (6.7) as a function of m^^; (b) As in (a), but as a function 
of mi^ . Rest of inputs as in Fig. 6.2. 



for diagram 3.1C|(a) , and with the help of the following combinations 
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(6.30) 



the contribution from diagram diagram 3.1C(b) is 



-,2 r 



F]' = -i^ 



- mt (Cii - Ci2)(M, ^(3) + ^^ ^(4) + ^^ ^(6) ^ ^jO ^(2)) 

+ M,Ci2(M,A(6)+mU^'^+mtA(6)+M°A(^)) + C'o] , 






(6.31) 



For the numerical analysis, we follow the directions given in section 
segment under consideration, namely 

20 < tan/3 < 40, 



In the relevant large tan/? 
(6.32) 



the bottom quark Yukawa coupling Ah is comparable to the top quark Yukawa coupling. At . Even though 
the extreme interval 40 < tan /3 < 60 can be tolerated by perturbation theor y, wc shall confine ourselves 
to the moderate range ( |6.32|) . This is necessary to preserve the condition ( 2.24 ) for the typi cal set of 
sparticle masses used in our analysis. We point out that the colour stability requirement ( 2.24 ) could be 
satisfied independently of tan (3 if the A-parameters would be chosen directly as a part of the set of inputs 
and then taken sufficiently small. Nevertheless this possibility is not so convenient in our analysis where 
the sparticle masses are the natural inputs that we wish to control in order to make sure that sparticles 



can be produced and decay at the Tevatron as explained in connection to eq.(1.4). 

The corresponding corrections J"' (pj) are shown in Figs. 3.11(a) and 6.11(b) as a function of the 
lightest stop and sbottom masses, respectively. The allowed range for the sbottom and stop mixing angles 



is conditioned by the upper bound on the trilinear couplings and is obtained from eqs.(3.1C) and (2.24). 
In the physical 9^ range, the variation of the correction ( Jf).T\ ) is shown in Fig. |6.12| (a). On the other hand 
the permitted range for the stop mixing angle, 9^, is much larger and we have plotted the corrections 
within the allowed region in Fig. |6. 12(b). Notice that the sign of the quantum effects changes within the 



dom ain of variation of 9^. F inally , we display the evolution of the SUSY-EW effects as a function of tan /3 



(Fig. 6.13 (a)) and of /i (Fig. 6.13 (b)) within the region of compatibility with the constraint ( 2.24 ) 
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Figure 6.12: (a) Evolution of the SUSY-EW corrections as a function of the sbottom mixing angle, 0^, 
within its all owed range; (b) As in (a), but as a function of the stop mixing angle, Qi. Remaining inputs 
are as in Fig. 6.2. 
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Figure 6.13: (a) The SUSY-EW correction as a function of tan/?; (b) As in (a), but as a function of /x. 
Rest of inputs and notation as in Fig. 3.2. 



A few more words are in order to explain the origin of the leading electroweak effects. One could 
expect that they come from the well-known large tan/3 enhancement stemming from the chargino-stop 
corrections to the bottom mass (Cf. Fig.[1.4|(b)). Nonetheless this is only partially true, for in the present 



case the remaining contributions (Cf. Figs. |6.8| , |6.9| and 6.10) can be sizeable enough. One can also think 
on the SUSY counterpart of Fig. 4.4(b), which we have depicted in Fig. 6.14, as an additional leading 
contribution, as, in addition to the tan/3 enhancement, has an Aq enhancement. However the addition 
of these two kind of contributions does not account for the total behaviour in all of the parameter space. 
To be more precise, in the region of the parameter space that we have dwelled upon the bottom mass 
contribution is seen to b e dominant only for the lightest sbottom decay and for the lowest values of 
tan/3 in the range ( 6.32 ). This is indeed the case in Fig. 6.12 (b) where tan/3 = 20 and therefore the 
bottom mass effect modulates the electroweak correction in this process and S^^ becomes essentially an 
odd function of the st op m ixing angle. This fact is easily understood since, as noted above, sbottoms 
are mostly chira l -eg . ( |6.19 )- and the 5_r is the only one with couples with A;, -eq. (2^). On the other 
hand, from Fig. 6.13| (a) it is obvious that the (approximate) linear behaviour on tan/3 expected from 
bottom mass renormalization becomes completely distorted by the rest of the contributions, especially 
in the high tan/3 end. In short, the final electroweak correction cannot be simply ascribed to a single 
renormalization source but to the full Yukawa-coupling combined yield. 

In general the SUSY-EW corrections to V{ba — > txr) ^'"^ smaller than the QCD corrections. The 
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Figure 6.14: Finite Feynman diagram contributing to the b mixed self-energy in the EW basis. 

reason why the electroweak corrections are smaller is in part due to the condition (p^24) restricting our 



analysis within the tan/3 interval (6.32). From Figs. 3.12 and 6.13(a) it is clear that outside this interval 



the SUSY-EW contributions could be much higher and with the same or opposite sign as the QCD 
effects, depending on the choice of the sign of the mixing angles. Moreover, since we have focused our 



analysis to sbottom masses accessible to Tevatron, again the theoretical bound (2.24) severely restricts 
the maximum value of the trilinear couplings and this prevents the electroweak corrections from being 
larger. This cannot be cured by assuming larger values of Mh and/or of fi due to our assumption that 
t — > iJ+ b is operative and because ^ directly controls the value of the (higgsino-like) chargino final state 
in our decay, so that basically we have \fi\ < m^^ — Mh- The restriction cannot be circumvented either 
if we assume larger values of mi , for it has been shown t hat t oo heavy stops are incompatible with the 



CLEO data on 6 ^ 57 both at low and high tan/3 [|60|, |l4§, |l64|-g70[ . We point out that the MSSM 
analysis of 6 -^ s 7 also motivated the sign choice ^ /i < in our numerical calculation |60| . Admittedly, 
the situation with radiative B-decays is still under study and there are many sources of uncertainties that 
deserve further experimental consideration. Still, we have used this information to focus on a limited 
domain of the MSSM parameter space. 

6.6 Conclusions 

In summary, the MSSM corrections to squark decays into charginos and ncutralinos can be significant 
and therefore must be included in any reliable analysis of top quark physics at the Tevatron within 
the MSSM. The main corrections stem from the strongly interacting sector of the theory (i.e. the one 
involving gluons and gluinos), but also non- negligible effects may appear from the electroweak sector 
(characterized by chargino-neutralino exchange) at large (or very small) values of tan/3. Failure of 
including these corrections in future studies of top quark physics at the Tevatron, both in the production 
and decay mechanisms, might seriously hamper the possibility of discovering clear-cut traces of SUSY 
physics from the identification of large non-SM quantum corrections in these processes. As already stated, 
we have mainly concentrated on the impact of these quantum signatures in the physics of the Tevatron, 
but important effects are also expected for experiments aiming at the production and decay of "obese" 
squarks at the LHC. The latter type of squarks could be free of some of the restrictions that have been 
considered for the present calculation. 

The present study has also an impact on the determination of squark parameters at the LC. The 
squark masses used in it are available already for a LC running at a center of mass energy of 800 GeV. 
The large corrections of both, the QCD and the EW sector (in the Yukawa approximation), makes them 
necessary, not only for prospects of precision measurements in the sbottom-chargino-neutralino sectors, 
but also for a reliable first determination of its parameters. 



Chapter 7 

Conclusions 



In this Thesis we have performed a study of some of the possible phenomenological consequences deriving 
from the interactions between the third generation matter supermultiplet and the the super- Higgs boson 
sector of the MSSM at the one-level, with especial emphasis on the implications for the top quark and 
Higgs-boson physic at the Tevatron collider. We have done so in the on-shell renormalization scheme, 
using a physically motivated definition of tan p. Our definition of tan (3 has the virtue of automatically 
incorporating the one-loop radiative corrections to the most plausible signature of the charged Higgs (if 
tan/3 > 2), namely the t Vt channel. Remarkably enough, this definition of tan/3 can also be extended 
to the situation when Mh± > rnt for a wide range of heavy charged Higgs masses to be explored at the 
LHC rather than at the Tevatron ^,^. For it turns out that the branching ratio of the charged Higgs 
into T never becomes negligible in that range. 



The effects of one-loop EW radiative corrections to the unconventional top quark decay mode 
t — > H^h are large in the moderate and specially in the high regime of tan/3, where they can easily 
reach values of 

6Ew{t^H+b)c:i+'iQ% (7.1) 



for negative ^ (and positive At) and a "light" sparticle spectrum (Fig. 4.6), and 



5Ew{t -^H+b)~ +20% (7.2) 

for positive /x (and negative At) and heavy sparticle spectrum (Fig. [4.19| ). In both cases we have sin- 
gled out the domain ^At < of the parameter space, which is the one preferred by the experimental 
data on radiative B-meson decays {b — > 57) . 

The previous results should be compared with the QCD and SUSY-QCD corrections which in 
previous studies ^,|48) were shown to reach values of 5qcd — —60%, 5susy~qcd — +80% and 
SsuSY-QCD — —40% for the same scenarios. In the perturbative regime of the calculation, the 
positive EW corrections attained for /j, > {jJ-At < 0) can be of the same order as they are in 
the ^ < {fJ.At < 0) case. Unfortunately, the EW effects never become huge enough so as to 
prevent the total MSSM correction from being highly negative for /x > {fJ.At < 0) - an unlucky 
fact which unavoidably leads to a severe suppression of the corresponding branching ratio in this 
case. Quite in contrast, in the fi < {fiAt < 0) situation there are indeed regions of the parameter 
space where the positive EW corrections could be perfectly visible; namely, in those places where 
the total QCD correction in the MSSM largely cancels out, e.g. around tan/3 = 30 in Fig.|4.6|(a). 
Negative corrections of the same order can be obtained provided fxAt > (Cf. Fig lTo|(b)). 



From these considerations it follows that, if there exists supersymmetric partners of the standard 
particles at a scale 100 — 500 GeV, the unconventional top quark decay mode t -^ H^ b has a partial 
width that differs significantly from the conventional QCD expectations. In this case the analyses of 
the Tevatron data could exclude a region in the tan /3 — Mjj± plane which is substantially modified 
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as compared to recent analyses from the Tevatron collaborations. Thus e.g. for a charged Higgs 
mass of 110 GeV, and using the conventional QCD corrected value for the decay width T{t -^ H^ b), 
present Tevatron data implies that the excluded region is 

tan/3 > 50. (7.3) 

If, instead, we assume that the charged Higgs belongs to the Higgs sector of the MSSM, then we 
find that the excluded values of tan (3 are 

tan/3> 35 , tan/?> 75 , (7.4) 



• 



for the two scenarios presented above respectively (see Figs 4.23| and 4.24). Remarkably there is a 



range in the parameter space (characterized by /it > 0, fiAt < 0) where no value of the tan/3— Mh± 
plane is excluded at all by present data on the charged Higgs decay of the top quark. Although in 
mSUGRA models in the literature one usually claims At < 0, we already emphasized in chapter ^ 
that this is not necessary the case, especially at large tan (3 and for sufficiently large values of the 
trilinear Soft-SUSY-breaking parameter, Aq, at the unification scale. Thus, also in specific minimal 
SUGRA models, one can have At > and so /i < 0, which is the most attractive possibility since 
one achieves large positive corrections to i -^ i7+ b compatible with BR{b -^ 57). 

The bulk of the EW corrections is given by the finite corrections to the bottom quark mass (or 



the bottom Yukawa coupling), which are proportional to —^lAt -see eq. (4.31). However, it should 
be mentioned that there are cancellations among other sources of significant corrections. This fact 
implies that the rest of the SUSY corrections was not obviously negligible from the very beginning. 



Thus, present CLEO data on the partial decay branching ratio BR{b -^ 57) |174| favours positive 
values of the EW corrections. However, the leading component of the quantum effects is the SUSY- 
QCD contribution which depends on the sign of /x (rather than that of Atjj). Therefore, in the 
end the sign that really matters for this process is that of jjl alone. The best possible situation for 
the charged Higgs decay of the top quark would occur for negative ji, since then the SUSY-QCD 
corrections are positive, and the EW corrections are also positive due to the 6^37 constraint. 
On the other hand, if /i is positive, then the total MSSM correction is negative (in spite of the EW 
component which must stay positive) . In this case the bounds on the (tan /3, Mfj±) space could 
disappear, as we said above. And in these circumstances, as explained in chapter ^ there is an 
alternative scenario with relatively light neutral MSSM Higgs boson which in combination with the 
negative t — *■ H'^ b searches could strongly point towards the SUSY nature of these Higgs bosons. 
As for the one-loop Higgs corrections sector of the MSSM, the over-all correction to the decay under 
consideration is very small due to huge cancellations triggered by the SUSY structure of the Higgs 



potential to this decay except for very low values of tan/3 (Fig. 4.6). Work is currently in progress 
to determine the effects of these large corrections to the top quark and charged Higgs associated 
production at the Tevatron and at the LHC. 



FCNC top quark decays into neutral Higgs particles have been reviewed. The correct inclusion 
of the left-right mixing between squarks implies an enhancement of the partial branching ratio 
BRit — > ch) of two orders of magnitude with respect previous estimates. We have performed a 
separate analysis of the SUSY-EW and SUSY-QCD effects, with different approximations. For the 
SUSY-EW sector we have used the super-CKM basis, whereby FCNC are produced through the 
charged sector, as in the SM. For the SUSY-QCD estimate we have supposed a non-flavour-diagonal 
mass elements in the left-chiral squark matrix. We have applied present bounds from EW precision 
data to these elements. The theoretical upper limits of these contributions to the decay width of 
the top quark are found to be (Figs.[5T5| and I 



(7.5) 
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and the typical values for this ratio are 10"^ and 10"^ - 10"* for the SUSY-EW and SUSY- 
QCD induced FCNC decays respectively. We have found that the SUSY-QCD induced FCNC 
decay widths are at least two orders of magnitude larger than the SUSY-EW ones in most of the 
parameter space, thus making unnecessary the computation of both the interference terms and 
FCNC EW induced effects through non-flavour-diagonal mass terms. If data improves the bounds 
on these mass terms they will become more strict, bringing the SUSY-QCD induced branching ratios 
down to the values of the SUSY-EW ones, and then a complete computation would be needed, if 
there would be any hope at all to see these effects at the small branching ratio predicted by EW 
corrections! The value found for BR{t -^ ch) is not sufficiently large to yield measurable effects at 
the Tevatron or at the LCR. There is, however, a possibility that this decay mode could be measured 
at the LHC. 



• 



If bottom-like squarks are heavy enough they could decay into a top quark and a chargino. This 
could serve as an unexpected source of top quarks at the Tevatron, at the LHC, or at the LC. The 
radiative corrections to the partial decay width b ^ tx~ are large, both the QCD and the EW-like 
in the Yukawa approximation. In the case of the QCD corrections, they are negative in most of the 
MSSM parameter space accessible to Tevatron. These corrections are of the order 

SqcDih^tx^) ^ -60% 

SqcDih^txI) ^ -20% (7.6) 



for a wide range of the parameter space (Fig. |6.5|) . In certain corners of this space, though, they vary 
in a wide range of values. EW corrections can be of both signs. Our renormalization prescription 
uses the mixing angle between squarks as an input parameter. This prescription forces the physical 
region to a narrow range when we require that colour breaking vacua is not generated. Within this 
restricted region the typical corrections vary in the range (Figs. |6T2 , 6.13 ) 



5Ew{bi -^txl) ^ +25% to - 15% 

Sswih -^ txi) ^ +5% to - 5% , (7.7) 

However we must recall that these limits are qualitative. In the edge of such regions we find the 
largest EW contributions. We stress that in this case it is not possible to narrow down the bulk 
of the corrections to just the renormalization of the bottom quark Yukawa coupling. Although it 
is true that for moderate values of the parameters (and for the lightest sbottom decay) the finite 
threshold corrections to the bottom Yukawa coupling yields most of the contribution, as soon as 
we take the parameters away of this central values the total corrections deviate significantly from 
the ones obtained using this only term. More work is presently in progress to generalize this results 
to the full SUSY-EW sector, by incorporating the neutralino-like decays, and also the inclusion of 
gaugino-higgsino mixing. 



Our general conclusion is that the supersymmetric strong and electroweak radiative corrections can be 
very important in the top/bottom-Higgs super-sector of the MSSM. Therefore, it is necessary to account 
for these corrections in the theoretical computation of the high energy physics observables, otherwise 
highly significant information on the potentially underlying SUSY dynamics could be missed. This is 
true, not only for the future experiments at the LHC and the LC, but also for the present Run I data 



(and the Run II data around the corner) at the Fermilab Tevatron collider |221 

^Note Added: See however note on section 5.5 (pg.|76|). 
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Appendix A 

Vertex functions 



In this appendix we briefly collect, for notational convenience, the basic vertex functions frequently 
referred to in the text. The given formulas are exact for arbitrary internal masses and external on-shell 
mome nta. Most of them are an adaptation to the g^i, = {+ — } metric of the standard formulae of 



^oM = / d'^q ,„, _ ^, , (A.l) 

Bo{p,mi,m2) = I d^'q—^ ^.21 r/ ' ^n2 ZT} ^ (^'2) 



Refs. [100, 101, 103|. The basic one-, two- and three-point scalar fmictions are 

1 

[g^ — m 

1 
[g2 -m?] [(q+p)^ -mj] ' 

Co(p, fc,TOi,m2,m3) = / rf"g-p^ 2TT7 — \ — \2 STTT — \ ~m2 2T ' (^•3) 

J [q-' - mf\ [(q + pY - m^\ [(17+^ + ky - m|J 

using the integration measure 

^„~= (4-n)^^ (A.4) 

^ ^ (27r)" ^ ' 

The two and three-point tensor functions needed for our calculation are the following 

[Bo,i3^,S^.]b,mi,m2)= /d"g^^ K^rl^'^d 21' (^-5) 

J [q^ -m(\[{q+py -m^\ 



[Co, C^,Cp^](p,A;, mi, 7712,7713) = 

[g^gA',gA'g^ 
[g2 _ ml] [{q + pf - ml] [{q + p + kf - ml] 



"- 1 r„2 „,2l W„ , „\2 _2l r/„ , „ , U\2 ^21 • VA.Dj 



By Lorentz covariance, they can be decomposed in terms of the above basic scalar functions and the 
external momenta: 

So (p, 7711,7712) = Ao("l2) +'71iBo(p,TOi,m2) , 

5^(^,7711,7712) = p^Bi(p, mi,TO2) , 

B^^(p, 7711,7712) = P,iPi'B2iip,mi,m2) + gf_,^B22ip,mi,m2) , 

C'o(p,fc, 1711,7712,7713) = Bo{k,m2,m-i) +mlCoip,k, nil, 1712, ms) , 

C^{p,k,mi,m2,mz) = P^Cn -I- fc^Ci2 , 

C^^(p, fc, mi, 1712, "I3) = PiiPuC2i + k^k^C22 + {PtiK + k^p^)C23 + 9^1^024 , (A. 7) 

where we have defined the Lorentz invariant functions: 

1 
2^2 



Bi{p,7ni,m2) = ■;^-j[ylo(mi) - ^o(m2) - /i-Bo(p, 7711, 7712)], (A. 



S2l(p, 7711, 7712) = ^7^7 7T [("■- 2)^0 ("I2) -27711^0 (p, 7711,7712) 

Zp'^in — 1) 
- nfiBiip, nil, 1712)], (A. 9) 

B22{p,mi,ni2) = — -T-[Ao(77i2)-f277iiBo(p, mi, m2) + /i-Bi(p, 7711,7112)], (A. 10) 

2(71 — Ij 
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Cii 

C\2 



= Y 



Bo{p + k,mi,m3) - Bo{k,m2,rnz) - fiC'o 
Bo{p,mi,m2) - Bo{p + k, 1711,1713) - /2C0 



C21 
C23 



= Y 



Bi{p + k,mi,m3) + Bo{k,m2,m3) - fiCu - 2C24 
Bi{p,mi,m2) - Bi{p + k,mi,m3) - /2C11 



(A.ll) 
(A.12) 



C2 



-{-pk[Bi{p + k,mi,m3) - Bi{k,m2,rn3) - /1C12] 



V[-Bi(p + fc, mi, mg) - /2C12 - 2C24]} , 



C24 = 



1 



2(n-2) 
the factors /i_2 and the matrix Y , 



[Bo{k,m2,m3) + 2to?Co + fiCn + /2C12] 



(A.13) 

(A.14) 



2,2 2 

p +mi -7712, 



/l 

/2 = fc + 2pfc + ?7l2 — JTlg , 

1 / fc2 — pfc 



The UV divergences for n — > 4 can be parametrized as 

e = fi — A, 
2 

A = - + 7b - ln(47r) 



(A.15) 



(A.16) 



being 7^ the Eulcr constant. In the end one is left with the evahiation of the scalar one-loop functions: 

(A.17) 



/ — 1 \ 777 



Boip,mi,m2) 



Co(p, fc,mi,m2,m3) 



Vl67r2 

—i 



167r-^ 



-xi In ■ 



p 



A + In i-2 - 2 + ln[(a;i - l)(a;2 - 1)] 



Xl 



xi — 1 
f -i \1 1 



2:2 In 



X2 



X2-\ 



Vl67rV 2 pk+p^^ 



with 



a:i,2 =a;i,2(p,mi,m2) = -H ^-^ — ± ^A ' (p , mi,TO2), 



(A.18) 
(A.19) 

(A.20) 



and where E is a bookkeeping device for the following alternate sum of twelve (complex) Spence functions: 



Y. ^ Sp 
— Sp 
+ Sp 



yi 



yi~z{ 
y2 

y2 ~ zi* 

yz - 2^1" 



Sp 



Sp 



— Sp 



yi -4 
y2 - zi 

^3 - 1 

2/3 - zl'' 



Sp 
- Sp 
+ Sp 



2/1 



2/1 -^2 

2/2 

2/2 - z|* 

2/3 



Sp 



Sp 



2/3 - z?* 



— Sp 



2/1-1 

2/1-4 
2/2-1 
2/2 - z^' 
2/3 - 1 



2/3 - Z2" 



The Spence function is defined as 



Sp{z) = - / ^ '-dt, 

Jo '- 



(A.21) 



(A.22) 
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and wc have set, on one hand: 



1,2 
Hi 
1,2 



Xi^2{p + k,m3,mi) , 

Xi^2{k, 1713, 1712); 



and on the other: 



where 



yi = yo + ^, y2 



vo _ _yo _ 1 9 + H 



g = —k"^ + m^ — rri^ , h = —p^ — 2pk — m^ + m\ 
and ^ is a root (always real for external on-shell momenta) of 

p^£^ + 2pki, + fc2 = . 



(A.23) 

(A.24) 
(A.25) 
(A.26) 



Derivatives of some 2-point functions are also needed in the calculation of self-energies, and wc use 
the following notation: 

d_ 

dp^ 



■^-^B^{p,mi,m2) = Bl{p,mi,m2) 



(A.27) 



We can obtain all the derivatives from the basic B'n. 



B'„{p,mi,m2) 



16tt' 



p2 A^/2(p2, m^,?7l2) 



xi(xi - l)ln ) - a;2(.X2 - l)ln 

\ Xi J \ X2 

which has a threshold for \p\ = mi + m2 and a pseudo-threshold for \p\ ~ \mi — m2\ 



(A.28) 



